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1 Introduction 


“I think there is a world market for maybe five computers." 


This famous quote from 1943 is allegedly ascribed to Thomas Watson (11, chairman of 
IBM, and represents — whether or not original - the expert’s view on the business of 
computers at that time. Filling whole floors with thousands of heat generating vacuum 
tubes, while still having less computation power than a today’s smart watch, only 
very special applications for computers were imaginable at that time. 


But already in 1965, the situation had changed completely, leading Gordon Moore, later 
CEO of Intel, to the prediction that the structure density on computer chips will double 
every year [2]. This self-fulfilling prophecy has been driving a tremendous development 
in computer science and technology for more than 50 years Bl. Decreasing structure 
sizes and shorter signal paths allow for increasing the computation frequency, but this 
race stopped in the early 2000s, as the dissipated power density exceeded the one of 
a hot plate [5]. In addition, the miniaturization itself starts slowing down [6] since 
computer industry approaches a size where quantum effects, i.e., the quantization 
of electron charges and tunneling effects, become relevant. When the level of several 
atoms per building block is reached, the miniaturization comes to a hard end and 
other ways to improve the computation power have to be found [7]. 


The demand for more computation power on the other hand is continuing to grow even 
stronger [8] and the number of known unsolved problems will not stop increasing. 
Among such problems are the understanding of fundamental biochemical processes 
like the cycle of photosynthesis or the biological effects of pharmaceutical drugs [9]. 
Including large molecules and the influence of quantum effects [10], these systems 
cannot be simulated classically on conventional supercomputers within reasonable 
time. Richard Feynman therefore proposed a different approach already in 1982: 


“nature isn't classical, dammit, and if you want to make a simulation of nature, you'd better 
make it quantum mechanical" |11 


This idea represents the birth of quantum simulation, where real quantum systems, 
so far inaccessible for a further investigation, are mapped onto controllable quantum 
systems, composed of artificial atoms. These man-made quantum systems can be 
tailored to the specific application, which in turn allows us to gain information 
about the real-world quantum systems. 
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Since then, the field of quantum simulation has been evolving and quantum algorithms 
like Grover’s search algorithm and Shor’s algorithm for prime factoring gave 
additional impetus to the field. They promise a faster and more intelligent way of 
finding the solution by transferring classical problems to the quantum world, where 
superposition states and entanglement can be exploited. Allowing applications beyond 
the original idea of quantum simulations, these algorithms created the research field 
of quantum computation. Here, many conceptual analogies to classical computing 
are drawn, such as quantum bits (qubits) for the smallest information unit and 
quantum gates for calculation operations. These analogies allow us to close the 
circle to the opening quotation of Thomas Watson. Now, instead of vacuum tubes, 
big cryostats and huge efforts to control single qubits are necessary and so far, 
consumer applications are again inconceivable. 


The computation power of current quantum computing devices is still very lim- 
ited, which is well represented by a recently demonstrated programmable quantum 
processor with 53 qubits in a superconducting architecture [14]. It is able to sim- 
ulate an academic problem in 200 seconds, for which a classical computer would 
take at least 2.5 days [15]. This achievement is a unique milestone in quantum 
computing and demonstrates the advance of superconducting quantum systems 
against other physical implementations. 


Already since the first realization of a superconducting qubit in 1999, one of the main 
challenges in the field has been the fight against decoherence [16]. Although the 
coherence times have greatly increased by the invention of different qubit designs 
and the improvement of coherence times is already seen in analogy to Moore's 
law [20], the threshold for useful fault-tolerant quantum computing is not yet reached. 
This means that despite the experimental demonstration of quantum error correction 
schemes [21128], the additional number of qubits and gate operations, which is 
required to perform this error handling, is still too large to be of practical use. 


Fighting decoherence is in particular challenging due to the great intrinsic sensitivity of 
qubits against a coupling to other systems. Since qubits are true quantum systems, the 
qubit state will be projected onto one of two states, as soon as somebody "looks", i.e., 
when the qubit couples to a different, observable system. Especially for macroscopic 
quantum systems, the freedom of designs and the possibility to choose their properties 
in a great parameter space come at the cost of increased parameter fluctuations, 
which results in additional noise and decoherence. One major goal is therefore to 
find the most prominent spurious coupling mechanisms and systems, and follow 
new paths to circumvent their influence. 


On the way to this goal, new systems and effects are found to be relevant and are studied 
in greater detail. Therefore, the vulnerability of a quantum state and the extreme 
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sensitivity of a qubit can also be seen as unique feature, enabling the investigation of 
other systems by means of a very sensitive sensor. Together with other mechanisms 
exploiting quantum effects, the basis for the field of quantum sensing is formed. 


Although many sensing experiments have been performed on other physical systems, 
superconducting qubits are again promising candidates, as they allow for an easy 
coupling due to their large magnetic and electric dipole moments [29]. Through 
techniques like dynamical decoupling or the evaluation of consecutive Ramsey 
experiments [31], they can be used to detect noise in a tunable frequency range from 
the millihertz to the gigahertz regime. It is also possible to discriminate between 
magnetically and electrically induced noise on the same device by changing the bias 
point of phase and flux qubits [32] 33]. Via relaxometry studies, superconducting 
qubits have been used to demonstrate the presence of single atomic fluctuators, 
whose frequency can be tuned by electric fields and mechanical strain, and which 
heavily influence the qubit coherence properties [54137]. These results of quantum 
sensing experiments therefore assist the development of more coherent qubits, 
creating a strong link between the fields of quantum sensing with superconducting 
qubits and quantum computation. 


This link is even more obvious in the development of a quantum computer, where 
many qubits have to be combined onto a single quantum chip, together with a 
vast amount of control and auxiliary structures. This variety of circuit details leads 
to an increased number of coupling partners and therefore loss channels to the 
qubit. Although microwave circuits can be classically simulated by numeric methods, 
these simulations do not include the nonlinearity of the qubits, imperfections in 
the fabrication process and experimental inaccessibilities. A precise experimental 
calibration of the cross-talk between individual components as well as a charac- 
terization of the microwave transmission from the control electronics to each qubit 
individually is therefore of great interest. 


To enable this calibration, a sensing scheme is derived in this thesis by utilizing the 
AC Stark shift of a multi-level superconducting circuit. This is then demonstrated by 
measuring the frequency-dependent transmission of microwave signals from room 
temperature devices to the quantum bit. This sensing scheme is additionally of great 
interest for the development of new hybrid quantum systems, where the harmonic 
resonance of the second system could be detected with the qubit. Such hybrid quantum 
systems have gained much interest in research as they are able to combine the 
rapid processing of superconducting qubits with, e.g., the long coherence times of 
paramagnetic spin ensembles or nanomechanical oscillators for ultra-sensitive 
vibrational sensing [43] 44]. In combination with micromechanical and magnetooptical 
devices, an interface to optical photons can be established [45] 46], which is in particular 
relevant for quantum information processing. Here, it would enable to coherently 
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transfer quantum states between remote quantum computing units and allow for 
a secure quantum communication over long distances [47] 48]. In general, hybrid 
quantum systems can be used to combine the benefits of two independent systems 149]. 


The border between the fields of hybrid quantum systems and quantum sensing is 
weak when a quantum bit is used to study a second quantum system, and the naming 
mainly depends on the research direction, even ifthe physical realization is identical. In 
quantum sensing, typically a known quantum system (i.e., the sensor) is used to study 
another, unknown (quantum) system. For a hybrid quantum system instead, the focus 
often lies on application, where benefits are combined and unique features appear 
[49]. Two prominent examples for a hybrid quantum system are used throughout 
this thesis: First, the qubit forms a hybrid system with a superconducting resonator. 
In this scheme, the resonator frequency exhibits a dispersive shift depending on the 
qubit state, meaning that the qubit state can be inferred by probing the resonator. At 
the same time, the resonator acts as band-pass filter for the qubit, protecting it from 
unwanted microwave radiation and excitation loss. Second, we study a qubit-magnon 
hybrid system, which can be used in future studies to create a link between stationary 
superconducting qubits and flying optical qubits, to enable quantum communication. 


Realizing the coupling between a qubit and the collective excitations in a ferromagnetic 
material would drive the highly active field of magnonics into the quantum 
regime, allowing for a more fundamental understanding of the relevant mechanisms. 
In contrast to electronics and spintronics, where the charge and spin of moving 
electrons are used as information carrier, magnonics relies on the magnetic exchange 
interaction between neighboring spins and hence realizes the transport of information 
without a physical movement of charged particles. Therefore, magnonic devices 
circumvent several problems of modern electronics, like heat dissipation by Ohmic 
losses, as demonstrated by the implementation in ferrimagnetic insulators [52]. With 
interference-based components [53], the footprint of computation elements can be 
significantly decreased compared to standard CMOS} technology, new multi-port gates 
can combine several logical elements into one [54], and frequency multiplexing up to 
the terahertz regime allows for a massive intrinsic parallelization of computation 
tasks [56]. With effects like the distortion-free propagation of magnon solitons 
and the Bose-Einstein condensation at room-temperature [58]|59], interesting physical 
phenomena exist, demanding for further research and investigation. 


Although magnon spin waves have been predicted in 1923 and have been subject 
to research and applications since the middle of the last century [61], most studies 
have been performed with high powers at room temperature. Magnon research at 
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low temperatures and low excitation power has only recently emerged, opening new 
possibilities to approach and study the spin system's ground state [62-65]. A true 
step into the quantum regime can however only be reached when implementing 
a nonlinearity in the system. This would then allow for magnon counting and a 
coherent manipulation of the magnonic state [66169]. 


Again, superconducting qubits suggest themselves for quantum sensing experiments 
with magnons in the low-excitation regime, as they are designable macroscopic 
quantum objects which operate at millikelvin temperatures in the gigahertz regime. 
In this thesis, we therefore place a ferromagnetic material in close vicinity to a 
superconducting qubit, to study their interaction and the properties of magnon 
excitations. However, this implies several experimental challenges, as the magnetic 
material needs an external bias field, while superconductivity is very susceptible 
to magnetic fields. State-of-the-art implementations therefore separate qubit and 
magnetic system spatially by adding a coupling cavity [66]. This only creates a cavity- 
mediated coupling instead of a direct interaction, which adds additional couplings 
and complexity. Therefore, we start to explore the experimentally challenging way 
of a direct qubit-magnon coupling in this thesis. Using a planar approach, the 
magnonic system is placed directly on top of the qubit, meaning that the qubit 
is operated in a magnetic field. 


Since no experimental data on the compatibility of superconducting qubits with 
magnetic fields are available, we investigate this dependency in the framework of 
this thesis, and find that a quantum coherent behavior persists up to magnetic fields 
far beyond what has been expected previously [70]. Besides these promising results, 
the acquired data reveal more insights about the qubit's internal material structure 
and the phase coherence of the qubit state. 


Based on these measurements, different magnetic materials are considered for the 
hybrid circuit and new theoretical descriptions and experimental methods to realize 
the coupling are investigated. Since the magnetic material of choice was not studied 
at low temperatures before, the mode spectrum and internal losses were unknown. 
To this end, the previously introduced sensing scheme can be applied, which enables 
the qubit to detect harmonic signals. 


In the following chapters, I first provide the theoretical background to understand 
the basic function of a superconducting qubit and the expected behavior of a 
superconductor in magnetic fields. After a summary on the measurement techniques 
in the microwave regime, the measurement setup and samples are presented. The 
subsequent discussion of experimentally acquired data is divided into three parts: 
First, the new quantum sensing scheme based on the AC Stark effect is described, which 
exploits the presence of higher excitation levels in the qubit. With that, the sensitivity 
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is increased and a simultaneous measurement of amplitude and frequency is enabled 
[71/72]. The application of this sensing scheme is demonstrated by characterizing 
the frequency-dependent microwave transmission from the measurement lines to the 
location of the qubit. Afterwards, the performance of a qubit-resonator system in a 
magnetic field is characterized and discussed. These data show promising results for 
further applications and reveal the relevance of details in the qubit's fabrication process 
[70]. After a summary of the theoretical background for collective magnetic excitations, 
the coupling scheme between qubit and magnetic material is elaborated and the 
behavior of the magnon-qubit hybrid system is characterized. Here, a rich mode 
spectrum is found, giving a promising outlook on future experiments and applications. 


2 Superconducting quantum bits 


In this chapter, I provide a summary of the theoretical background necessary to under- 
stand the experiments with superconducting quantum bits presented in this work. 
Since many topics and effects have to be discussed and mentioned here to enable the 
reader to understand the experimental data, I tried to simplify the following explana- 
tions as much as possible. Therefore, some statements might look a bit over-simplified 
for the experienced reader, but they retain the key features of our quantum systems. 


This chapter starts with a general theoretical discussion of quantum bits, which is 
valid for all kinds of physical realizations and broadens up to multi-level quantum 
systems. After an introduction to superconductivity, these two topics are merged into 
a section about the superconducting transmon qubit, its features and limitations. 


2.1 Theory of quantum bits 


2.1.1 Qubit states and the Bloch sphere 


In analogy to a classical bit, a quantum bit, or qubit for short, has two possible 
eigenstates, |0) and |1). Due to its quantum nature, it can also be in a superposition 
of these states, resulting in an arbitrary qubit state of 


|) = a |0) +b/1) 


(2) 9 (m Q) 


The complex amplitudes a and b have to fulfill |a|? + |b|? = 1, as |¥) has to be 
normalized with (Y|Y) = 1. The global phase of | V) has however no physical meaning, 
so we can chose 0 < a < 1, a € R without loss of generality. These restrictions are 
intuitively implied by spherical coordinates, with azimuthal angle 0 and polar angle 4: 


0 
|Y) = cos 5 


0 
(0) 4 e? sin 5 |1). (2.1) 
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Figure 2.1: Bloch sphere representation. An arbitrary qubit state |) on the Bloch sphere can be described 
by the two angles 0 and $ known from spherical coordinates. For a pure quantum state, the state is always 
on the surface of the Bloch sphere, ||V)| = 1. The axes are chosen such that (|0) + |1)) 2 is on the x axis 
and (|0) 4- i|1)) V2 points along the y axis. 


As these points lie on the surface of a sphere, this gives the basis for representing 
an arbitrary qubit state in the Bloch sphere, as illustrated in Fig. The common 
choice of aligning |0) along the positive z-axis may seem counterintuitive at first, but 
has historical reasons related to experiments with a spin in a magnetic field. There, 
the bias field is intuitively chosen to be aligned in positive z direction resulting in 
the spin in the ground state also points upwards. Although a great number of qubit 
architectures do not require a biasing field nowadays, this convention persists. 
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2.1.2 Qubit operators 


Any operation on a qubit state can be expressed by a linear combination of the 
Pauli rotation matrices ôy, Oy and 6;. In the Bloch sphere representation, each matrix 
operation corresponds to a rotation around the respective axis by an angle of 7. One 
common combination of rotation operators is the creation/annihilation operator 


Di = 


(ôx +id,), (2.2) 


NIA 


which excites the ground state c4. |0) = |1) and deexcites the excited state o_ |1) = 
|0). The excitation probability of the qubit can be acquired by measuring the 
expectation value of the 0; operator: 


(6;) = (¥|6,|¥) = — |a|? + ||? = —cos8 = —1...1. (2.3) 


For a qubit transition energy of E = fico, this results in a qubit Hamiltonian of 


wol, 
Hg =h 5 Oz. (2.4) 


With the time evolution operator given by U = exp -iH (t — to), the evolution 
of the qubit state is given as 


|¥(t)) = U(t) [F (to = 0)) = cos d |0) + e(ó 716) sin £ 1). 


2 
This means that the qubit state always rotates around the z-axis with a speed of 
Wo1, corresponding to a Lamor precession. For a more intuitive understanding, this 
rotation is typically removed by describing the qubit state in a frame rotating at 
the frequency of the qubit. Since the tool we use to investigate the qubit in our 
experiments, i.e., a microwave drive, rotates itself with a frequency wp, we naturally 
shift our perspective to a rotating frame. Consequently, we only see a remaining 
rotation at the speed of wo — wp. 
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2.1.3 Driving a qubit 


To drive the qubit transition, we apply a microwave tone of frequency wp, amplitude 
A and phase offset 9, which results in the drive Hamiltonian 


Ap/h = Acos (wpt + 9) ox. 
With the unitary transformation 
Û = exp (Sete) , (2.5) 


we transform the drive and the qubit Hamiltonian Eq. to a frame rotating 
with the drive frequency 


Ah = 0 (Ag+ Ap) Ot /h — id 


= whut + 5 CMT +e %o_ perio, + een (2.6) 
^ shuts 4 5 Gia eio.) (2.7) 
= Aus + A (cos(p)öx + sin(q)6y) 
1 A cos ọ Ôx 
= 5 Asing |:| ôy |, (2.8) 
Au 6; 


where the detuning Aw = wq — «p between drive and qubit frequency was intro- 
duced. In a rotating wave approximation, fast rotating terms have been neglected 


from Eq. to Eq. (2.7). 


From Eq. we see that the phase ¢ of the microwave signal sets the axis around 
which we rotate in the xy plane. As this phase can only be defined in relation to 
a second pulse, and the absolute phase of a quantum state |¥) is not of interest, 
the first microwave pulse on the qubit defines the phase reference. We also see that 
the last line containing 0; does not depend on the amplitude of the drive tone. 
This is consistent with the qualitative argument given in the previous section, that 
we only see the frequency mismatch A,, between drive tone and qubit frequency 
rather than the qubit frequency itself. 


For many experiments it is sufficient to consider only a rotation around the y axis, 
resulting in g = 0 and Eq. simplifies to 
A 


" Aw. 
H/h= 59x s mL. (2.9) 
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with the eigenenergies 


125.3. 1 
E,/h = £z / A? 1 42, = ; CR. 


Here, Qe is the generalized Rabi frequency. From a more exact analysis, we get that 
for a non-vanishing detuning Aw, the effective rotation axis is no longer the ôy axis but 
a tilted axis. This leads to an increased rotation frequency Og > A and a decreased 
amplitude, meaning that the state |1) is never reached by a detuned drive [73]. 


2.1.4 Dispersive readout of a qubit 


For different qubit architectures, different readout techniques can be used, where 
e.g. the flux close to a flux qubit is measured or the tunneling current in a 
charge qubit [16]. However, there is a more general approach, similar to the Jaynes- 
Cummings model used in quantum optics. The qubit (an atom) is coupled to a 
resonator (a cavity) and the system is probed by applying photons to the resonator. 
The Hamiltonian for this system is given by 


A A A 


H= Hres Tr Ag pu Hint 


"s Ws REIT 
= hw, (a'a | 5) ; 906; hg (ea 6-4). (2.10) 


Here, w, denotes the resonator frequency with its photon creation and annihilation 
operators ât and A and g is the coupling strength between qubit and resonator. 


Through the mutual coupling, the eigenfrequencies of the system change in dependence 
of the number of photons in the whole system. In the dispersive limit, where 
A, = Wr — wo > g, we can simplify the expression to 


A cx. | " 
Hege/ hz (wr + X0z) at + 5 (wo + X) % (2.11) 


where terms of the order O(g?/A,?) are neglected in a perturbative approach. 
With the dispersive shift x, given for true two-level systems as x = g*/A;, we 
get an effective resonance frequency of the resonator of w, + X0;. Measuring the 
resonance frequency, we perform a projective measurement of 0, i.e., we measure 
the qubit state. Since the wavefunction |V) of the qubit collapses into one of the 
eigenstates of the 0; operator, we get for each measurement one of the two possibilities 
w = wr + x. To measure the probability P, (Y) = |(1|¥)|°, we need to repeat this 
projective measurement and apply statistics. 


The expectation values (ôx) and (ôy), i.e., the projection of the Bloch vector onto 
the xy plane cannot directly be accessed for our qubit types. However, by apply- 
ing a 0/2 or 04/2 rotation pulse directly before the dispersive readout, these 
values also become accessible [77]. 


11 


2 Superconducting quantum bits 


2.1.5 Mixed states 


So far, we only considered pure states with |e) |? = 1 which are located on the 
surface of the Bloch sphere. Under the influence of spontaneous relaxation via the 
ö_ operator or dephasing via random fluctuations of wo, our statistically measured 
state is no longer a pure state. This mixed state can no longer be described by a 
single ket, but by means of a density matrix: 


b= nre) (Hl, (2.12) 


where p; is the fraction of the ensemble state being in state |Y;). These mixed 
states are represented by points within the Bloch sphere, which have a shorter 
distance to the center of the sphere. 


More detailed descriptions of the relaxation and dephasing processes are given in 
Sec. where the measurement procedure is explained. 


2.2 From qubits to qutrits and qudits 


In the theoretical model of a quantum bit, only two qubit eigenstates, namely |0) 
and |1) are considered, disregarding the existence of further energy levels. In the 
vast majority of physical realizations for quantum bits however, more energy levels 
exist and are treated very differently: If these energy levels are far away in terms 
of transition energy, they usually are ignored, e.g. in quantum bits using atomic 
transitions or a spin degree of freedom. In other systems, like nitrogen vacancy (NV) 
centers in diamond, auxiliary energy levels are used to drive conditional transitions 
which favor the population of one of the qubit levels or they are used to read out 
the state of the qubit by photoluminescence [78]. 


Additional energy levels can also be used for special computation or simulation algo- 
rithms, if the transitions between these levels can be addressed individually. Grover's 
algorithm for the search in a large unstructured (quantum) database for instance 
was theoretically modified to operate on multiple energy levels of a single physical 
system instead of many independent qubits. This method was recently demonstrated 
on the four energy levels of a single magnetic molecule's nuclear spin [80]. 


Since the name qubit is reserved for two-level quantum systems, the names qutrit 
for a three level quantum system and qudit for general multi level quantum systems 
have been introduced. However, not every multi-level quantum system can be seen 
as a valid qudit as all levels need to be individually addressable and controllable. 
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This is possible if the system features an anharmonicity, which makes the transition 
frequencies between the levels distinguishable. 


Within certain limitations it is possible to ignore the higher energy levels of a qudit and 
treat it as a true qubit, which simplifies the mathematical expressions and theoretical 
description. To follow common terminology, throughout this work multi-level systems 
will be called “qubits” as long as only the first two energy levels are actively used, 
while higher levels are only considered to explain spurious effects and limitations. 
When more than two levels are required to describe an effect, the term “qudit” is 
employed, even if the same system was previously referred to as qubit. 


Many different physical systems are used to realize qubits or qudits, each having its 
own benefits and drawbacks. For this work, superconducting qubits have been chosen 
as they are macroscopic quantum objects that are well designable and simulatable, and 
offer different coupling mechanisms. Since there are schemes to integrate several qubits 
on one chip and their fabrication methods are closely related to the well-established 
procedures used in the semiconducting industry, they are one of the most promising 
systems to build a scalable quantum computer. So far they are the first and only 
platform, where quantum supremacy, the benefit of quantum computing against 
classical computation, was demonstrated [14]. 


In the next sections, the fundamentals of superconducting qubits are introduced. 


2.3 Superconductivity 


Superconductivity was first discovered by the Dutch physicist Heike Kammerlingh 
Onnes in 1911 [81] and described as a sudden drop in resistivity of pure mercury 
at very low temperatures. Superconductors are therefore ideal conductors and can 
carry a current without any dissipation. Several years later, it was discovered that 
superconductors are not only ideal conductors but also expel any magnetic field 
from their inside as soon as they are cooled below their critical temperature, the 
so-called Meißner-Ochsenfeld effect [82]. Hence, superconductors also are ideal 
diamagnets. These unique properties are however subject to limited validity, resulting 
in a critical temperature Tį, a critical current density je and a critical field He, above 
which superconductivity breaks down. 


Apart from many phenomenological descriptions, no fundamental theory for the 
underlying reasons to superconductivity was found until Bardeen, Cooper and 
Schrieffer developed the BCS theory [83]. The key to this theory is a phonon-mediated 
weak attractive force between two electrons, which leads — despite their Coulomb 
repulsion — to the formation of electron pairs, known as Cooper pairs. Pairing 
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two spin 1/2 particles results in a particle of integer spin, which is therefore of 
bosonic nature. Bosons in turn do not follow Pauli’s exclusion principle and 
the Cooper pairs condensate into the BCS ground state, which is described by a 
single macroscopic wave function 


vr) = |F| ei, (2.13) 


with the macroscopic phase 0(7). This wave function spreads over the whole conden- 
sate, i.e., the whole superconductor. As it does not decay spatially, for any path on 
a closed superconducting loop we can only accumulate a phase of integer multiples 
of 2zt to maintain a consistent wavefunction. In other words the wavelength of the 
Cooper pairs has to be an integer fraction of the ring's circumference. 


The wave vector |k| = 27t/A for a charged quantum particle does not only depend 
on the kinetic energy or momentum of the particle, but also on its surrounding 
vector potential A with V x A = B. We get 


hk = mi 4- qÅ. 


We now demand that the integration of the wave vector around a closed ring 
results in an integer multiple n of 27 


ann= $k di= 7 dö-dr+ t gA ar (2.14) 


and can replace $ A-d7 = f V x A-dS = [B-dS = © using Stokes’ theorem, 
where ® is the magnetic flux. This gives 


we = $ je dF +, (2.15) 
q Ns 


where the supercurrent density js = qns? was introduced together with the density 
of Cooper pairs ns. For a superconductor with thick walls, we can assume that 
the shielding currents are only flowing in a thin layer below the surface and that 
a path with js(?) = 0 exists: in this case the integral vanishes and we finally 
find the flux quantum @®p to be 


d, = 


^ h h 
„E (2.16) 


where q = 2e for the Cooper pairs was used. This means that the flux flowing through a 
hole in a superconductor is quantized and can only be present in integer multiples of dy. 
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If the superconductor however is free of holes, only n = 0 can be the solution to 


Eq. (2.14), as k(?) is continuous. From Eq. (2.15) we hence get 


an, fb oF = —@ =~ [Has 


Using again Stokes’ theorem and comparing the integrand, this results in 


B= ml x fs. (2.17) 
S 


Taking the curl of both sides in Ampére's circuital law VxB= Hojs, where a 
non-magnetic superconductor (u = 0) and j = js is assumed, we get 


wy x p= Vx (Vx B) = V(V.B) - WB = -V38, (2.18) 


where Maxwell’s equation V - B = 0 was used. Combining now Eq. (2.17) and Eq. (2.18) 
gives 


2 
vag = Pm g 
m 
or, regarding a superconducting semispace in x > 0 and vacuum outside 


d?B 2 
Ps) = ar Z Bz(x) 


resulting in the solution 


m 


Hoq?ns (2.19) 


Bz(x) = Bz(0)e At with Ap = y 
This means that the magnetic field does not discontinuously drop to zero but decays 
exponentially from the surface of the superconductor on the length scale of the London 
a da depth AL as we go deeper into the superconductor. This is depicted in 


Fig. Z.Z a). With V x B = pojs the screening currents decay on the same length scale. 


2.3.1 Type I and type II superconductors 


Expelling the magnetic field from the inside of the superconductor costs energy 
and is therefore only possible up to certain fields: Once the energy needed to start 
the screening current and provide the compensation field is higher than the energy 
gained by the condensation of the Cooper pairs, the superconducting state is no longer 
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B.) 


Jy (x) 


Figure 2.2: Boundary effects of superconductors. Magnetic field distribution at the surface of a 
superconductor. a) An infinte superconducting semispace for x > 0 in an outer magnetic field. The field 
decays over a length of AL and screening currents flow at the surface. b) A thin superconductor in a 
parallel field. The magnetic field inside the superconductor never approaches zero, which requires less 
energy for the expulsion of the field. This configuration can therefore be superconducting up to higher 
fields. c) A type II superconductor in a perpendicular magnetic field. Normal conducting cores are 
created, through which flux can flow. They repel each other and align in a hexagonal lattice to have the 
largest distance. 


favorable and superconductivity breaks down. Therefore, a critical magnetic field] 
B. exists, beyond which the superconductor turns normal. 


For different geometries and materials it is possible that not the whole superconductor 
turns normal at the same time, but only parts get normal conductive. In this regime, a 
third energy comes into play: the energy of the interface between superconductor and 
normal conductor. To describe this interface, we use the formalism of the Ginzburg- 
Landau equations and focus on a qualitative explanation of the results. Detailed 
derivations can be found in any standard physics textbooks about superconductivity 
[86488]. To describe the interface, the London penetration depth À; from Eq. 
is relevant together with the Ginzburg-Landau coherence length gL, which can 
be intuitively understood as the spatial extent of a Cooper pair. Due to the strong 
correlation of the Cooper pairs it is clear that any variation in the density of Cooper 
pairs ns can only occur on the length scale of gL. In the normal conducting area 
it has to be ns = 0. In simple words, the superconductivity builds up over the 
length scale cL from the boundary, meaning that the total gain in condensation 
energy is decreased with increasing c and increasing area of the interface. On 


1 The terms “magnetic field" and “magnetic flux density” are used synonymously throughout this thesis, 
following common usage. Whenever a clear differentiation is needed, it will be stated by the used 
quantity symbols, H and B. 
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the other hand, the energy needed to expel the flux the superconductor decreases, 
when Aj is high and the interface is large. 


Eventually, this means that for AL/&cı « 1 the creation of boundaries costs energy 
and is not favorable, and such a superconductor turns normal at once. This is a 
type I superconductor. For the case Àj; /Gar, > 1, boundaries become energetically 
favorable at a certain magnetic field B, ; and magnetic flux starts to enter the type II 
superconductor. To maximize the number of boundaries, the flux inside the normal 
conducting area has to be as small as possible, which means after Eq. that flux 
vortices containing one flux quantum 4 are created, until the second critical field B.» 
is reached, where the superconductivity eventually breaks down. The state between 
B., and B, is called the vortex state and depicted in Fig. [2:2] 9. 


The ratio Àj /&cr can be influenced by the mean free path length ! of conduction 
band electrons in the normal conducting material. While A; slightly increases with 
decreasing l, ČgL decreases strongly. It is therefore easily possible to decrease the mean 
free path length by introducing defects or doping with other metals and with that turn 
atypelinto a type II superconductor. This relation is usually expressed in the form 
that a superconductor in the "dirty limit" is always a type II superconductor. 


2.3.2 Thin films in a magnetic field 


So far, only the influences of magnetic fields on large superconducting samples have 
been discussed. If we now apply a magnetic field in parallel to a thin superconducting 
film (d ~ AL), the field enters from both sides into the superconductor, the screening 
currents cancel out in the center and the superconductor is not at any point free from 
magnetic fields, as shown in Fig. [2.2|b). This behavior can again be described and 
modeled by the Ginzburg-Landau equations and in-depth descriptions can be 
found in original research publications as well as in standard physics textbooks 
[86488]. Due to a large number of degrees of freedom in the experiments presented 
in Ch.|6land Ch.E] where a superconductor in an in-plane magnetic field is studied, 
the derivation of an exact formula here would bring little additional insight. It shall 
hence be sufficient for us to stick to the following qualitative argumentation: 


With the magnetic field entering from both edges of the superconductor, a considerable 
magnetic field is still present at the center of the sample. As already discussed, expelling 
the magnetic field from the inside of the sample is (disregarding other interactions) 
energetically not favorable. This intuitively becomes clear when comparing the 
superconductor to a classical diamagnet where all molecular magnets have to align 
opposite to the outer magnetic field. This means that the magnetic moment per 
unit volume is reduced for a thin film compared to the bulk case and a thin film 
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(i.e. a smaller magnetic needle) is more stable in an outer magnetic field and can 
withstand higher fields, before it “flips”. Eventually, this leads to the experimentally 
proved conclusion that the critical parallel magnetic field for a thin superconducting 
film is much higher than for the bulk material. 


For the case of a thin superconducting film in a perpendicular magnetic field, the 
relevant length scale is not Ay, but again the Ginzburg-Landau coherence length óc; . 
While thicker films with d >> ĉc form characteristic patterns of normal conducting 
areas, the so-called intermediate state, thinner films of type I superconductors, where 
d < Ča, behave just like regular type II superconductors [90] with the creation 
of flux vortices containing a single flux quantum. The qualitative argumentation 
here is similar to the limit of a dirty superconductor where now the thickness of 
the sample restricts the mean free path and coherence length. In the intermediate 
regime d ~ aL however, this approximation is very crude and more sophisticated 
thoughts have to be taken [90]. For our experiments, we can still restrict our thoughts 
to the limit of thin superconductord?] 


2.4 Josephson Effect 


If two superconductors are connected by a weak link, e.g., a thin normal metal 
or insulator, the Cooper pairs can tunnel coherently through this barrier and the 
wavefunctions of the two superconductors overlap. Despite the insulating barrier, 
the current can still flow without resistance, which is called the Jopsephson effect 
[22]. Therefore, no voltage drop occurs, but the superconducting phase changes, i.e., 
there is a phase drop across the junction. The gauge independent phase difference 
between a position before (1) and after (2) the junction, 


27T 2o 
=6,—&- = f Adi, 2.20 
p = 0, -b Oi (2.20) 
is related to the current flowing through the junction via 
I(p) =I. sing, (2.21) 


where k is the critical tunnel current of the junction. From the derivative of the 
gauge invariant phase, we get the equation 


h. dg. 


2 For aluminum, é = 1.6 um was measured (01). 


18 


2.4 Josephson Effect 


meaning that the voltage across the junction is proportional to the change in phase and 
therefore dependent on the change in current. Compared to the voltage at a classical 
inductance V = LI, we can define the Josephson inductance 

Po 1 


eg (2.23) 


A Josephson junction (JJ) can therefore be seen as a nonlinear inductance, where 
the inductance depends on the current flowing. 


This nonlinear element enables us to build an anharmonic quantum oscillator and use 
it as a qubit, which will be introduced in Sec. For the moment, we will have a short 
look on a different application and properties of Josephson junctions. 


2.4.1 The SQUID 


If we interrupt a superconducting loop by two JJs, as depicted in Fig. 2.3]a) it is clear 
that after Eq. 2.15), the flux ® enclosed in the loop can only take integer multiples of 
Po. This means that screening currents in the loop start to circulate when an external 
magnetic flux ®, is applied. If we now drive an external current I, across the junctions, 
their effective critical current is reduced. This setup is commonly called a SQUID - a 
superconducting quantum interference device. It can either be used as very sensitive 
sensors for magnetic fields by increasing the area of the flux loop, or as one effective JJ 
whose critical current I, eff can be tuned by an external magnetic field. 


2.4.2 A Josephson junction in a magnetic field 


When placing a single JJ in an in-plane magnetic field, depicted in Fig. [2.3|b), the 
field enters preferably in the regions of lower critical current, i.e., the junction. This 
also becomes clear from Aj, in Eq. when taking into account that the density 
of Cooper pairs ns is reduced in the JJ. The magnetic field in the junction leads to a 
gradient in the vector potential A, giving a contribution to the gauge invariant phase 
(Eq. (2.20)). We can split the junction into an array of parallel junctions, demonstrated 
in Fig. [2.3]c). Increasing the magnetic field until a total flux of $9 has entered the 
junction area leads to a phase difference between the ends of the JJ of 277, which 
eventually means a complete suppression of any net current across the junction. This 
experiment is comparable to the single slit in optics, where the slit is the junction and 
the magnetic field corresponds to the viewing angle. Continuing this analogy, we know 
from optics that the transmission I(®) is the Fourier transformation of the slit's spatial 
transmission I(x) and follow for the transmission function of the junction I; (d / Po) 
that it is the Fourier transformation of the junction's transmission j.(x). This relation 
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Figure 2.3: SQUID and JJ in a magnetic field. a) A DC SQUID consists typically of two JJs in parallel, 
where the flux ® enclosed in the loop (gray area) induces a screening current Iscreen and hence tunes the 
effective critical current for an additional bias current Ij across the SQUID. b) Cross-section through a 
single JJ with an insulating barrier (I, shaded) between two superconductors (S, white). In an external 
magnetic field, a screening current flows, but the penetration depth A, is larger in the insulating area, 
leading to a non-uniform current distribution. c) We can model a single junction by a parallel circuit of n +1 
individual JJs, forming n SQUID loops. The phase drop across each junction is 9; = $4.4 + 27t + Pn / Po. 
This means for );; 6; = Po we get Pn — Po = 27t and we can find for each junction with 9; a junction with 
Pi + 7, whose currents cancel out each other, leading to a vanishing net current. If the field is further 
increased, this blockade is partially lifted, resulting in a Fraunhofer pattern for I. (B). 


can also be derived in analytically exact ways and was demonstrated experimentally 
[93]. For the experiments presented in later chapters, we usually assume a rectangular 
shaped junction, leading to a sinc shaped critical current of 


o 
[0] sın TE @ 
I, = pP 9| = IP sinc—. 
ICONE pe | rene. 


2.5 Thetransmon qubit 


Combining the findings that a qubit does not exactly have to be a two level system, but 
can also be an anharmonic multi-level system (Sec. [2.2} and that a superconducting 
Josephson junction offers a current-dependent inductor (Sec. (2.4), we intuitively 
find one of the most accessible realizations of superconducting qubits: the transmon 
qubit [17]. Although historically, the transmon qubit was derived as an improvement 
to the Cooper pair box or charge qubit, we follow here the more intuitive way of 
understanding the transmon as a nonlinear quantum oscillator. The most simplified 
electromagnetic oscillator consists of a capacitor and an inductor being connected 
to each other. The parabolic potential of this harmonic oscillator gives equidistantly 
spaced energy levels which cannot be used as a quantum bit. If we now replace the 
conventional inductor by a JJ, we find an anharmonic potential where the resonance 
frequency depends on the current flowing through the junction and therefore on 
the excitation number of the system. 
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For a more mathematical description, we define the degrees of freedom with which 
we want to describe the system. Instead of voltage and current, which are used in 
classical electrical engineering, we choose flux and charge here [94| 95]: 


We assume that the system is at rest in the beginning, #(—œ) = Q(—oo) = 0. 
With Eq. (2.22) we directly see 


o 
= 27° —, 
9 dy 
Expressed in these variables, the energy stored in the capacitor is now 
Q? 
We = = 2.24 
es (2.24) 
and the energy in the Josephson junction 
Do Do 27D 
Wy = | L(IJIdI = — == f . 
y= f iaa = 55 f 1(p)dg = - 73 cos, (2.25) 
This gives the transmon Hamiltonian as 
Q? o 2 
H= I 
mn 
= 4EcN? — E; cos ọ (2.26) 


where we introduced the difference in number of Cooper pairs on the islands 
N = Q/2e and defined the charge energy Ec = e?/2C and the Josephson energy 
Ej = 1.®9/2rt. From the fact that N appears in this equation we see the transmon's 
close relation to the Cooper pair box, meaning that offset charges can still influence 
its eigenenergies. However, when increasing the capacitance, the charging energy Ec 
goes down and with that the contribution of the offset charges. 


To describe the transmon as an anharmonic oscillator, we expand the cosine up 
to fourth order in ọ, resulting in 


a E E 
Happrox = 4Ec Ñ? Ej > p? 2i g^. (2.27) 


Recalling that [9, N] = i, we can compare the quadratic terms with N?, 9? to a standard 
quantum harmonic oscillator, and express the operators as 


p= (2c) (b +5), (2.28) 
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with bt, b being the creation and annihilation operators. This results in the Hamiltonian 


Ecc aie d Ec (p+. 0M 
Happrox = \/BECE) (i | z) Ej — =e (#45) 


x V SEcE; bth — E (8*6 A 1) + const, (2.29) 


where in the lower line, all terms without pairs of b and bt have been neglected, 
as they do not contribute to the transition energies when treating the quartic 
term in perturbation theory [26]. 


From this we easily find the energy of the m-th excited state to be 


Em  ,/8EcEjm < (m? m) - const. (2.30) 


The transmon in this approximation hence has an anharmonicity of &m = Ej 5,1 — 
Em-1,m = —Ec and a fundamental qubit transition frequency fw, = i 8EcE] — Ec. 
Typical values of Ec/h ~ 200 MHz for the samples used in this thesis are enough 
to prevent the excitation of higher levels by the bandwidth of the pulses used 
in the experiments. We can therefore take the transmon as an effective two-level 


system and hence as a qubit. 


We note here that the approximation in Eq. is only valid if 9 is small, meaning 
that we are in the transmon regime with E; > Ec and are only considering the 
lowest excited states. The full solution to the transmon Hamiltonian in Eq. can 
be calculated by using Mathieu functions. However, as it is numerically easier to 
implement the diagonalization of the Hamiltonian instead of using Mathieu functions, 
we only refer here to Ref. for a full discussion of this solution and will use other 
methods to describe the exact eigenstates and eigenenergies. 


With a decreasing Ej / Ec ratio, the influence of N gets stronger which we can probe 
by inducing offset charges by means of a gate voltage N — N — ng. The variation of 
the transition energies Ej; (rg) with respect to the gate charge is given in Fig. [2-4]a)-c), 
where a decreased sensitivity on offset charges with increasing Ej/Ec is apparent. 
A comparison between the exact solution and the approximated transition energies 


from Eq. (2.29) is demonstrated in Fig. D.4| d). 


In the experiments, we can excite higher levels of the transmon by subsequently 
driving the individual transitions, i.e., perform the transition |0) — |1) by applying 
an wg photon, then the |1) — |2) transition by the corresponding w12 photon, and so 
on. Alternatively, we can directly drive the |0) — |m) transition by the simultaneous 
irradiation of m photons of frequency wom / m and the help of m — 1 intermediate virtual 
energy levels. Due to multi-photon processes, these transitions start to appear only at 
higher drive powers compared to the single photon processes. The dispersive shift of 
the readout resonator, i.e., the readout signal, increases for the higher excited states [7]. 
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Figure 2.4: Transmon potential for different Ej/Ec ratios. a)-c) The normalized transmon energies for 
increasing ratios of Ej / Ec. With increasing capacitance and therefore decreasing Ec, the influence of single 
offset charges on the islands gets smaller and the energy levels get flatter. However, the level structure 
gets more harmonic, meaning that Eo» / Ep & 2. d) To resolve these differences better, transition energies 
between adjacent levels are displayed. The solid lines and shaded areas highlight the range of transitions 
when varying the gate charge and therefore their spread shows the susceptibility against charge noise. 
Dashed lines show the approximation of Eq. for comparison. It can be seen that this approximation 
is only applicable for the first transitions and deep in the transmon regime. 


2.5.1 Concentric transmon 


For a transmon qubit, basically only two ingredients are needed: a Josephson junction 
and a comparatively large shunting capacitance. Despite this simplicity, a variety 
of different designs has been proposed in the last decade and many improvements 
could be found. Starting from the first transmon design [17], depicted in Fig. 
a), the capacitance was fabricated as interdigited finger structure, which enables 
a compact design but at the same time increases the electric field strength in the 
gap between the fingers. With this, the coupling to defects on the surface of the 
superconductor is increased, leading to higher losses and reduced coherence times. 
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This disadvantage could be evaded by the introduction of a microstrip transmon [98], 
where the electric fields are drawn into the low-loss substrate by increasing the spacing 
between the capacitor islands and introducing a metalized backside of the chip. Using 
inert TiN instead of the commonly used Al, the existence of surface defects could also 
be reduced. Besides the large size of this design, one drawback is the antenna-like 
structure which increases the qubit's radiative loss as well as spurious coupling of the 
qubit to other elements and signals on the chip. To counteract these disadvantages, the 
dipole moment was reduced by wrapping one electrode around the other, resulting 
in the concentric transmon [99], displayed in Fig. D.5|c) and d). This design offers a 
reduced dipole moment and therefore reduced radiation as well as a confined field 
distribution. For this work, two different designs of the concentric transmon are used: 
one with a single JJ and therefore fixed frequency and one with a gradiometric design, 
featuring 2 JJs. As the gradiometric design is only sensitive to the difference in flux 
between the two inductance loops, it is — in contrast to the majority of other tunable 
superconducting qubits — ideally insensitive to homogeneous magnetic fields. This 
feature is of unique advantage for sensing applications in magnetic fields. 


The concentric transmon was adopted and modified by other groups [100], offering a 
direct access to the qubit by moving the microwave feed lines to the third dimension. 
While this design might be beneficial for scaling up the number of qubits per chip, 
it reduces the flexibility and the ability to place other samples directly above the 
qubit chip, which we will later use in Ch. 


Besides the concentric transmon, other transmon designs exist and are commonly 
used. Here, in particular the "Xmon" design has to be mentioned, which offers 
multiple arms for capacitive coupling and therefore enables error correction schemes 


like the surface code 101]. 


A detailed description of all samples used throughout this thesis and their char- 
acterization can be found in Sec. 


2.6 Loss mechanisms of superconducting qubits 


The excitation of a qubit can decay into a variety of different channels and it is not 
always easy to identify the reasons leading to an experimentally observed relaxation 
time T1. In the following, different decay mechanisms are listed, where the list makes 
no claim to be complete. Again, only qualitative argumentations are given, as the 
experiments later do not allow for a detailed differentiation. When several decay 
mechanisms are present, the resulting decay rate is the sum of the individual decay rates 


s=} i= LT (2.31) 
i i 
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Figure 2.5: History of the concentric transmon design. a) The first transmon design features a SOUID 
loop (gray area) and interdigited finger capacitances [17]. b) Microstrip transmon with capacitive islands 
induces mirror charges in the ground plane (dashed frame) and concentrates the electric field lines into the 
substrate [98]. c) Concentric transmon with low dipole moment in the non-tunable version. d) Concentric 
transmon in gradiometric design. Only the flux difference 6; — ®ı creates a current across the junctions 
and biases the qubit 199]. In a) and d), the lower graphs show an effective circuit diagram of the design 
with the JJs in red and flux-sensitive areas in gray. Figure adapted from [77]. 


2.6.1 Purcell and radiation loss 


When the qubit is coupled to a harmonic resonator, we see from the gö_b* term 
in the Jaynes-Cummings Hamiltonian Eq. that the excitation can swap from 
the qubit to the resonator, where it can eventually decay, dissipate or leak into 
a third system. It is clear that besides the coupling strength g also the detuning 
A = wq — WR between the two systems and the bandwidth x of the resonator influence 
the losses. For the bandwitdh, the full width at half maximum of a Lorentzian-shaped 
resonator is assumed. With this, we get [102]: 


2 x/2 
(x/2)? +42 


I1 purcell = 2g (2.32) 


The Purcell loss can be estimated already at the design stage and is therefore usually 
optimized with respect to a good readout signal, given by the dispersive shift x = g?/A. 
The resonator bandwidth is limited by fabrication methods and materials. 


This single-mode Purcell decay is however not only present for the designed readout 
resonator, but for every spurious mode, the qubit can couple to, e.g., modes of the 
sample box or other modes on the chip. Apart from that, the qubit excitation can also 
relax into non-resonant systems like transmission lines or other control lines by the 
multi-mode Purcell effect. Then, the Purcell loss is described by the impedance of the 
environment and the coupling capacity or inductance between qubit and environment 
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[102]. Due to the enhanced spurious coupling of a qubit with large dipole moment, 
these losses are commonly referred to as radiation losses. 


To avoid these losses, one can either firstly design the environment to give little 
possibilities for spurious coupling, which eventually leads to placing the qubit 
inside a three dimensional cavity. Together with a very low x of these polished 
superconducting cavities, this is the basis for the great coherence times of the 3D 
transmon [103]. Secondly, one could reduce the dipole moment of the qubit, leading 
to a very confined electric field distribution and therefore a decreased coupling to 
other structures. This approach was chosen for this thesis, as, unlike in the first 
approach, it is here possible to place other elements like wide-band transmission 
lines and lossy dielectrics close to the qubit. 


2.6.2 Dielectric and TLS loss 


As already mentioned, one goal for the design optimizations of the transmon 
was to concentrate the electric fields into the substrate, as the interfaces between 
superconductor and substrate can be controlled better than the surface of the 
superconductor to air, and the substrate itself is very pure. The contribution of 
the chip's substrate can usually be neglected: For intrinsic silicon, the loss tangent 
tan ô = Im(e) /Re(e) is below ô < 1077, where e is the complex dielectric constant of 
the material. This would limit the qubit coherence to I; = dwo1 < (0.3 ms)! [104 
105], and is not a limiting factor for state-of-the-art transmon qubits. 


Since for the JJs of superconducting qubits the most common choice of material 
is Al/AlOx/AI, also the qubit's capacitance pads are often fabricated from Al for 
simplicity. One major problem is that aluminum oxidizes quickly at air or in wet 
processing steps, leading to the formation of amorphous aluminum oxide on the 
surface of the superconductor. In this amorphous material, mechanical degrees of 
freedom with an electric dipole moment exist, which hence couple to electric fields, 
absorb energy and increase the dielectric losses in the material. These so-called two 
level fluctuators (TLF) or in general two level systems (TLS) have been studied in 


great detail [34 and are still subject to active research. 


Despite the efforts taken to reduce the field strength in the gap of the capacitor 
and concentrate the field lines into the substrate, the loss tangent of Ajo, = 3: 1073 
for amorphous aluminum oxide is so large that still the small portion 
of the electric field energy in the surface and interface oxides has to be taken into 
account. A more detailed study on an all-aluminum transmon concludes that 
the contribution of the surface oxides is on the order of magnitude of current qubit 
coherence times and hence influences the measured T; values. 
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Although it is possible to fabricate the capacitance pads from a more inert material, 
the insulating barrier is exclusively fabricated from AIO; and commonly, at least a 
short connection to the junctions is made from Al. Apart from that, also residues of 
processing materials, like resist or glue, or even water and bound oxygen molecules 
can lead to dielectric losses on the surface of the superconductor - independent of 
the material in use. Observing fluctuations in the qubit's coherence times, we can 
conclude that atomic fluctuators on the surface of the superconductor contribute 
considerably to the overall losses [31]. 


2.6.3 Quasiparticles 


Quasiparticles (OPs) are excitations of the superconductor, e.g., split pairs, where 
one partner of the Cooper pair is missing. They can be described to be in the range 
from hole-like to electron-like, do not participate in the many-body wavefunction and 
therefore increase the energy of the system [86]. Hence, QPs cannot tunnel through the 
JJ without resistance and a QP current is in principal dissipative [93]. We note at this 
point that this explanation is very simplified, and in a more sophisticated analysis, 
the phase 9 across the JJ has to be taken into account. In the transmon operation 
regime however, p is small and we can simplify our analysis. 


For a superconducting transmon qubit, the influence of QPs on the qubit's relaxation 


rate was in detail derived [108] and further calculated [109 to be: 
x 
Iop Rd —v 2wgAo/h => Tgp xX XQP- 


The QP concentration xop is here normalized to the number of Cooper pairs and 
Ag is the energy gap of the superconductor. 


At finite temperatures, QPs are thermally excited, resulting in a QP concentration of 
XoPtherm(T) = V27kgT/Ao exp (—Ao/kgT) [108]. For a temperature of T ~ 20mK 
and a thickness-dependent gap of thin aluminum films Ao > 200 peV [111], we would 
expect Xoptherm(T) < 10792, resulting in coherence times of thousands of years. In real- 
ity however, QP concentrations of xop ~ 1077 are observed [103], giving 1/Iop ~ 1ms. 


The origin of these excess QPs is still subject to research, and commonly, the following 
factors are considered: With the decreasing thermal conductivity at low temperatures, 
it becomes challenging to thermalize the qubit to the cryogenic environment and the 
chip temperature is usually assumed to be above the temperature of the base stage. 


Also, any kind of radiation with hw > Ag absorbed by the superconductor will break 
Cooper pairs and eventually create a cascade of OPs. To this end, filters for infrared 
light are typically included in the microwave wiring to the sample for high-coherence 
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measurements. Impacts of high energy particles like cosmic radiation have also 
been observed to increase the QP density [112]. 


Apart from that, OPs can be excited by driving a strong current I > I, across the 
junction or applying a voltage. This for instance happens when strongly populating 
the cavity coupled to the qubit, which induces an oscillating voltage across the JJ, 
splitting Cooper pairs and hence creating QPs [114]. 


Eventually, also oe gne fields lead the creation of QPs, where for parallel fields 
it holds xop ~ Bj [1 and the QP recombination rate, i.e., the relaxation rate of 
out-of-equilibrium QPs goes back with increasing magnetic fields [116]. 


Apart from shielding the samples against unwanted radiation and operating at the 
correct power levels, QP traps have been introduced, where a combination of materials 
with different gap energies are used such that the gap of the trap is lower than the 
one of the device, Aop > Aor [117] [118]. From the fact that any excitation of the 
condensate needs to be above the condensation energy, it follows that Egp > Ap. 
Coming from material 1 to 2, the OP will relax, e.g. by interaction with phonons, until 
Ao,p > Egp > Aor. Due to its now decreased energy, the QP cannot go back to material 
1 and is trapped. With the reduced gap however, it is easier to excite QPs in the trap, 
for which reason the traps should be located with appropriate spacing from the JJ. 


In the extreme case of Ao = 0, a normal metal is used as trapping material [1191121]. 
This can also be achieved by suppressing superconductivity locally, e.g. by normal 
conducting flux vortices in the superconductor [113/122] [123]. It has been demonstrated 
that the coherence of a QP limited fluxonium qubit can be enhanced by applying a small 
field, which in turn creates flux vortices at intended locations far from the junctions. 


2.6.4 Movement of flux vortices 


For type II superconductors, it was experimentally observed that their RF dissipation 
increases abruptly when Bc is reached [124], leading to the conclusion that flux 
vortices contribute to the losses. In a gedankenexperiment we apply B; > B.ı to a 
type II superconductor to create the vortex state. Due to the mutual repulsion of 
the vortices, they arrange in a regular hexagonal lattice in the equilibrium state to 
have the largest possible spacing. If we now apply a current along the x direction, 
we immediately get from Ampere's circuital law V x B = poj; 


d . 
ay = Hox, 
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which means that the flux density increases in the y direction, leading to a gradient in 
the density of flux vortices. This means the vortices get pushed to positive values of y, 
or in other words, a Lorentz force acts on the vortices and they start moving. 


As described in Eq. (2.16), on a contour around the flux vortex, the phase changes by 
27. Moving flux vortices therefore lead to a change in the phase of the surrounding 
superconductor and therefore via Eq. to a voltage. The electric field correspond- 
ing to this internal voltage is orthogonal to the movement of the flux vortices and 
orthogonal to the magnetic field, hence it points along the direction of j. With E || f 
however, dissipative processes take place and heat is generated. 


We can come to the same conclusion from a microsopic analysis: Looking at a moving 
flux vortex, we see that it destroys Cooper pairs in front of it, as the inner of the 
vortex is normal conductive and QPs recombine behind the vortex, resulting in 
Cooper pairs. For the non-stationary case, we remember that the flux vortices move 
towards areas with higher B, meaning that the magnetic field is higher in the front 
of the vortex compared to its back. This means that splitting the Cooper pairs in 
the higher field takes less energy as will bee released by the condensation in the 
lower field area, leading again to dissipation. 


As a simple conclusion we find that moving flux vortices contribute to overall losses 
for RF applications [124], which was experimentally proved several times [127-129]. 


To reduce this dissipative influence of the vortex motion, a very common method 
is the creation of pinning centers, where the flux vortices preferably locate. These 
pinning centers can be created by structural defects or the doping of other materials, 
but too small pinning centers still allow the finite-sized vortex to wiggle around, 
leaving a contribution to the overall losses. Alternatively, for thin film superconducting 
structures, flux trap holes are structured into the design of the chip, where their 
size is smaller than the structures they should protect and corresponds to one flux 


quantum at the expected field strength [129-131]. 


When flux vortices are present in the superconductor and the field is ramped 
down again, flux vortices are not directly destroyed. Instead, vortices with opposing 
orientation start to enter the superconductor from the sides and only in an interface 
region between up and down vortices, the annihilate. This description was introduced 
by C. Bean and is hence called the Bean model [132133]. The mesoscopic magnetization 
ofthe superconductor is hence hysteretic and the configuration of flux vortices depends 
on the magnetic history of the sample. To decrease the amount of flux vortices in the 
superconductor, it is therefore not only sufficient to decrease the magnetic field or 
apply a counter field, but a demagnetization sequence has to be applied, where the 
magnetic field is slowly ramped between positive and negative values with decreasing 
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amplitude [131]. Due to the multiple slow ramping processes, it is usually faster to 
switch off the field, warm up the sample above T. and cool it back down in zero-field. 


A good overview over magnetization effects and hysteresis of superconductors 


is also provided by Ref. [134]. 
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Like a large number of quantum systems, superconducting quantum bits have 
excitation energies in the range of several eV to few meV. The transition frequen- 
cies of interest for coupling and sensing experiments are hence in the microwave 
regime. This means we can use commercially available microwave technology 
for our measurements and are able to describe the systems with well-established 
formalisms of microwave theory. 


In this chapter I shortly recapitulate the theoretical description of microwave systems 
with a focus on the means and measures we will later use to characterize our systems. 
The chapter starts with general microwave networks and leads via a single resonator 
to two coupled oscillating systems which form hybridized modes. 


3.1 Scattering parameters for a microwave network 


To specify a network of microwave components with N ports, one of the most intuitive 
descriptions from an experimental point of view is the description by an N x N 
scattering matrix. The elements of this matrix correspond to the ratio of measured 
voltage at output port i to input voltage at port j: 

out 

Sj = aw with V2; = 0, 

] 
where the restriction means that we only drive the system at port j. For V}, 
Vout the complex voltage amplitudes of the oscillating signals are taken, where 
the angle between the signals, arg(S;;) = arg(V?"' IVE); is the phase difference 
between incoming and outgoing signal. Writing the incoming and outgoing am- 
plitudes as vectors, we can write 


jy out _ êy} in 
V =SV", 


or for the most commonly used case of a two-port network 


yout $4 $12 yon 
= à 2 1 
Gs a $25 ym (3 ) 


31 


3 Microwave theory 


If the network under test does not contain any systems that can break time reversal 
symmetry, e.g., ferrites or active components, the network is called reciprocal, resulting 
in a symmetric scattering matrix $T = S. This means that exchanging ports 1 and 
2 does not alter the behavior of the system [135]. 


Typically, the parameters of the scattering matrix have a characteristic dependence on 
the frequency. For active components and in the high power regime, the $ parameters 
also become power dependent, resulting in $(w, P^). The parameters of this 2 x 2 
matrix can be experimentally acquired with a two port vector network analyzer (VNA) 
which allows for frequency sweeps at variable input powers. If the network has more 
than two ports, the full matrix can be measured by terminating all other ports with a 
matched load to avoid reflections and repeating this procedure for all combinations. 


Due to the high dynamic range of the parameters, they are typically given in 
logarithmic units, standardized as decibel. One Bel (i.e., ten decibel) corresponds 
to a ratio of ten in powers, or in general 


P 


re 


Two points in this definition are commonly disregarded, which can lead to errors, 
misunderstandings or simply imprecise phrasing: First, it is important to note that the 
unit Bel always represents a relation of powers, no matter which size of quantity is 
actually displayed. So if we want to refer to a voltage ratio in terms of dB, we have to 
consider the power at the characteristic impedance Z given by P = V?/Z and hence 


P V? V 
Lp = 101 — = — = 
p = Wlogyo Ba 10 logy) vi 201ogio Ver 


which means that a prefactor of 20 appears instead of 10 if we relate voltages, 
currents, amplitudes, and other quantities. 


Second, the Bel always requires a reference value P4. For the $ parameters, this 
is obviously the input power, but it is also possible to take a fixed reference 
value to state absolute power values. For the power levels used throughout this 
thesis, a very common unit is decibel-miliwatt (dBm), where the reference value 
is Prep = 1mW. A power level of Lp = -30dBm then corresponds to an absolute 
power value of P = 10'?/mWw = 1pW. 


Using logarithmic units becomes especially handy when handling attenuation and 
amplification factors of microwave components, as they can be simply added to the 
output power of the microwave source in dBm, as long as all factors are given in dB. 
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Figure 3.1: Schematic of a resonator. A transmission line with impedance Zo is capacitively coupled to a 
resonator with impedance Z4. The interaction or energy exchange rate between the systems is given by Ke. 
Inside the resonator, losses to other systems and dissipation can occur, represented by a rate xj. 


3.2 Microwave resonators 


When the characteristic impedance of a microwave transmission line Zp changes at one 
point abruptly to a value of Z4, a fraction of the wave is reflected at this point, with [135 


| 24— Z9 


Su = ss 
" Z1 + Zo 


(3.2) 


where we note that for the case Z4 < Zo a phase difference of 7t can be observed. This 
scenario is depicted in Fig.B-1] If we now introduce a second change in impedance at a 
distance | form the first change, the wave is reflected at both impedance mismatches 
and we can observe standing waves at a frequency of 

Wn Up n C 


254. An ALLER 


for mode number n, where we assumed that we have a A/2 resonator, i.e., the 


impedance outside the resonator is on both sides higher (or on both sides lower) than 
the impedance inside the resonator. Here we used the phase velocity vp = c/ \/€eg¢ in 
the medium, depending on the vacuum speed of light c and the effective dielectric 
constant in the medium €g¢. The mode-number dependent wavelength is given 
by An = 2l/n, with n € Nyo. 


A resonator can also be created by connecting a capacitor C with an inductance L, 
where a resistor R is typically taken into account to sum up the losses. If the resistive 
losses are low, the oscillation frequency is given by wọ = 1/ VLC. Note that this 
resonator ideally has no higher harmonics, whereas the transmission line resonator 
has an infinite number of modes wy, = nw . The following discussions however apply 
to both types of resonators and L, C, and R denote the corresponding inductance, 
capacitance and resistivity of an equivalent lumped-element circuit [135]. 
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3.3 Quality factors 


Besides the resonance frequency of a harmonic oscillator, an important quantity is its 
quality factor Q, describing the ratio between the total energy stored in the oscillator 
Etot and the energy dissipated per cycle Eqiss = Plosg/ wo: 


where x is the full width at half maximum (FWHM) of a Lorentzian shaped frequency 
response and corresponds to the inverse photon lifetime [135]. In other words, 
the quality factor describes 7t times the number of oscillations the resonator can 
perform, until its amplitude has dropped to 1/e. 


Photons in a resonator can be lost either by coupling into a different system, i.e., 
the transmission line, or by internal losses due to ohmic dissipation and dielectric 
losses. The quality factor is therefore divided into a coupling and an internal 
quality factor, Oc and Qr: 


ld 
QL Qi Q 
Or KL = Kj + Kc, 


where the loaded quality factor Qr, accounts for the total losses. 


While Qe in our case is mainly given by geometric values and therefore fixed for a 
given design, the internal quality factor for a superconducting resonator depends on a 
more sophisticated set of variables. Amongst others, these are the choice of materials, 
details of the fabrication process and the interface cleanliness. It can also be altered 
in an experiment, e.g., by applying a magnetic field. Here, the same mechanisms as 
previously discussed for the superconducting qubit apply, in particular losses due to 
quasiparticle excitations (Sec. and flux vortex movement (Sec. 2.6.4). 


3.4 Reflection spectrum of a resonator 


In the experiments, we observe a transmission line resonator whose one end is 
capacitively coupled to a transmission line of impedance Zo, while the other end is open. 
If we assume that the losses are generated by the oscillating current in the resonator, 
we can assume the form of a serial RLC circuit, which has an input impedance of 
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close to its resonance frequency. Plugging this impedance into the scattering pa- 
rameter from Eq. (3.2), we get 


_ Z-Z 24 2QL/ Qe 
Z1 + Zo 1+ 2iQ, 42" 


S11 (3.3) 


where we introduced the detuning Aw = w — wọ . From this formula we see that 
for large detuning, i.e., in the off-resonant case, we always are close to SOR = 1. At 
resonance however with Aw = 0, we get SR = 1 — 2QL/ Qe, being again real. If we 
shift the result of Eq. by the mean of S981 and SỌR1, we get 


& OR , cR Qr 2 . 6 
S11 (S5 + Sn) / S11 Qe ( 62 +1 152 + z) 


with ô = 2Q;Aw/wo. This gives 


[S11] = QL/Qe 
arg($1) = arctang(—26,1 — 6?) = 2arctan a, 
where arctan2(y,x) is the function providing the complex argument to z = x + iy. 
From this result we see that Eq. represents a full circle in the complex plane, 
centered around x = 1 — Q;,/Qc. For the case Qj, > Q./2, and therefore Q; > Qe, the 
origin of the complex plane is inside the circle and the phase of 51; rolls off over 271. 
Since more energy is lost by the coupling than by internal dissipation, this is called the 
over-coupled case. For the under-coupled case, with Qi < Qe, most energy is dissipated 
internally, the origin is outside the circle, and the phase only changes by less than 7t. 
The case of Q; — Qc is called critical coupling. This behavior is demonstrated in Fig.B.2] 


We also see that we can access the ratio Q1,/Q., and with that the individual values 
of Q; and Q. when fitting a circle to the complex scattering data. To this end, we 
utilize the circle fit routine [136], where more parameters are introduced to correct 
for imperfections. Now, Eq. reads 


S11 = Aet e-t (1 — 2017 Qe dip) (3.4) 
1+ 2iQ, e 


where A is the maximum amplitude of the signal and accounts for overall losses 
and « corresponds to an overall phase offset. t denotes the time the signal takes 
to travel through the setup. This means vpT is the effective length of all wiring, 
which gives a frequency dependent phase offset. Finally, p accounts for impedance 
mismatches and gives a rotation of the circle around the off-resonance point, which 
can be seen by a skewed shape of the amplitude signal. Measured data fitted by 
this equation is presented in Fig. 


35 


3 Microwave theory 


m 

B 

~ -101 — a=5 

a — O/Q=1 
7154 _ 0/0 205 


arg(Si1) (rad) 


r T 
—1.0 —0.5 0.0 0.5 1.0 
Re(S;:) frequency Aw/wo 


Figure 3.2: Visualization of $4; for different Q;. In the complex plane (left) all calculated traces for 
S11(Aw/ wo) are on a circle starting at the off-resonance point 51; = +1. For an internal quality factor 
much higher than Q., the circle has a radius of 1 and no amplitude signal is visible (top, dashed blue line). 
In the phase signal (bottom), a 27 rolloff is visible. The greatest change in the amplitude signal is visible 
for a critically coupled resonator with Q; = Q. (green line). Measurement data can be found in Fig. 5.3] 


3.5 Average photon number 


For the observation of quantum phenomena and for the investigation of power 
dependent effects, it is necessary to know the average number of photons (n) in the 
resonator. This is related to the energy stored in the resonator by Etot = (11) hwo, where 
we can use the definition of the quality factor in Eq. to obtain 


(n) hw = Pos Z. (3.5) 


In the steady state, the dissipated power equals the difference between incoming and 
outgoing power Poss = Pin — Pout = Pin(1 — |S11 2). Since the losses due to coupling 
are accounted for in Pout, the quality factor in Eq. only has to cover internal losses. 
At resonance, we can take SR = 1 — 201/Q.,, finally leading to 


E QiPin (4 ~ Suð) _ 4QiPn [ QL Qr 4 Qr 


m hw? hw; VQ. Qi 


(3.6) 


For constant quality factors, the mean photon number in the cavity is hence proportional 
to the applied power (n) x Pin, when excited on resonance. 
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Figure 3.3: Measured data of over- and under-coupled resonators. A resonator with two different internal 
quality factors caused by different surrounding magnetic fields. The amplitude data are normalized prior 
to the fit to account for the transmission through the microwave setup. The solid lines represent the 
results of the circle fit, where the dashed lines show the “axis” of the circle, i.e., the connection between 
off resonance point and resonance. The measurement with high Q; (blue) shows a well visible tilt of this 
axis, indicated by the angle $ in Eq. (3-4), caused by impedance mismatches. This can also be seen in the 
[51 |? data as an asymmetric curve shape, where the right shoulder is higher than the left one. The quality 
factors extracted from the fit are Qi = 17300, Q. = 8200 for the resonator shown in blue and Q; = 2300, 
Qe = 8100 for the orange one. 


3.6 Coupling resonant systems 


If we not only take one microwave resonator but use two oscillating systems with 
a;(t) = a;e! and let them exchange energy at a rate g, their oscillation modes will 
hybridize when their resonance frequencies are close, |w1 — w2| ~ g, leading to new 
resonance frequencies. We can calculate these eigenfrequencies by solving 


C — 04A g ay —0 
g w— 05] V2 ` 
This gives two solutions for w: 
Q1 +w 1 
— s 2 + 5 V (an w2)? + 492, (3.7) 


which is commonly called avoided level crossing or anticrossing as the resonance 
frequencies “repel” each other. The two eigenfrequencies never approach closer than 
Q4. —w_ > 22, where the equal sign holds for the resonant case, w1 = w2. To resolve the 
anticrossing experimentally, we need an external tuning parameter which changes the 
difference between the two bare eigenfrequencies. This behavior is depicted in Fig. B.4] 
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Figure 3.4: Avoided level crossings with constant coupling strength. Solution of Eq. for a constant 
coupling factor 2g (red arrow) and different behavior of the bare eigenfrequencies (gray dashed lines) to 
the tuning parameter. We see that from the experimentally measured data we can fit and calculate back to 
the bare frequencies, the coupling and the dependence of the bare frequencies on the tuning parameter. 


While Eq. only gives the resonance frequencies, we do not obtain an expression 
for the linewidths of the resonances depending on the frequency. For a full expression 
of the $1; parameters, we could use the input-output formalism [137], which would 
require a more in-depth introduction here. Instead, we follow a more argumentative 
deduction here, based on the previous findings, which comes to the same result. 


For the following argumentation, it is useful to rewrite the reflection coefficient 
of a single resonator from Eq. using the linewidth x instead of the quality 
factors. We hence find 
1 2Ke 

KL + 21Awı : 


Su = (3.8) 


When we now add a second system, we have to account for the losses to this system, 
represented by x2, by increasing the loaded linewidth £j, = xy, + Xo. The amplitude 
not reflected by the second resonator is given by 1 — S22 = 4g/ (xi2 + 2iAw2), where 
Ki» is the internal linewidth of the second resonator. The additional factor of 2 comes 
from the fact that g represents the half width at half maximum in contrast to x. This 
quantity is mediated by the coupling with xo = g(1 — S22) so that Eq. results in 


2 
Syed = 


(3.9) 


Ag? , i 
Kardon + KL + 2iAw 


This very hand-wavy argumentation should not be understood as a full and exact 
derivation but rather as an explanation for the result, which is exactly derived by 


input-output theory in Ref. 


1 We note again that x in this work is always understood as full width at half maximum linewidth 
(FWHM). In the cited reference, the half width at half maximum linewidth is used instead. 
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Figure 3.5: Spectrum of [S]? for different internal linewidths. Calculation of avoided crossings 
according to Eq. (3.9). The resonance frequency of the first resonator, which couples to the transmission 
line, is fixed and the detuning of the second resonator is given on the x axis. The coupling of 
g/2n = 50MHz is the same for all images. The linewidths of the first resonator are chosen to be 
Ke/27 = 10 MHz = 2x; 1 /27t, giving an overcoupled resonator and therefore a good visibility in the 
amplitude data. The internal linewidth of the second resonator is chosen as xj» = 2/5 (left), ia = g 
(center) and £j» = 5g (right) to visualize the transition from strong to weak coupling. 


If we calculate the |S11| spectrum for different detunings Awı, Ac» and vary the 
participating linewidths, we see a transition of the frequency response depending 
On 2 2 Kmax, Where Kmax is the largest of the included linewidths. For the variation 
of ki», this is depicted in Fig. In the case xi» < g (left), the avoided crossing 
is well visible, but it blurs out when increasing the losses of the second system. 
For x1» > g (right), only a weakening of the first resonator’s resonance line can be 
observed when the systems are in resonance. In this limit, the excitation in the second 
system decays faster than it can swap back into the first system. In other words, 
the systems cannot couple resonantly but the second resonator is just another loss 
channel for the first one. This also explains the line broadening of the first resonator 
when tuned to the same frequency. As the anticrossing is only visible for g > Kmax, 
we call this regime the strong coupling regime. 


It is important to note here that the solution of Eq. is not affected by these changes 
of Kmax and always represents the energies of the avoided crossing, meaning that 
Eq. can only be taken into account for the strong coupling regime. 


The previous considerations only studied steady-state solutions when exciting one 
system stronger than the other one, such that they continuously exchange energy. 
Then, a beating, .i.e., a repulsion of the resonance frequencies as shown in Eq. 87) is 
observed. If the two systems are equally excited by microwave pulses, they oscillate 
in phase (or anti-phase) and no energy is exchanged. We studied this topic in-depth 
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by pulsed measurements on a microwave cavity coupled to a magnetic oscillator 
and published the results in Ref. [138]. 
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In this chapter, I give an overview of the experimental methods used throughout 
this thesis. I start with the cryogenic environment, which includes microwave 
and DC wiring and then go to the microwave measurement setups, which 
enable continuous wave and pulsed measurements. A remark on the home- 
grown software package QKIT concludes the measurement part, after which 
the qubit samples are presented. The chapter ends with a short summary of 
relevant fabrication methods. 


4.1 Cryogenic setup 


To set the ambient conditions necessary for the operation of superconducting 
quantum bits, we need to cool the quantum chip down to several tens of mil- 
likelvin and provide microwave and DC connections for control and readout. 
Such low temperatures are not only necessary to reach the superconducting 
state, but also to reduce the thermal population of the qubit and even more 
important the quasiparticle density in the superconductor, which would otherwise 
strongly limit the qubit coherence. 


4.1.1 Dry dilution refrigerators 


To reach the regime of about 20 mK in continuous operation, the choice of cooling 
systems falls almost exclusively on dilution refrigerators. These cryostats operate 
with a mixture of ?He and “He, where the two isotopes differ in their behavior at 
low temperatures. Below 2.2K, the bosonic “He atoms can fall into a condensate 
state, which is known as superfluidity. 3He atoms however are of fermionic nature 
and do nof] reach a condensate state, but stay in the liquid state with normal 
viscosity. Below this transition temperature, the two isotopes undergo a phase 


1 Actually, it was found that ?He reaches a superfluid phase at 2.5 mK by the formation of Cooper 
pairs . This is however not of practical use and not relevant here. 
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separation, one with a high and one with a low concentration of 3He. Due to 
their different mass densities, they also separate spatially. 


The small amount of *He atoms in the dilute phase can move in the superfluid 
liquid resistance free and therefore almost as in vacuum. The transition of 3He 
atoms from the concentrated to the dilute phase therefore equals an evaporation 
process, where the transition energy is taken from the thermal energy of the en- 
vironment. Since there is an equilibrium concentration of about 6% of ?He atoms 
in the dilute phase, this equilibrium would be quickly reached and the cooling 
power would vanish. Therefore, the ?He concentration has to be artificially de- 
creased to push the system away from equilibrium and get a continuous cooling. 


This is done by connecting the dilute phase to a still chamber, which can be heated 
and where a turbomolecular pump is connected. Due to their different partial 
pressure, mainly ?He evaporates in the still, which is then cleaned in a liquid ni- 
trogen coldtrap, compressed, and liquefied on its way to the dilution unit, where it 
is inserted into the ?He-rich phase. A schematic of this is flow is given in Fig. [4.1] 


Although this method allows in principle for arbitrarily low temperatures, 
heat input from the cryogenic circle and via mechanical connections limit the 
base temperature of bare ?He/^He dilution refrigerators to about 10 mK. With 
additional wiring for microwave and DC connections, this typically increases to 
a value of 10 — 20mK for multi-purpose cryostats. 


The pre-cooling to reach the phase separation temperature can either be done by 
reservoirs of cryogenic liquids like nitrogen and ^He, or by a closed cycle pulse 
tube refrigerator. The latter has the advantage that no infrastructural facilities 
for the liquefaction and storage of liquid nitrogen and helium are necessary and 
that the cryostat can be operated for weeks without user interaction, which is why 
it was used for this thesis. The pulse tube cooler acts analog to an inverted Stirling 
motor, where the working gas itself is used as a second piston, avoiding moving 
parts in the cold area. Requiring only electric power and cooling water, a two-stage 
pulse tube cooler can cool down to typically 3 K. To bridge the gap to the phase 
separation at 2.2 K, the ?He/^He mixture is compressed to 3 bar, precooled at the 
pulse tube stages and then expanded via a Joule-Thompson (JT) valve, leading to 
a sudden drop in temperature. In normal operation, the pressure is reduced to 
about 500 mbar and additional cooling power of the JT valve is not required. 
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Figure 4.1: Cryogenic measurement setup. A full photograph of the Oxford cryostat is shown on the 
left, where the according temperature stages are marked. On the right, the cryogenic microwave setup 
is illustrated, where the colors of the lines correspond to the choice of material for the coaxial cables: 
stainless steel (red), CuNi (blue), superconductor (green), and regular copper coaxial cables (black). 


In this thesis, experiments were carried out in two different cryostats of different 
manufacturers, one BlueFors LD 400 and one Oxford Triton 200. They differ in 
some engineering details and in their inner configuration but rely on the same 
physical principles for cooling. Great care was taken in all cryogenic setups to 
reduce the heat input on the temperature stages, since the minimum working 
temperature is limited by the cooling power. Therefore, any wiring between 
two stages was carefully selected in terms of material and wire cross section, 
to reduce the thermal conductivity. Also, Ohmic dissipation of current-carrying 
components has to be taken into account. 
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4.1.2 Cryogenic microwave setup 


Having a large metallic cross-section, coaxial wires can contribute significantly to 
the heat input at low temperature stages. Therefore, each of the three components — 
outer conductor, dielectric, and inner conductor — needs special consideration. Being 
metallic, the inner and outer conductor are the most obvious crucial parts in cryogenic 
microwave cables. While regular cables are based on copper, we use more resistive 
alloys like copper-nickel or stainless steel in our cryostats to reduce the thermal 
conductance. This increases at the same time the electrical resistance and gives rise to 
microwave losses. The increased electrical resistance is beneficial as its basis is formed 
by an increased electron-phonon-interaction, which helps to thermalize the electrons. 


Besides the electronic and the phononic temperature, we also take care of the noise 
temperature and power of the signal we send to the sample in the cryostat. The noise 
temperature of the input signal is approximately room temperature and commercially 
available electronics create signal strengths which are orders of magnitude stronger 
than what is needed for few photon operation in the cryostat. We therefore have to 
attenuate the signal, where it is obviously beneficial to put the attenuators at the cold 
stages of the cryostat, since every resistor creates Johnson noise corresponding to its 
temperature. If the signal is attenuated at N stages by an attenuation factor a; > 1 
at temperature Tj, the final noise temperature at the sample is given by 


E N Ti 
mcus c 
i-1 1lj-14j 


where Ty+ı = 300K is the signal input at room temperature and Tı the base 
temperature of the cryostat. Although it seems beneficial to pack all attenuation 
to the coldest stage in order to get a low noise temperature, this is experimentally 
not feasible since the attenuators would become very warm due to the dissipated 
signal and hence increase Tj. Typically, the total attenuation of around 60dB is 
distributed amongst the three coldest stages and additional 10dB are associated 
with the loss by the resistive coaxial wires. 


For both fridges used in this thesis, the attenuation was distributed onto the base, 
still and second pulse tube stage at temperatures of about 30 mK, 800 mK and 
3.5 K, respectively. For the BlueFors fridge, a combination of 20/10/30dB was 
used and for the Oxford fridge 20/20/20dB. This gives a noise temperature of 
TaB = 41mK and Tyo = 38mK. 


While this high attenuation is intended for the input signal path, this scheme is not 
applicable to the output signal path, where the signal power is in the low photon 
regime. Even the attenuation of the resistive coaxial wires would lead to the loss of 9 
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out of 10 photons. Instead, we use a coaxial cable made of superconducting niobium 
which provides ideal electric conductance while at the same time having no thermal 
conductance. This wire is used from the base plate to the 4K stage, where a HEMT?| 
amplifies the signal by about 40 dB and stainless steel cables can be used afterwards. 


The third component of coaxial wires is the dielectric, which is typically PTFE, also 
known as Teflon. Since PTFE is transparent in the infrared regime and infrared 
raditation from room temperature can create quasiparticles and therefore limit the 
coherence, some effort was spent in the field and in our group to develop filters 
for infrared radiation. While a significant improvement of coherence times was 
absent, some filters had mechanical problems leading to visible cracks and improper 
microwave transmission after a temperature cycle of the cryostat. The IR filters have 
therefore only been used in the BlueFors cryostat. 


At the base stage of the cryostat, we use a circulator to be able to perform a reflection 
51; measurement of our sample with our setup of dedicated input and output lines 
and a second one to shield the sample from the electronic noise coming from the 
HEMT. Additional band-pass filter are used to shield the sample from any unwanted 
microwave signals which could eventually pass the circulators out of their working 
band. Depending on the number of samples in the cooldown, also a 6-port microwave 
switch was used between circulator and sample, to be able to measure multiple 
samples without warming up and opening the fridge. 


4.1.3 Magnetic field bias 


Both fridges feature a superconducting solenoid, which can create a homogeneous 
magnetic field. In the BlueFors fridge, a commercially fabricated solenoid with 
5796 windings creates 47.5mT/A and in the Oxford fridge, a self-made solenoid 
of 10330 windings creates 81.5mT/A. 


In the BlueFors cryostat, the coil is mounted to the still stage such that the sample 
can be fixated to the base stage (see Fig. (4.2). This causes less heating on the base 
stage when ramping the magnetic field due to eddy currents and is beneficial in 
the case of a quench of the coil. We took great care in this case that no part of 
sample or wiring touch the solenoid body to avoid a thermal short. Also, the sample 
temperature is assumed to be slightly higher due to the influence of radiation from 
the solenoid. To this end, the coil in the Oxford fridge is mounted to the base stage 
and the ramping rate of the magnetic field was reduced here. 


? HEMT: High electron mobility transistor 
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Figure 4.2: Solenoid configuration in both fridges. a) The solenoid in the BlueFors fridge is mounted 
to the still stage with the sample reaching in from the base stage below. b) In the Oxford fridge, both 
solenoid and sample are fixated at the base stage. c) The sample is aligned by eye such that the chip is at 
the very center of the diameter. 


The samples in the solenoid are aligned by eye, which leaves a small probability 
for an off-center position and a tilt against the field lines. Although the pure qubit 
coherence as presented in Ch. [6] might be improved by a better alignment and 
the use of compensation coils, we believe that this is not a realistic scenario as 
soon as magnetic materials are introduced on top of the qubit, which will alter 
the magnetic field configuration. 


To supply the comparatively high currents down to the base stage, different materials 
are used: Between room temperature and the 70K stage, massive copper wires are 
taken which increase the thermal load on this stage. For the lower cooling power of 
the 3.5K stage, this is not feasible and ribbons of high temperature superconductor 
on a stainless steel matrix are used here. Below, a NbTi superconducting wire in a 
Cu or CuNi matrix is used. The technical realization of these connections is in detail 
described in the PhD thesis of my colleague Marco Pfirrmann [140]. 
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Figure 4.3: Filter box for DC current bias. Schematics (left) and photograph (right) of the DC bias filter 
box at the 3.5K stage. The filter is divided into 3 x 8 wires, where the ones on the left only have a RC 
filter, the ones in the center have an RC filter and a 1: 1 current divider and the ones on the right simply 
pass through. Note that the image only shows the upper side of the layout and the bottom side looks the 
same. For the RC filters, multiple capacitors are used in parallel, as the ones with a high capacitance have 
a bad high frequency transmission. 


4.1.4 DC bias 


To be able to tune the transition frequency of the tunable sample described in Sec. [4.3] 
a current has to be provided down to millikelvin temperatures. The focus here lies not 
on the high current regime but on additional filtering, which removes high frequency 
noise and provides a constant current as stable as possible. The wiring between 
room temperature and the 3.5 K stage is made from BeCu, the rest down to the base 
stage by NbTi. All wires are multi filament braided looms, where adjacent lines are 
twisted to reduce coupling of electromagnetic noise. To filter out noise, a filter box 
was designed for the 3.5 K stage, shown in Fig. {4.3} This box features three different 
filter configurations, to allow for a flexible usage. While one part is only equipped 
with RC filters, the second features an additional current divider to divide the current 
and with that the current noise on low temperatures. 
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4.2 Microwave measurements 


4.2.1 Spectroscopic measurements 


For many characterization experiments of superconducting resonators and qubits, it is 
sufficient to probe the system in the steady state. Using commercially available devices 
which are designed for this task, a standard spectroscopic setup can be assembled 
from few parts. Coming along with this structural simplicity, also the control and 
data acquisition is very simple with this setup. 


A typically used setup for spectroscopic measurements is shown in Fig. It 
consists of a vector network analyzer (VNA) as readout device which probes the 
resonance spectrum of the resonator and additional microwave sources to apply 
probe and drive tones. To reach a high drive amplitude despite the high attenuation 
on the input lines, we use a high power microwave amplifier at room temperature 
for the measurements presented in Ch. |5| The choice of attenuation at the output 
of VNA and microwave source shows a typical configuration and is adapted to 
the needs of the experiment, respectively. 


4.2.2 Time resolved measurements 


We use a time domain setup to measure the characteristic timescales of the qubit, 
which is described in the next section. Beforehand, we will shortly get an overview, 
which physical quantities can be measured by which pulse sequence. 


In general, all time domain measurements consist of a sequence of pulses, after which 
a readout tone is applied. The readout projects the qubit state onto the z axis in the 
Bloch sphere and therefore on one of the states |0) or |1). This pulse sequence is then 
repeated for a number of different pulse lengths or delay times, depending on the 
protocol. To acquire statistics of the measurement, this protocol is then repeated several 
thousand times and the result is recorded. It is important to repeatedly loop through all 
pulse configurations instead of directly repeating each configuration multiple times, to 
ensure that fluctuations and slow drifts affect all acquired data points in the same way. 
After the readout, a waiting time before the next pulse is implemented to let the qubit 
and resonator relax back to their ground states. The three-dimensional representation 
of the complex qubit state was performed by quantum state tomography [77]. 
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Figure 4.4: Spectroscopic setups. A two-tone (left) and three-tone (right) microwave setup for spectroscopy 
measurements. Both setups include a vector network analyzer (VNA) to probe the resonator and a probe 
microwave source to probe the qudit level transitions at different powers and frequencies, Pp, wp. In 
the three-tone setup, an additional drive microwave source with an high power amplifier is included. 
The attenuation indicated in the scheme represents typical values and is adjusted depending on the 
experiment. 


Rabi measurement 


The Rabi measurement consists of an excitation pulse at the qubit frequency, with 
increasing pulse length t, followed by the readout. We recall from Eq. that 
this corresponds to applying ö, for the time t and the qubit state starts to rotate 
around the x axis. This can be seen in Fig. [4.5] a) as a cosine-shaped oscillation in 
the readout value. Rabi measurements are frequently performed to measure the time 
tz, for which the pulse of distinct power has to be applied to rotate the qubit state 
from |0) to |1). The readout signal is then expected to be 


(02) =e7*/TRabi cos(27tT/2tz), (4.1) 


with Trabi = (TI + T; 1). 


Relaxation time measurements 


To measure the relaxation time Ti, a zt pulse is applied to bring the qubit in the |1) 
state. After waiting for a time At, the remaining qubit excitation is read out. Since the 
qubit does not acquire any phase information in the relaxation process, the trajectory 
on the Bloch sphere goes directly through the origin of the sphere, see Fig. [4.5|b). 
If the initial 7t pulse did not perfectly match in length or frequency, the initial state 
is different, but we still observe an exponential decay of 


een, (4.2) 
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Figure 4.5: Typical time domain measurements. Each subfigure shows the pulse sequence at first, a 
tomographic image of the qubit state evolution in the center and the typical measurement outcome on the 
right. The time in the data points is color coded, from red to turquoise. a) A Rabi experiment which rotates 
the qubit state around the x axis. b) For the relaxation time measurement, the qubit is prepared in the |1) 
state and the time evolution is recorded. It is well visible that no additional phase information is acquired 
by the decay. c) A Ramsey experiment with a detuning between drive and qubit frequency. The Bloch 
sphere in the center shows the evolution before the second 77/2 pulse is applied, which enables to measure 
the value of (0). d) Demonstration of an all-axis control of a fast tunable sample. All measurement data 


and plots are taken from [77]. 
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4.2 Microwave measurements 


Ramsey measurement 


For a Ramsey measurement, the qubit is prepared in the state (|0) +|1)) / 2 by 
applying a t/2 pulse. After a time At, a second t; /2 pulse is applied and the state is 
read out. If the qubit state is unperturbed in the time At, we will always measure the 
|1) state. From Eq. we know however that even if no drive is applied, the operator 
A0; effectively acts on the qubit state, where Aw is the detuning between drive 
and qubit frequency. Also, any relaxation or dephasing will affect our measurement 
outcome, resulting in |0) before the last t; /2 pulse. This measurement scheme is 
depicted in Fig. c) and follows the following formula: 


(0;) = —e-^'/T? cos A. At, (4.3) 


where we explicitly keep the index R on the characteristic timescale I to avoid 
any ambiguity and confusion of the star notation. 


To measure the Ramsey dephasing time, it is typically beneficial to not set Aw = 0, 
but chose a small detuning such that several oscillations are visible. This gives 
a much clearer envelope and the possibility to distinguish between frequency 
mismatch and dephasing. The detuning A, is technically not the actual detuning 
between drive and qubit frequency, but the rate Ag/ At at which the phase between 
the two t„/2 pulses is adjusted. 


Due to the oscillating nature of this measurement result, Ramsey measurements are 
commonly used when a sensitive measurement of any change in qubit frequency is 
needed. If the TR time is long enough, the time within which the qubit can accumulate 
phase and therefore the resolution of the Aw measurement grows. 


Hahn echo measurements 


If we assume that the dephasing is caused by slow changes in the qubit transition 
frequency, this means that Aw is different for each measurement and the averaged 
vector on the Bloch sphere has zero length after short At when averaging over 
some seconds. This can be compared to averaging the outcome of multiple Ramsey 
experiments with long IS times but different frequency, which will result in a very 
fast decay of the resulting envelope. 


If we assume that the fluctuations are so slow that Aw does not change between the 
first pulse and readout, we can counteract these fluctuations by applying a t pulse in 
the middle of the waiting time At. This pulse flips the qubit state from (|0) + |1)) / v2 
to ([0) — |1)) / V2, which effectively changes the sign of the acquired phase. The phase 
acquired after the t pulse then compensates this and independent from a Aw, the 
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Figure 4.6: The time domain setup. To enable pulsed measurements, a four channel digital-to-analog 
converter (DAC) is used, where each pair of channels creates the pulse envelopes for the IQ mixers. By 
using the same microwave source for up and down conversion of the readout signal, we ensure a common 
frequency and phase reference. The signal is recorded by a two-channel analog-to-digital converter (ADC). 


state will end up in (|0) — |1)) / V2. Faster fluctuations in Aw cannot be compensated 
by this type of experiment. Due to the independence from the absolute value of Aw, 
the measurement signal shows no oscillation but a decay following 


(0;) = e AUTE, (4.4) 


The name “echo” measurement comes from NMR experiments, where a spin ensemble 
is used. In this case, all spins oscillate in phase again after At/2 after the tz pulse, 
such that they emit a measurable microwave signal. This signal can be seen as 
echo of the initial excitation pulse. 


4.2.3 Time domain setup 


To perform the measurements introduced in the last section and therefore measure 
the characteristic timescales of the qubit, we use a time domain setup, which is 
depicted in Fig. The pulse envelopes are generated by a four-channel arbitrary 
waveform generator (AWG) and fed to the I and Q inputs of two IQ mixers, one for the 
qubit manipulation pulses and one for the readout pulse on the resonator frequency. 
The two signals are then combined by a directional coupler. The RF inputs for the 
mixers are generated by two local oscillators (LO). Since IQ mixers are nonlinear 
and very non-ideal components, the technique of single sideband mixing is utilized, 
where the pulse envelopes are multiplied by a sine and cosine signal and fed to 
the I and Q inputs of the mixer, respectively. The oscillating signals typically have 
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a frequency of wyg/27 = 63 MHz for the readout channel and wyg/27 = 75 MHz 
for the manipulation channel. This shifts the generated tones in frequency, which 
means that we can shift the frequency of the LO by the IQ frequencies and any 
leakage of the LO due to non-ideal mixers will not be resonant with any transition 
on the qubit chip. Constraints for these IQ frequencies are the output bandwidth 
of the AWG, the bandwidth of the IQ mixers and the anharmonicity of the qubit, 
since we do not want to excite any other transition with a leaking LO. Also, the IQ 
frequencies should be different from each other to account for possible cross-talk 
of the AWG channels. The generation and optimization of phase-controlled pulses 
with this setup was in-depth described in Ref. [77]. 


The reflected readout signal from the resonator is amplified and down-converted by 
an IQ mixer with the same reference LO as for the up-conversion. After a low-pass 
filter to eliminate high-frequency mixing products, the signal is amplified and fed 
into a two-channel digitizer. The acquired data for multiple iterations of the pulse 
sequence are averaged on a computer and the amplitude and phase value of the 
complex fast Fourier transformation at the IO frequency of the readout channel 
is calculated. This result corresponds to the S11 signal of the resonator at a fixed 
frequency and can be used to detect a shift of the resonator resonance frequency 
and therefore the qubit level population. 


4.2.4 Data evaluation and measurement software 


For the experiments throughout this thesis, a lot of effort has been spent to enable 
and facilitate the fast and easy acquisition of measurement data, alongside with a 
configurable data viewer, fast evaluation and fitting tools and a flexible database 
to organize all measurement data. This software package, named QKIT, is written 
in python and freely available under Ref. [141]. 


This toolset enables to repeatedly acquire spectroscopic or time domain data, where 
almost any imaginable automation can be included. Allowing for an automated fit 
after the measurement, the calibration data of a Rabi experiment can for instance 
be automatically fed to the next T; measurement, which requires the t; time for an 
excitation pulse. Especially for the magnetic measurements, it was from time to time 
necessary to thermally cycle the sample, i.e., bring it over the critical temperature 
for a short time to eliminate all trapped flux. This procedure can be easily included 
into the measurement scheme and the measurement automatically continues, when a 
satisfactory base temperature is reached again. Finally, the measurement sequence 
automatically logs and saves all settings of the measurement devices, leading to a 
semi-automated lab book and reproducible measurement schemes. 
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Figure 4.7: The measurement database and quickplot feature of qKIT. Above, an excerpt of the 
measurement database is shown which lists after the unique ID of each measurement the starting 
timestamp, comments, a user defined rating, relevant device parameters and fit results of the data. AII 
user-defined data can be directly changed in the table and sorting and filtering can be enabled to restrict 
the range of listed measurements. On the bottom, two windows of different types of measurements are 
shown. To the left, the anti-crossings of the tunable transmon with the resonator frequency are shown 
in amplitude and phase as color coded data. Any constant background is eliminated in this view by 
subtracting the average over the x direction. While this coarse method is not directly applicable for 
pictures intended for publications, it increases the contrast of the measured data and allows to clearly see 
even small changes. To the right, a reflection measurement of a resonator is shown. 


When acquiring lots of data by varying different tuning parameters, it is extremely 
helpful to be able to quickly screen and review the recorded data to see relevant 
features and be able to scan interesting regions in more detail. A recently developed 
quickplot function enables to plot the relevant data for each measurement in the 
database within fractions of seconds (depending on the size of the acquired data) and 
automatically apply filters for background removal, phase wrapping or elimination 
of runaway values. Each measurement can be identified with a six-digit unique 
ID and the meta data of all measurements is stored in a measurement database. 
This database allows to filter for samples, measurement type, fitted values, or any 
other parameter and enables to give ratings for the quality and relevance of a 
measurement. In sum, these tools help to keep an overview over all acquired data 
and find relevant measurements again quickly. A demonstration of this database 
and the quickplot interface can be seen in| Fig. 4.7. 
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Figure 4.8: The two qubit samples used throughout this thesis. Micrographs of the two types of 
transmon qubits. On both images, from left to right, the microstrip microwave transmission line, the 
readout resonator and the qubit is shown. Left: The non-tunable transmon made from TiN, fabricated at 
NIST. Right: The gradiometric concentric transmon is in principal tunable by a current flowing through 
the flux bias loop on the right. This loop induces unequal fluxes in the qubit’s tuning loops and therefore 
tunes the qubit, as described in Sec. [2.5.1] The flux bias loop is executed as coplanar waveguide structure 
with a ground plane directly next to the conductor to reach a 50 Q microwave impedance. To ensure a 
proper grounding, two bonds are connecting both sides of the ground plane. 


4.3 Qubit samples 


For the experiments in this thesis, two different versions of the concentric transmon 
have been used, which also differ in choice of material. A micrograph of the 
two samples is given in Fig. 


The first chip was fabricated at the NISTP|by the group of Martin Sandberg and David 
Pappas, and gratefully delivered to us for measurements. Besides the suppression 
of the dipole moment as described in Sec. the focus for this sample was the 
use of titanium nitride, a material engineered for a high critical temperature 
and low microwave losses [143]. TiN has the additional advantage to be stable in 
humid air at room temperature and is therefore not likely to form surface oxides. 
The Josephson junction is however realized by a Al/ AIO,/ AI structure, which has 
been found to give the best coherence. This two-step process brings additional 
features which are explained in Ch. 


From spectroscopic measurements, we find the fundamental qubit transition frequency 
as Wo, /2zt = 4.755 GHz and determine from the higher level transitions a charging 
energy of Ec /h = 197.7 MHz and a Josephson energy of Ez/h = 15.5 GHz. This results 
in a ratio of Ej/ Ec = 78.4. The full spectrum of this qubit can be found in Fig. [4.9] 


3 National Institute for Standards and Technology in Boulder, Colorado 
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Figure 4.9: Transition frequencies of the TiN transmon sample. Results from a two-tone experiment to 
see the qudit level transitions under increasing qudit probe power. Since multi-photon transitions start 
to occur at higher powers, the transitions can be identified and Ej and Ec can be extracted. Due to an 
increased sample temperature, also transitions starting from the |1) state can be seen (dash-dotted lines), 
which fit nicely to the values expected from the qubit parameters. Adapted from in]. 


Capacitively coupled to the qubit isa A/2 resonator at wr = 8.573 GHz by a coupling 
strength of g/27 = 71.5 MHz. On its other end, the resonator couples to a 50 
microwave transmission line in microstip design. 


The second type of sample was fabricated at KIT by Jochen Braumiiller and is in-depth 
described in Ref. [99]. As explained in Sec. this design features two JJs in a 
gradiometric SQUID geometry, which makes the qubit transition frequency sensitive 
to field gradients but not to homogeneous fields. For the experiments in Ch. 
at higher magnetic field values however, we found that the inhomogeneity of the 
magnetic field is large enough to significantly tune the qubit frequency. Together 
with a fluctuating current of the solenoid’s power supply, the qubit transition was 
strongly broadened and slowly drifted with time. This made the qubit unusable for 
magnetic field sweeps with the given magnetic setup. The experiment was therefore 
executed with a non-tunable transmon sample, where one of the JJs was open due 
to fabrication issues. This qubit sample had a fixed frequency, which only weakly 
changes in the magnetic field. Since the exact same qubit design was used as for 
the tunable qubit, all experiments in Ch. [7| can be easily repeated with a tunable 
transmon once the bias magnetic field is stabilized. 


The transition frequencies of this sample are shown in Fig. from which we 
can calculate a charging energy of Ec/h = 182MHz and a Josephson energy of 
Ej/h = 32.7 GHz. For this sample, we get Ej/Ec = 179.7, being even deeper in 
the transmon regime than the TiN sample. In addition to a A/2 readout resonator 
at węg = 8.78GHz which is also coupled to a microstrip transmission line, this 
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Figure 4.10: Transition frequencies of the all-Al qubit sample. Results from a two-tone experiment 
to see the qudit level transitions under increasing qudit probe power. We again see the multi-photon 
transitions starting at higher probe powers. Due to the high EjEc values extracted from this data, this 
sample is even deeper in the transmon regime than the TiN sample, which can be seen in an almost 
constant frequency difference between adjacent multi-photon transitions. 


design includes a 50€) impedance matched flux bias line to be in principle able 
to tune the qubit by fast pulses. 


4.4 Sample fabrication 


4.4.1 Qubit samples 


Since for this work, a variety of samples generously provided by other people 
was available, there was no need to fabricate own qubit chips. To understand the 
interference effects presented in Ch.{6|better, the fabrication process is sketched shortly 
here. The fabrication procedure of the JJs with shadow angle evaporation is shown 
in Fig. For the TiN chip, a two-stage process is used, where first the resonator, 
transmission line and qubit islands are structured in TiN by optical lithography. In 
the Al chip, all structures including the junctions are integrated in one layout. The 
resist system for electron beam lithography consist of two different resists, which 
have different sensitivity to electrons. It is therefore possible to create an undercut 
under a free standing bridge. In a first step, aluminum is evaporated under an angle, 
where no material is deposited in the shadow below the bridge. In the second step, 
a controllable partial pressure of oxygen is applied to the chamber, which oxidizes 
the surface of the aluminum. Next, Al is evaporated under a different angle, which 
leads to a small overlap contact of the two Al layers. Finally, the resists are removed, 
which takes away all unwanted metalization. 
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Figure 4.11: Fabrication process of Josephson junctions. In a preceding step, the large metal parts are 
created from TiN and the two resists (light and dark red) are structured. The resists form a free standing 
bridge which is being hold in the third dimension not shown here. First, aluminum is evaporated under an 
angle (blue arrows), whose surface is oxidized in the second step by an oxygen atmosphere. Afterwards, a 
second layer of Al is deposited (blue arrows), where a small overlap contact between the two layers is 
formed (yellow area). After removing the resists, the device is ready to be used and can be contacted. 
The other Josephson contacts (black areas in step 4) are commonly neglected due to their large size but 
become relevant in Ch. [6] 


4.4.2 Magnetic samples 


The magnetic material used in Ch.|7|was grown by sputter deposition with synchronous 
co-sputtering of Iron and Cobalt. To get a stoichiometry of about Feo.7Coo., the sputter 
guns are adjusted to have a power of Pre = 100 W and Pc, = 37 W, given by measured 
individual deposition rates and known material densities. After this process, a 
protective coating of AIO, was sputtered to prevent a fast oxidation of the cobalt-iron 
alloy in the further processing steps and measurements. The fabrication parameters 
for the optical lift-off process are given in Appendix 


After structuring the ferromagnet, a second optical process was applied to structure 
rectangles of resist at the edges of the individual chips. These structures are baked on 
a hot plate at 230°C for 6 minutes to glass the resist, making it insoluble in solvents. 
Additional parameters to this process are also given in Appendix [B] 
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superconducting qubits 


Sensing a microwave signal directly on the quantum chip at millikelvin temperatures 
is a challenging task, especially when the photon energy is widely undefined. In this 
chapter, I present a method to perform these measurements based on the AC Stark 
effect on a multi-level superconducting quantum circuit. After a general introduction 
to quantum sensing, I present the AC Stark effect first in general and then in detail for 
a multi-level system by means of analytic and numeric theory. After a comparison with 
measured data, I detail on the sensing scheme and present a first proof-of-principle 
experiment, leading to more advanced measurements. The results of the experiments 
in this chapter lead to the publication of Ref. [71]. 


5.1 Quantum sensing 


Although quantum technologies are used for sensing and measurement devices 
already for several decades, the broad field of quantum sensing is still highly 
subject to active research and development. An outstanding example in the field of 
superconductivity is the concept of a SQUID (see Sec. [2.4.1), which was developed 
over 55 years ago and is still very frequently used in various fields of research, 
manifested by its appearance in more than one publication per week! on average. 
Similar to semiconducting sensors which had applications long before semiconducting 
computers became relevant, it becomes more and more apparent that quantum sensors 
are of general use already before quantum computers are considered to be state of 
the art. The amenities of the development of new quantum sensors compared to 
quantum computers is that all spurious and unwanted effects in quantum computing, 
leading to noise, decoherence, and gate errors, can in principle be used to build 
a quantum sensor which exactly measures these decoherence effects and spurious 
influences. The weakness of quantum computing can therefore be seen as the strength 


1 Source: Preprints published on arXiv.org, based on the data for Jan-Nov 2019. 
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of quantum sensing. This means, on the other hand, that the emerging topic of 
quantum sensing can be and frequently is used to reuse spin-off results from research 
on quantum computers, where a restriction in the applicability of quantum computing 
was found and turned into a feature for quantum sensing. Being a valid approach 
for the development of quantum sensors, this motion in the field blurs out the name 
^quantum sensing". Although being not clearly defined, usually this term is used 
when one of the following three features is used to measure a physical quantity [29]: 


1. A quantum object, 
2. quantum coherence or 
3. quantum entanglement. 


While strictly speaking only the third point really exploits the quantum nature of the 
sensor and opens the possibility to reach a sensing resolution beyond the classical 
limit, the other types of sensors also have their right to exist and are commonly entitled 
as quantum sensors. Despite not being "super-sensors", they are still able to beat 
the resolution and flexibility of classical and established sensors. Especially when 
examining other quantum systems, they can eventually be the only sensors applicable 
in terms of size, energy scale and coupling mechanisms. 


The list of quantum sensors which are nowadays routinely used is very long, and 
some of them are not even recognized as quantum sensors any more, since they are 
so commonly used. One very basic example is the usage of reference substances 
for ESR and NMR?|experiments, which have a characteristic frequency-to-magnetic- 
field relation and sharp resonance features. This makes them useful as magnetic 
field reference sensors, frequently used to calibrate the magnetic field created by 
permanent magnets or a magnet coil. 


Speaking of NMR experiments, the whole human body can be turned into a quantum 
sensor in the context of magnetic resonance imaging (MRI). Here, the NMR properties 
of the hydrogen atoms in the body are measured and spatially resolved, giving 
insight to the chemical environment of the H atom. Measuring T; and T? makes 
it possible to draw conclusions on the type of tissue, water and fat contents, and 
more. By using special techniques it is also possible to characterize motion and 
flow as well as metabolic processes in the body in real time, making MRI to a 
valuable examination tool in medicine [145]. 


One of the most broadly used types of quantum sensors with respect to sensing 
quantities and application areas are nitrogen vacancy (NV) centers, which represent 


? ESR: electron spin resonance, NMR: nuclear magnetic resonance 
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electronic spin defects in the crystal structure of diamond. They have been demon- 
strated to be able to measure magnetic fields, electric fields, pressure, and temperature 
and can be used as bulk material or in the form of stable single nano crystals. As 
nano crystals, they can be added to a scanning probe, which allows for sub-100 nm 
resolution in magnetic field resolution of surfaces. Or they can even be inserted into 
living cells where they allow for an in-vivo temperature measurement [29]. 


Similar to that, SQUIDs can be mounted onto a scanning tip where they are used 
for measurements of magnetic field, current and temperature with ultra high spatial 
resolution. By fabricating the nanoSQUID directly onto the tip, SQUID loop diameters 
of less than 50nm have been demonstrated [146] 147]. 


Exploiting not only quantum effects but also the sensitivity and vulnerability of a 
coherent quantum state, superconducting quantum bits can also be used as quantum 
sensors and tools for investigation. Compared to the previously presented nuclear or 
electron spins, superconducting qubits are macroscopic quantum objects that can have 
large magnetic and electric dipole moments and are therefore particularly suitable 
for coupling and sensing experiments. In a great variety of experiments, they have 
been used to study for example environmental noise input or spurious coupling to 


defects by means of T; relaxometry or resonant coupling 148]. 


To study an oscillating physical system of interest by means of a superconducting qubit, 
usually strong coupling is desired, which enables a coherent exchange of excitation 
and therefore information about the system. The necessary coupling strengths can 
however be challenging to achieve and some effects are only visible when the two 
systems are close in resonance, like the level repulsion described in Sec. 


Instead of a coherent excitation exchange, it would be beneficial to be able to use the 
qubit to sense the oscillation of the second system, which would reduce the severity of 
the second system's linewidth and therefore the need to be in the strong coupling limit. 


Being able to detect arbitrary microwave signals by a qubit would also open new 
ways to characterize and optimize environments for quantum experiments, especially 
with respect to quantum computation. This sensing scheme would allow to directly 
measure the cross-talk from all other elements on a quantum chip to the qubit, 
which could otherwise only be characterized in the resonant case or by means of 
simulations. Likewise, it is practically impossible to measure the full frequency 
dependent transmission from the microwave source to the qubit position by other 
means, as the wiring in the cryostat as well as the transmission on the quantum chip 
change when cooling down from room temperature to several millikelvin. 
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To construct such a microwave sensor from a superconducting qubit, we will in the fol- 
lowing study the influence of an oscillating field on a quantum bit. This is generally de- 
scribed by the AC Stark effect, which will be explained and detailed in the next sections. 


5.2 The AC Stark shift 


In analogy to the Zeeman effect which splits atomic transitions in a magnetic field, 
Johannes Stark studied the influence of a static electric field on these transitions 
and Schrédinger later derived a theoretical treatment for this effect [150]. It 
is since understood as the splitting and shift of energy levels from the interaction 
of a molecule’s dipole moment with the electric field. To observe this effect, strong 
electric field strengths are necessary, being experimentally hard to achieve and 
control. The DC Stark effect was therefore not as extensively studied as its magnetic 
counterpart, the Zeeman effect. 


With the introduction of lasers however, it was possible to strongly increase the limits 
of experimentally achievable AC electric field strengths and hence study the AC Stark 
effect. Experimentally, the shift and splitting of the D absorption line in potassium 
under a strong laser drive could be observed [151], and many other experiments 
followed including atoms or molecules in an optical cavity. 


With the formation of circuit quantum electrodynamics (OED) in analogy to cavity 
OED, it was obvious that also an artificial atom, i.e., a superconducting qubit, should 
show a shifting of its energy levels when strongly populating the coupled readout 
resonator [76], which was then also experimentally observed [152]. 


In both atomic QED and circuit QED however, the term “AC Stark effect" is sometimes 
generalized to all experiments including a strong drive tone or field [153], resulting in 
the description of Rabi or Autler-Townes splittings [27]. In this thesis, we will focus 
instead on non-resonant drives and study their influence on the qubit transitions. 


5.2.1 Resonator-induced AC Stark shift 


Following the analogy to cavity QED, the AC Stark shift is commonly studied as the 
influence of a coupled resonator. In many publications, it is therefore seen as the coun- 
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Figure 5.1: AC Stark shift caused by resonator population. We populate the readout resonator with 
different photon numbers by increasing the power Pr at the VNA to read out the resonator. The red 
line is a fit of the qubit transition frequency to the data after Eq. from which we expect a linear 
power dependence. This can be seen in the right plot, where the power is given in linear power units 
(mW). The left plot with the power given in logarithmic units (dBm) however renders the low photon 
number regime in greater detail and is therefore typically used. The photon number indicated in the 
upper axes is found by Awg = (ä'A) g?/Ar, where g/27 = 71.5 MHz and A;/27 = 3.818 GHz is known 
from other experiments. This means a single photon shifts the qubit frequency by Awg = 1.3MHz. 
We also see that the single photon regime can be reached by setting the output power of the VNA to 
P, ~ -21 dBm = 7.5 iW. Adapted from 


terpart to the Lamb shift, presented in Eq. (2.11) for the qubit readout [76]. By simply 
shifting the assignment of the x0: 'á term to the qubit part, the effect is directly obvious: 


^ "m A 
Her/hz (wr + x02) ata + z (won +x ) ôz (2.11%) 
1 à 
~ (wr —— )0 à (wort x + xâ â) às. (5.1) 


^ 


Now, a shift of the qubit frequency Awg = xá'áà proportional to the number of 
photons in the resonator can be seen, where we recall that the dispersive shift 
is given by x = g?/A,. 


This behavior is experimentally demonstrated in Fig. where the power of the 
readout signal was steadily increased, resulting in an increased number of photons 
in the resonator and eventually in a shift of the qubit transition frequency. We recall 


2 
from Eq. that the population in the resonator is given by (1) — Pin zug 7 This 
dependence is commonly used to adjust the power level for the readout tone to 
be in the single photon regime and to reassure that the qubit transition frequency 


is not influenced by the readout. 
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For real sensing applications however, the effect in this form is not of great use, as we 
can only detect the photon number of a coupled harmonic oscillator and the coupling 
constant g as well as the detuning of the systems have to be known. It also depends on 
the internal quality factor which is known to be power-dependent [154], giving an ad- 
ditional degree of freedom and no usable relation between power and photon number. 


Remembering the initial description of the AC Stark effect however, we realize that only 
an alternating electric field acting on the atom is needed. The fixed frequency cavity 
was solely introduced in atomic QED to increase the field strength and compensate 
for the low dipole moment of single atoms. 


5.3 Multi-mode AC Stark shift 


To access a more generalized approach to the AC Stark shift than the cavity-atom- 
model, we will expand our analysis in the rest of this chapter by the following factors: 
Instead of a resonant cavity system, we will treat a more universal drive field, which 
experimentally does not limit us to pre-defined elements or frequencies and lets us 
describe a greater variety of systems and processes. Also, the atomic part of the system 
will no longer be restricted to a true two level system, but we allow a more universal 
multi-level system with a small anharmonicity. This gives a better description of the 
system we are using experimentally and finally makes more observables accessible, 
which can later be used for an enhanced sensing scheme. 


In the following, we will access this regime using two different theoretical approaches: 
First, we start with a theoretical analysis of the system. Using a simplified model for 
our anharmonic system and applying the approximations of perturbation theory, we 
find an analytic formula for the qudit transition frequencies under the influence of the 
drive. This relation provides additional insights and allows us to derive simplified 
expressions for special cases. Due to the implied approximations however, these 
solutions cannot fully describe our experimental observations. 


For this reason we utilize numerical simulations in a second step, including a more 
detailed model of our qudit system and excluding the approximations of perturbation 
theory. Yielding good results on the exact problem we want to solve, this method 
however comes with the drawback of only calculating a solution for the given starting 
parameters, without any gain in physical understanding. Both paths therefore have 
their own benefits and will be presented in the following. 
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5.3.1 Approximate analytical solution 


To derive an analytical expression for the qudit transition frequencies under the 
influence of a drive, we start with the approximated transmon Hamiltonian of 
Eq. (2.29), where a constant anharmonicity is assumed: 


Hg/h = wab*b — 86 (bb +1). (5.2) 


For simplicity, we introduced the harmonic fundamental transition frequency Wg = 
y 8EcE;/h and the anharmonicity y = 2Ec/h and neglected the constant terms. We 
model our drive field with amplitude Ap and frequency wp by the Hamiltonian 


Hp/h = Ap (b + b*) COS Wpt. 


The amplitude of the drive Ap is here given in units of an angular frequency, which 
corresponds to the usage of the Rabi drive amplitude A in Eq. (2.9). On resonance, this 
is the frequency of the observed Rabi oscillations, meaning that the qubit changes from 
[0) to |1) within the time At = 71/ A and the energy of one fiw is then stored in the qubit. 
To satisfy an experimental physicist’s understanding, a simple power-to-amplitude 
relation can hence be given by P = Ahw/r. 


The beauty of choosing the amplitude for our analysis is that it describes the 
power directly applied to the qudit, independent of any coupling or frequency- 
dependent transmission function. As we will later use the readout resonator only 
to probe the qudit transition frequencies, we can neglect it here for our analysis. 
The full Hamiltonian then reads: 


H/h = wabtb- Toto (bb +1) + Ap(6 + 6") cos wp. (5.3) 


We now transform the Hamiltonian into a frame rotating with the drive frequency 
Wp by applying the unitary transformation 


u(t) — eiwpb'bt 


resulting in 


where fast rotating terms again have been neglected in a rotating wave approximation. 
Before we can treat this Hamiltonian by perturbation theory, we have to apply a 
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displacement operator Dy = gib —a*b with DDD, = b + x in order to make the drive 
term small compared to the other terms. This results in 


Ég/h = D'HpD,/h 
= (wa sige H = lal?) p*b T (66)? - (z*b + biz) 


-i (a*b*66 + ab*6*6) T (g2pt2 + «**5) + const, 
where the parameter ¢ was introduced as 
- Y AD 
PE wp 5 (1 [el] ees (5.4) 


We can now chose the parameter « in the displacement operator such that we get 
Ç = 0, which eventually cancels the terms linear in b and 6+ in Eq. (5.4). With 
this choice, the Hamiltonian reads 


Hg/h = (5.5) 
= a e 2\ pte "Y gto "Y, 2042 | 92524 Y (x ptp2 „+27 
c wp - 1 lal?) 6*5 4 (b*6)? — 1 (ab? + 0°76?) > (wb + ab b), 
—— Mmmm 
AS /h Ay/h 


where we already divided the Hamiltonian into the unperturbed part H9, whose 
eigenenergies E? can be easily calculated, and the perturbation Hz, which we will 
treat with perturbation theory. 


Evaluating the individual contributions, a closer analysis shows that higher orders 
of perturbation theory have to be included for strong drive amplitudes to get 
better results. Therefore, we evaluated the calculations up to fourth order, which 
can be found in Appendix 


(5.6) 


The full expansion of these products would become lengthy and gives no addi- 
tional insight. It is therefore not presented here. For the later presented analytical 
data however, the whole expression was included in the calculation. In order 
to understand the physical meaning of this findings, we restrict our analysis to 
corrections up to first order. We get 


Ei /n = (a — wo - T - yat?) k- Te, (5.7) 
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as the first order corrections vanish with E} = (k| Hj, |k) = 0. If we recall that we 
set € = 0 after Eq. (5.4), we find for the parameter a: 


ya|a|? + 20Ap — Ap = 0, (5.8) 


where we introduced the detuning between the drive and the first qubit transition 
Ap = wp — wg + Y/2. This gives in general three solutions for a. Since we limit our 
studies to the non-resonant case, it makes sense to restrict the drive frequencies to 
a regime where it cannot excite level transitions. When the higher level transitions 
are below the fundamental transition in frequency (y > 0), we want the drive to be 
above the first qubit transition (Ap > 0) and vice versa. It therefore makes sense to 
study the case Ap/Y > 0, which is true for the later presented measurements. We 
find that « has only one real solution in this regime: 


Se 2/3 
21/3 ( /81Aby +964} +9Ap V7) -4034p 
a= A 


(5.9) 
1/3 
(36 [37° (274b y + 3245)] ^ +3244p7?) 


For the other case, where Ap/y < 0, the number of real solutions for « depends 
on the drive strength Ap. Below a critical amplitude Ap, three real solutions 
exist, where two of them vanish above Ap;crit- This case is however unpractical for 
sensing applications, as already mentioned. 


With this theoretical analysis, we can make a sanity check by examining Eq. in 
the case of very large detuning, Ap > Ap, y, which eventually corresponds to the 
dispersive limit previously regarded in Eq. (2.11). In this case, we get a ~ Ap/2Ap 
which implies that « is small and we can neglect higher order terms in perturbation 
theory. With this and Eq. (5.7), we get in the laboratory frame 


2 
k. (5.10) 


E? /ħ = wak Tk(k +1) | 


Jp 


Comparing this with the case of the AC Stark shift caused by a resonator coupled to 
the qubit, as discussed in Sec.[5.2.1| we find again that the shift of the qubit frequency 
is proportional to the resonator population by the relation (âtâ} « Pp « Ap?, where 
Pp is the drive power. From this and the result of Eq. (5.1), that is Awg = (0*4) g?/ Ar, 
we additionally see that the drive amplitude of a coupled resonator is given by 
Ap? = 4$? (ata) in the limit of large detuning. 


A second conclusion can be drawn from the unperturbed eigenenergies given in 
Eq. 6-7), where we see that a only acts on terms with k and not with i2. This means 
that all levels shift in parallel, as long as « is small and the higher order perturbation 
terms can be neglected. As said, this especially holds true for the dispersive limit. 
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Figure 5.2: Isolines of the displacement parameter « for different combinations of detuning and 
amplitude. The solution of Eq. is calculated for different combinations of drive amplitude Ap and 
detuning between drive and qubit frequency Ap (solid lines). By taking the parameters relative to the 
anharmonicity y, this parameter can be eliminated. The dashed lines represent the approximate solution 
& ~% Ap/2Ap, which renders a valid approach for large Ap. This plot also demonstrates the validity of the 
assumption that « is small for Ap = Ap. 


To show the range and validity of these approximations, « was calculated for different 
values of Ap and Ap, and is depicted in Fig. (5.2|as contour plot. This shows that « 
is small for large detunings Ap = Ap and in the general case of small « the linear 
approximation is justified (dashed lines). Apart from these extremum considerations, 
any further analysis needs to take into account the higher orders of perturbation theory 
and we have to calculate the full solution of Eq. (5.6). For a better comparability with 
the later acquired experimental data, we calculated the energy of the multi-photon 
transitions from the ground state |0) in the laboratory frame and displayed them for 
different values of detuning Ap and anharmonicity y in Fig. 5.3| against the drive 
amplitude. For this graph, a base qubit frequency of wq = 5.0GHz was assumed, 
similar to the later experiments. The solid lines in Fig.|5.3]a) and b) correspond to a small 
detuning Ap, while the dashed lines represent the analytic data for an increased Ap. 
This demonstrates clearly that the bending of the individual qudit transitions does not 
only depend on Ap/24Ap, as suggested by Eq. (5.10), but on both values independently. 


With these findings, the groundwork for using an anharmonic circuit as sensor for 
detuned microwave signals is set, giving us the opportunity to detect the amplitude 
and frequency of the signal by measuring the shift of the qudit's transition frequencies. 
A comparatively small anharmonicity y, which is usually considered to be unwanted 
and a drawback of devices like the transmon, can be advantageous for this sensing 
scheme: As a depends on the relations Ap/y and Ap/Y, the quantity 1/y can 
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Figure 5.3: Analytically calculated level shifts. The multi-photon transition frequencies are calculated 
using the analytic formula with perturbation theory up to fourth order. It can clearly be seen that for an 
increased detuning (dashed lines) the level shift is less pronounced and rather parallel, whereas for a 
small detuning (solid lines) the levels closer to the drive shift more strongly and eventually cross. Part a) 
represents a weakly anharmonic circuit with an anharmonicity of y = 400 MHz, comparable to a standard 
transmon qudit, whereas part b) demonstrates that the sensing scheme is in principle also applicable to 
other qudit systems with higher anharmonicity. Adapted from (71). 


be seen as an amplification factor for a and hence increases the sensitivity of the 
sensor. This also becomes clear when comparing Figs. a) and b), where y is 
increased by a factor of 3 from a) to b). 


5.3.2 Numerical solution 


In the previous analytical treatment of the driven qudit system, several approximations 
have been applied, leaving room for an inadequate description of our experimental 
system. To get a more exact simulation of our physical system, we can perform 
a numerical analysis, which can also be executed for more complex systems that 
are hard to describe analytically. 


To resolve the level transitions of our anharmonic system under a drive and to do 
the simulations in closest analogy to the experimental conditions, we simulate not 
only the influence of the drive tone, but also of the weaker probe tone which we use 
to excite the level transitions. The simulation Hamiltonian then reads: 


Hnum = j B li) (i| + Aph(B + Bt) coswpt + Aph(B + Bt) cos pt, 
i 


where Ap and wp are the probe tone amplitude and frequency, respectively. In contrast 
to the previously used anharmonic approximation for the qudit, we now leave the 
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individual level energies E; completely unrestricted and take the anharmonic raising 
and lowering operators ĝt and f instead of the harmonic ones. To calculate the 
anharmonic operators, we have to consider the cross-coupling-constants g;j & (i|N]j) 
in the basis of the initial transmon Hamiltonian Eq. (2.26), which create additional 
off-diagonal terms in the operators ß = Y ys) gu |k) TA 


One limitation of numeric simulations is that the Hilbert space of the problem has to be 
represented by a matrix. This means that only a limited set of degrees of freedom can 
be represented in the simulation and every additional degree strongly increases the 
computation time. To this end, the Hilbert space of the charge operator N was restricted 
to ((N) = -15...+ 15) and hence 31 levels, and the transmon excitations are considered 
up to the |9) level, as higher levels are beyond the Josephson potential barrier. 


Master Equation and Lindblad Formalism 


To calculate the time evolution of the system described by this Hamiltonian, we 
use the master equation in the Lindblad form, which describes the evolution of a 


mixed state with density matrix 6 by [155,156]: 
ô i ^ " ^ Tf 1 ats x 
di = E rn] + m (Lut; 2 (E10) 4 (5.11) 


with {a,b} = ab + ba being the anticommutator. The Lindblad collapse operators 
Lai represent the different channels through which decoherence acts on the 
quantum system at a rate 7,. Relevant for our analysis are here the qudit decay 
operator Lı = B with the energy relaxation rate yı = I, and the dephasing operator 
L, = B'B/A/2 with a rate of y» = Ig. 


The third process we can take into account is the spontaneous excitation by thermal 
energy with the operator Î; = pt and a rate of 73 = Ife, Where ñtherm is the 
equilibrium thermal population of the system. If we take into account a nonzero 
thermal population, we also have to correct the decay rate to y1 = T] (Atherm + 1) [71 : 
For our simulations, we assumed a thermal population of fttherm = 0.1, corresponding 
to an electronic chip temperature of T — 100 mK. This not only accounts for possible 
dissipation caused by the strong drive tone but is also in agreement with experimental 


findings from other groups [157]. 


Simulation software and environment 


Knowing all input parameters of Eq. (5.11), we can numerically calculate the time 
evolution of a starting state 6 iteratively for small time steps dt. After a steady state is 
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reached, we calculate the expectation value of the number operator (BB) = Y k |k) (k| 
for the final state. This corresponds to our experimentally observed quantity, the 
shift of the readout resonator's frequency. 


To calculate the time evolution, the software package QuTIP was used, 
providing a Quantum Toolbox in Python. This toolbox comes, amongst other things, 
with predefined functions and operators to easily create a representation of pure 
states, density matrices, operators, and more and allows to calculate expectation 
values and matrix elements in a very accessible manner. It also features a master 
equation solver which takes finum, ĝo, Èa and a list of time steps dt as inputs and 
calculates the evolution of 6 for these time steps. 


For each combination of Ap, wp and wp, a full time evolution of the master equation has 
to be calculated. This is a computation intensive task for a detailed set of combinations, 
even with the subsequently presented method to decrease the computation time. As the 
time evolution depends strongly on the choice of parameters, no intermediate results 
can be exploited for other calculations and each combination has to be calculated 
independently. On the other hand, this means that the computation task can be 
massively parallelized, and all combinations could be computed at the same time 
- if the computation power and memory is available. 


To this end, we made use of BWUNICLUSTER in the framework of BwHPC, the high 
performance computing cluster of the state of Baden-Württemberg. This cluster 
features 512 computing nodes] each equipped with 16 cores operating at 2.6GHz and 
a main memory of 64 GB. Being a scientific computation facility, the operating system 
of the cluster is already equipped with QuTIP, providing easy access to simulations. 
Since the parallel computations do not need to be done simultaneously nor interact 
with each other, it was possible to split the computation task into several small jobs 
which can run on spare computing capacities and hence need less time in the job queue. 


Coherence and computation time 


In the experiments, we do spectroscopic measurements, meaning that we apply all 
microwave tones for several milliseconds and are certain that a steady state has 
evolved. The measurement time per parameter combination is here typically given 
by hardware limitations and the data acquisition time. For the simulations however, 
every time interval dt by which we prolong our analysis takes additional computation 
time, meaning that we try to reach the steady state as fast as possible. 


3 In total, more than 890 computation nodes contribute to the cluster. The additional nodes are however 
specialized to certain kinds of jobs and therefore not considered here. 
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Figure 5.4: Numerically calculated qudit spectrum for different collapse rates. Simulation of the qudit 
spectrum with a drive amplitude Ap/27 = 0.3 GHz at wp/27 = 4.95 GHz. To speed up the convergence 
of the master equation simulation, the decay and dephasing rates I1 = 1/T; and I, = 1/T, are modified, 
resulting in a shorter overall simulation time Tsim until the solution was in a steady state. As a good 
compromise between simulation time and visibility of all levels, the last combination (black crosses) with 
T; = 0.25 us, Tọ = 0.25 ps and Tsim = 0.5 ps was chosen for all other simulations. The lines are guide to 
the eye and do not represent simulated data. Adapted from [71]. 


The steady state under a drive is brought about by the Lindblad collapse operators 
and the corresponding Y, rates. For an all-spectroscopic steady state analysis it 
hence suggests itself to chose higher I1, Tọ rates for the simulation compared to the 
measured rates in order to reduce the overall simulation time. 


To verify that this modification does not influence the statement of our analysis, a 
spectrum of qudit transitions was simulated at fixed input parameters for different 
I1, Ig rates and an adapted overall simulation length. The results of this analysis 
are shown in Fig. where for each relaxation rate combination three peaks are 
visible, corresponding to the single photon w12, wo, and two photon w transition. 
It is apparent that for greater I; rates and hence reduced Ty times the lines get 
slightly broader, whereas a very sharp and high peak can be observed for the longest 
dephasing and simulation time (gray). The key point is however that the position 
of the peaks and their general shape does not depend on the chosen rates, leading 
to the conclusion that it is a valid approach to increase the rates for a spectroscopic 
analysis and balance between computation time and feature visibility. 


To cancel out small oscillations in (ß'ß) which are caused by the drive even in 
a steady state, the result is always averaged over a period of T,/4 at the end 
of the simulation time. 
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5.4 Experimental observation of the 
multi-mode AC Stark shift 


For the experimental measurements of the AC Stark effect, the spectroscopy setup 
with additional drive line was used, as described in Sec. To be able to study 
the influence of the drive also for stronger drive amplitudes, a high power amplifier 
was used after the drive microwave source. Since the output power of the amplifier is 
limited, the amplification ratio decreases at high powers. In a separate measurement, 
we calibrated the output power of the amplifier depending on the power setting of 
the microwave source to be able to correct our data. In the following, all specified 
drive powers correspond to the power at the output of the amplifier. 


As qudit sample, the single-junction transmon chip with capacitance islands made from 
TiN was used, described in Sec. [4.3] There, we found a fundamental qubit frequency 
of Wo, = 4.755 GHz together with an anharmonicity a) = —y/2 = -190MHz 
from spectroscopic measurements. 


To see the effect of a drive tone on multiple qudit transition frequencies, we apply 
a drive tone at a fixed frequency of wp = 4.95GHz and increasing values for the 
power and hence the amplitude Ap. This results in a detuning of Ap = 195 MHz. To 
probe the transition frequencies, we apply an additional probe tone with a power Pp, 
strong enough to excite multi-photon level transitions, and scan its frequency to excite 
different transitions. The result is then measured by a weak readout tone applied to 
the resonator frequency wr with a power in the single photon regime. 


The dispersive shift and with that the qudit population is depicted in Fig. a) 
together with the analytical solutions from Eq. (5.6). We see a good agreement for the 
first three level transitions of the qudit (solid lines) and get data as expected from 
Fig. where the wo; and wo2/2 transitions approach each other. 


However, a variety of additional lines appears in the spectrum, which require further 
argumentation. These lines can be either explained by a thermal population of the |1) 
level of the qudit, resulting in the w12 and w23/2 transitions (dashed lines). Or we 
have to take into account a multi photon process with virtual energy levels, where 
the two photons originate not only from the probe tone, but at least one of them is 
supplied by the drive. This means we can drive the |0) — |2) transition with the two 
photons wp + wp = wo (blue dotted line). Due to the fact that the drive frequency is 
fixed, this curve has a steeper slope compared to the wọo2/2 transition (blue solid line). 
With combinations of multiple photons from drive and probe tone, we are eventually 
able to describe all observed lines in the spectrum (see legend). 
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Although the measured spectrum appears to be well represented by the analytic 
formula, we see deviations starting for Ap > 0.4GHz, corresponding to Ap > y and 
hence an increased parameter « in our analytical treatment. We note here that the 
proportionality constant in the relation Ap « Pp cannot be found experimentally 
(besides with the AC Stark effect) and therefore was chosen to give a best fit between 
experiment and numerical simulations, where Ap is directly provided. 


To see whether these deviations come from the approximations made for the analytic 
solution or from experimental imperfections, like wrongly calibrated amplifiers or other 
elements in the setup that start to get nonlinear at high powers, we simulated the exper- 
iment with numerical methods as described before. The results of these measurements 
are given in Fig.5.5]b), again with the analytically calculated lines for comparison. 


Comparing the numerical data to the experimentally found values, we see a per- 
fect agreement, representing all lines found in the experiment and following the 
same deviations from the analytic solution. We can therefore conclude that our 
experimental data are correct and the deviations are caused by the approxima- 
tions in the analytical solution. 


In the numerical data, we see additional transitions, which are not visible experi- 
mentally. These transitions include the higher excited |4) and |5) states, meaning that 
the amplitude of the probe tone was chosen too high compared to the experiment. 
We also see a general increase in background qudit population with increasing drive 
amplitude, suggesting that the drive starts to populate the qudit. This background was 
calibrated away for the experimental data to increase the visibility of the features. 


5.5 Sensing scheme 


After having demonstrated that the experimental outcome matches to the analytical 
calculation and even better to the numerical simulation, we are able to invert the 
scheme and use the qudit as a sensor for off-resonantly applied signals by measuring 
the shift of the qudit transition frequencies and deducing the parameters leading to 
this shift. This means we can measure the amplitude of any drive field that couples 
to the qudit by means of the qudit itself, leading to a new in-situ measurement 
method. In doing so, the qudit can be under the exact same conditions as it will 
be used for other experiments later. 
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Figure 5.5: Measured and numerically simulated AC Stark shift together with analytical transitions. 
a) The measured transmon level transitions for a drive at wp = 4.95 GHz and increasing drive amplitude 
Ap. The analytic solution of Eq. is displayed as colored lines. Especially for the level transitions 
(solid lines), a good fit to the measured data is observed. Other transitions including thermally excited 
states or drive photons can also be modeled by the analytic solution (dashed /dotted lines). Deviations for 
higher drive amplitudes are explained by the use of a simplified transmon Hamiltonian in the deduction 
of the analytic solution. The data is column-wise normalized as the drive also influences the resonator 
frequency. b) Numeric simulations in close analogy to the experiment. The data is well comparable to 
the measured data and shows the same deviations from the analytical solution, since the full transmon 
Hamiltonian is used for the simulation. Due to a slightly higher probe power in the simulation compared 
to the experiment, more transitions are visible which include the |4) state of the qudit. Adapted from nl. 
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Figure 5.6: Graphical representation of the FLT. From the measured data on qudit base frequency and 
anharmonicity, each in reference to the unperturbed qudit parameters, we can pick a point on the 
frequency lookup table (FLT) and resolve the combination of drive amplitude Ap and frequency wp which 
leads to this shift. While this FLT was calculated using the analytic formula in Eq. (5-6), it is also possible 
to fill the table with data from numeric simulations, resulting in more accurate results. Adapted from [n]. 


5.5.1 Frequency lookup table 


To operate the qudit as a sensor, we have to calculate back from the measured 
frequency shift to the originating drive field. This is impossible to do for the numerical 
simulations and increasingly difficult when taking into account the higher order terms 
in perturbation theory. Instead of iteratively approaching and fitting the parameters 
for each time we do a measurement, it is more convenient to calculate once the 
transition frequencies for a set of different combinations using the unperturbed qudit 
parameters. We store these transition frequencies in a frequency lookup table (FLT) 
and only have to find the closest matches in the table when doing a measurement. 


For a graphical representation of the FLT, we chose to plot the table contents 
against a reduced anharmonicity y* = (2w 91 — wo2)y and the normalized qubit 
frequency «B /w(,, where wm is the measured qudit frequency and «0, is the 
unperturbed frequency. Plotting into this coordinate system lines of constant Ap and 
Ap, respectively, brings us to a representation like in Fig. 5.6] With a measured y* and 
w® / w0, we can then directly read the values of Ap and Ap responsible for this shift. 
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5.5.2 Measuring the transmission function 


Putting the sensing scheme to the test, we use it to measure the transmission function 
T(w) of our setup. This includes the whole way from the microwave source via the 
cryostat wiring to the sample box and even inside the box from the transmission line to 
the qudit. Being able to measure only parts of it when warmed up to room temperature, 
the whole transmission function is usually completely inaccessible. It has to be noted 
here that the proportionality constant between drive amplitude at the chip and applied 
power at room temperature is now frequency dependent with Ap « T(wp)Pp. 


For the measurement, we fix the applied drive power Pp and scan the drive frequency 
Wp in a second experiment. Scanning the probe tone, we can identify the qudit 
transitions and hence calculate y* and Wot, which we use to look up the drive 
amplitude Ap. Since the frequency wp of the drive is known, we can either search 
for the best matching drive amplitude for the parameter combination (4^, wij, Wp) 
or we leave the drive frequency as open parameter as well. It is of course expected 
that the results for the case of a given wp are more accurate, but the sensing 
scheme should nevertheless be able to determine the frequency of the drive tone. 
For increasing detuning, the influence of the AC Stark shift decreases, resulting 
in a decreasing resolution and accuracy. 


The results of this measurement are given in Fig. where the raw experimental 
data, i.e., the dispersive shift at different probe powers and drive frequencies is given 
in subfigure c). The precision of the measurement is indicated by the shaded green 
and gray areas in Fig. [5.7]a) for the amplitude and Fig. 5.7]b) for the drive, which 
visualize the spread in results when varying the fitted «jj and wo2/2 values by 
Aw/27 = X1MHz. This is a worst case estimate based on the rather coarse step 


size when scanning the probe tone, Awp/27 = 1 MHz. For the results, this gives on 
average AAp = +92 MHz and Awp = +124 MHz, if wp is taken as a free parameter. 
As the sensitivity for the drive frequency and therefore the resolution in this parameter 


gets very poor for a detuning of more than 5 GHz, no errorbars are shown in this 
regime in Fig. |5.7|b) and the extracted points are plotted in gray. For the case of a 
fixed wp, the error in amplitude reduces to AAp = +12.5 MHz. 


This error estimation is very crude and also the measurements themselves have 
to be seen more as a proof-of-principle demonstration than as a highly optimized 
measurement. The resolution could be greatly improved by different means which will 
be discussed in the next section. Nevertheless, the principal shape of Ap(wp) makes 
clear that this method can be used to make statements about the transmission function 
of the cryostat. For an ideal setup at constant input power, we would expect an increase 
of Ap close to the qudit transition frequencies due to the frequency-dependent dipole 
moment of the qudit and for higher frequencies a monotonous decrease converging 
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Figure 5.7: Proof-of-principle demonstration of the sensing scheme. For a constant drive power Pp at 
room temperature, the drive frequency is scanned and the qudit transition frequencies are recorded. With 
the help of a FLT, we can calculate back to the amplitude a) and frequency b) of the drive. As explained in 
the main text, the accuracy of the amplitude detection can be greatly improved by taking the known 
drive frequency as additional input (green circles). For the blue crosses, both Ap and wp are free. The 
shaded areas illustrate the uncertainty of the sensor when taking into account an accuracy of +1 MHz for 
both wo /27t and wo2/47. For more than 0.5 GHz detuning between drive and first qudit transition, the 
frequency determination fails in this demonstration, therefore the data points are grayed out and no error 
bars are shown. Subfigure c) shows the measurement data, i.e., the dispersive shift of the resonator in 
gray scale. To verify the validity of taking wp as input parameter, the transition frequencies are calculated 
from the data shown in a) and b) and plotted as blue crosses (green circles) when having wp as free 
output parameter (fixed input parameter). For both methods, we find a good agreement to the measured 
data within the resolution of the measured data. Adapted from (71). 
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Figure 5.8: Demonstration of increased sensitivity and resolution. By taking Ramsey measurements 
instead of spectroscopic measurements as in Fig. [5.7] the resolution in measuring the qudit transition 
frequencies is improved, leading to a greatly increased accuracy of the sensor. In addition, the FLT was 
generated by numeric simulations and not by the analytical approximation as before. These measurements 
are taken on a different chip and in a different cryostat than the ones in Fig. so the values are 
not directly comparable. Also, due to the fact that no input amplifier was used for the drive, the drive 
amplitude is much lower. With reduced drive power and increased minimal detuning, the shifts in woı 
and wo» /2 are much smaller and barely visible in subfigure c). The overall accuracy in both amplitude 
and frequency is however greatly increased, visualized by the much smaller gray band around the data 


points. Courtesy of M. Kristen [72] [160]. 
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towards a constant value. The measurement however only shows the onset of this 
increase towards wo, and oscillations of Ap for higher detuning. These oscillations can 
be explained by a non-ideal transmission of the microwave lines through the cryostat, 
being soldered and connected at each temperature stage of the cryostat. 


5.6 Further advances 


After the publication of the paper associated with this chapter [AL demonstrating a 
proof-of-principle experiment, the method was further studied and improved, in par- 
ticular by the associated Master's thesis of Maximilian Kristen which also resulted 
in a joint paper [72]. The results of this work are briefly summed up in the following. 


5.6.1 Drawbacks of approximate theory 


To increase the accuracy of the sensor, it is important to create a reliable FLT. With 
the deviations seen in Fig. P-7]it is reasonable to assume that the analytic solution 
is not perfectly suitable for the calculation of a good FLT. The origin of these 
deviations was found by performing numerical simulations with the approximated 
transmon Hamiltonian Eq. (5.3), resulting in a very good agreement between the 
analytic solution and the simulation result. This leads to the conclusion that the 
approximation of the transmon Hamiltonian is too inaccurate and dominates the 
discrepancies to the experimental results, whereas the fourth order perturbation 
theory can be seen as a valid approximation. 


To this end, the FLT was calculated by numerical simulations, taking more computation 
time but finally giving more accurate sensing results. 


5.6.2 Time-domain approach 


As an alternative to spectroscopic measurements which require fine-grained frequency 
sweeps, pulsed measurements can be used. Under certain conditions they allow for a 
very fast determination of the qubit transition frequency by utilizing an optimized 
detection scheme [161]. This scheme is based on Ramsey oscillations explained in 
Sec. 4.2.3land can use advanced methods to minimize the number of time steps At 
necessary to get a value for wg]. The same technique can be applied to measure Ramsey 
oscillations between the |1) and |2) states and hence infer the w17 frequency. 


These measurements have been performed on a different qudit sample in a different 
cryostat, so the results, shown in Fig. are not directly comparable to the ones 
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previously found. However, it is obvious that the sensing accuracy is now greatly 
improved based on an enhanced FLT and finer resolved measurement data. With the 
improved FLT, the data for the extracted drive frequency now matches the applied 
drive frequency much better, especially in the regime of low Ap, where the analytically 
calculated FLT has major inadequacies. With an average error of Awo, = +6.7 kHz 


and Awı2 = +10.5 kHz, the accuracy of the input data was also improved by about 
two orders of magnitude compared to the first demonstration experiment. 


5.7 Conclusion 


In summary, we have demonstrated a sensing scheme for microwave tones by a 
superconducting anharmonic multi-level system, based on the AC Stark effect. This 
sensing scheme can be of great use for the optimization and engineering of quantum 
circuits as it enables to measure the cross-talk of microwave structures to a qudit on the 
exact same chip as intended for other quantum experiments. Using directly the qudit 
as sensor eliminates the need to fabricate additional chips for pre-characterization of 
spurious couplings and helps to validate the results of microwave simulations. 


With the possibility to measure and characterize oscillating signals in frequency 
and amplitude, this sensor also lends itself to the field of fundamental research, 
where it can be used to detect resonant frequencies of other systems in a quantum 
environment. As semiconducting structures typically stop working at cryogenic 
temperatures and dissipate a lot of heat, many commonly used and commercially 
available sensors for microwaves are not usable at millikelvin temperatures. A 
simultaneous detection of amplitude and frequency exploiting the characteristic 
changes in the qudit’s transition frequency fingerprint can be seen as clear advantage, 
which outperforms other qubit-based sensing schemes. 


For all kinds of cryogenic experiments involving microwave signals, this sensor can 
be used to measure the frequency-dependent microwave transmission through the 
cryostat and the quantum chip, which is otherwise not accessible at the operation 
conditions. We are able to characterize the transmission for a broad frequency range 
of more than 1 GHz above the qudit transition frequency, where the resolution starts 
to decrease with increasing detuning. 


For the sensing scheme, technologically simple spectroscopic measurements are 
sufficient, but pulsed measurements can be used to increase resolution and decrease 
the measurement time. By performing Ramsey experiments, we can exploit the high 
sensitivity of the qudit phase towards changes in the transition frequency. With this, 
the sensor accuracy is greatly enhanced and we demonstrated an average amplitude 
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resolution of AA/27 = 3.4 MHz as well as an average frequency uncertainty of 
Aw /27 = 25 MHz [72]. The accessible sensing range and resolution can be adapted 
in time-domain measurements by adjusting the reference frequency until one of the 
following limits is reached: For low signal amplitude and large detuning, the frequency 
shifts decrease until the phase accumulated in a Ramsey measurement within the 
decay time of the qudit states is no longer measurable. Typically, the lifetime of the |2) 
state is the limiting factor here, as the shift of this transition is smaller and the lifetime 
of this state is shorter than for the |1) state. In the other limit of low detuning and 
strong drives, resonant effects like level transitions and the Autler-Townes splitting 
have to be taken into account, which can finally disturb the sensing protocol. 


Finally, our theoretical description can be seen as an extension of the Jaynes-Cummings- 
Hamiltonian for multi-level systems which is not restricted to the dispersive regime 
but can be applied also for low detunings. 
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in magnetic fields 


Using superconducting transmon qubits for hybrid systems and as sensors in magnetic 
applications requires a proper qubit function under the influence of magnetic fields. In 
the following chapter, I present our investigations on a transmon qubit in a magnetic 
field to fathom these limits of application. After an introduction covering existing 
research in this field, I will present spectroscopic measurements, which reveal addi- 
tional properties of the qubit. To characterize the losses of oscillating superconducting 
systems, the quality factor of the resonator is studied, directly followed by time 
domain measurements on the qubit. Here, different loss mechanisms are explained and 
evaluated, before the chapter concludes with measurements on the pure dephasing 
rate. The findings presented in this chapter lead to the publication of Ref. [70]. 


6.1 Introduction 


To use a superconducting qubit for quantum sensing applications and for hybrid 
systems, it is of great use to know its limits of application. Especially the critical 
parameters of superconductivity are here of interest, as they disturb the qubit's 
working scheme. These are in particular the critical current, critical temperature 
and the critical magnetic field. 


Limitations in terms of current can be seen by either inducing a constant current 
through a SQUID loop, as it was shown for the tunable sample, or by creating an AC 
electric field across the JJ, which was in depth explained in the last chapter in terms of 
the AC Stark shift. Both methods lead to a controllable decrease in qubit transition 
frequency, which has been studied and used in application [17,71]. 


The behavior of a qubit at elevated temperatures is known to be limited by a strongly 
increased number of quasi-particles (QPs) and long before the eventual breakdown 
of superconductivity [162]. Both relaxation and coherence times have been observed 
to rapidly decrease above a temperature of 100mK [163] 164]. 
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Besides that, the influence of high frequency phonons has been studied by the 
application of ultrasonic pulses to superconducting quantum bits [165], as well as 
the dependence on stray infrared light [166]. 


For magnetic fields however, no in-depth study was made which investigates the 
coherence properties of a superconducting qubit under increased magnetic fields, 
until the coherence breaks down. So far, only the influence of a magnetic field up 
to the order of 100uT was studied, where a slight improvement of the coherence 
properties was found due to the creation and pinning of flux vortices which act as 
OP traps [113]. Other research groups assume superconducting qubits to break down 
at very little fields, although they do not provide any measured data to support 
this claim. Instead, they propose other, more stable junctions, like the proximitized 
semiconducting nanowire [167]. 


The general consent is however to screen magnetic fields as good as possible and 
meanwhile, a multilayered shielding based on permalloy and superconductors is 


considered to be state of the art [168; 169]. 


There are many applications like quantum sensing of magnetic excitations [29], 
the creation and harnessing of Majorana fermions or research fields like 
cavity quantum magnonics [64] [66], which all require the qubit to be in or close to 
magnetic fields. It is therefore inevitable to study the influence of magnetic fields to 
superconducting qubits. For all these applications, it would be beneficial to maximize 
the magnetic field in which the qubit can be operated. As we know from Sec. 
we should then align the magnetic field in parallel to the chip surface. At this point 
it should be however noted that a parallel field orientation over the whole chip 
is hard to achieve, because the structured superconductor on the chip bends the 
field lines and leads to a flux focusing for the regions free from superconducting 
material. For specific applications, the inclusion of structured magnetic materials will 
also heavily influence the local magnetic field distribution. Therefore, no additional 
effort was spent to assure a parallel alignment but the chip was straightened in 
the solenoid by eye, as described in Sec. 


6.2 Qubittransition frequency in an in-plane 
magnetic field 


Before we perform pulsed measurements, we firstly have to measure the dependence 
of the qubit transition frequency on the magnetic field. As described in Sec. the 
critical current of a JJ in an in-plane magnetic field periodically goes down, where the 
shape follows a Fraunhofer pattern. With wo, « VIe, a similar behavior is expected 
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Figure 6.1: Measured transmon frequencies in a magnetic field. The dispersive shift of the resonator is 
shown color-coded to visualize the transmon transition frequencies under the influence of a magnetic 
field. For each field value, the readout frequency was adjusted to follow the resonator frequency. The 
absolute value of the normalized dispersive shift is plotted here to account for these offsets. In the original 
data, the dispersive shift changes its sign due to a changing quality factor of the resonator, which would 
distract the colorscale. According to our expectations, we see a decrease in frequency and a periodic 
revival. Compared to the expected Fraunhofer pattern of a single JJ in a magnetic field however, the top at 
B — 0 is too flat and the side maxima are too high in frequency. Adapted from [70]. 


is expected when placing the qubit junction in an in-plane magnetic field, such that 
flux can enter the insulating barrier. This would include a pronounced maximum at 
o / dy and periodic side maxima with decreasing amplitude. 


6.2.1 Measured data 


When we place the TiN qubit introduced in Sec. (4.3|inside the solenoid and apply 
a magnetic field, we observe a periodic decrease in wo; following this behavior, 
as shown in Fig. If we compare this data to the expectation of a Fraunhofer 
interference pattern, the measurement data shows a flat top at  — 0, very steep 
drops at ® ~ dy, and side maxima which are almost as high as the main maximum. 
The pattern can hence not be explained by a simple rectangular junction cross 
section. Since the junction was fabricated in a shadow angle evaporation process, 
deviations from a perfectly rectangular cross section are very likely. Therefore, we 
Fourier transformed different possible junction geometries but did not find any I.(B) 
respectively «9; (B) relation resembling the measured data. A different process is 
hence needed to explain the measured results. 
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Figure 6.2: Illustration of all created JJs. Following the fabrication scheme in Fig. we find two 
additional JJs being created. While JJ3 only shunts the lower Al; layer and is therefore negligible for low 
kinetic inductance films, the qubit excitation current flows through JJ which explains the additional 
contribution. The current path is indicated by the white arrows. Adapted from [70]. 


6.2.2 Revisiting junction fabrication 


If we revisit the junction fabrication process for the TiN sample from Sec. 4.4.1] 
we see that besides the intended small Josephson junction JJı in the center of the 
structure, two other junctions, JJ? and JJ3 exist, which are commonly neglected due 
to their large junction area and hence high critical current. These junctions are 
illustrated in Fig. Due to the multi-step fabrication process in which first the 
TiN capacitance pads are fabricated and then the junctions on top of this, the current 
has to flow across JJ; and JJ; to go from one electrode to the second, as indicated 
by the white arrows. JJ; only shunts the lower Al layer, which is a negligible effect 
as long as the inductance of this layer is low. 


If we qualitatively analyze the influence of the parasitic junction JJ}, which is in series to 
the intended qubit junction JJı, we see the following behavior: As long as the parasitic 
junction is not close to being blocked by a flux quantum, its critical current is very high 
and not limiting or influencing the qubit current The qubit transition frequency can be 
therefore assumed to be rather flat and unimpressed by the magnetic field. In the case 
of $/ Dy z n (n € N \ 0) for the parasitic junction, we will see a steep decrease of the 
qubit frequency. This goes along well with our measured data, and we will deduce a 
quantitative formula for a transmon qubit with two junctions in series in the following. 


Transmon Hamiltonian for two JJs in series 


To derive the transmon Hamiltonian for two JJs in series, we start with the general 
Hamiltonian introduced in Sec. which has contributions from the charging energy 
in Eq. (2.24) and the Josephson energy in Eq. (2.25). In the general case, we have: 


H-x TEOL (6.1) 
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To solve this, we have to find the current-phase-relation I(q) for the total system. 
We see that the phase acquired by the individual junctions, 91 and 2, adds up toa 
total phase, and the current through the junctions has to be the same: 


9x = 91+ G2 
I(gx) = Ia sin 91 + 2 sin qa. 


We introduce the ratio r = I,2/ I< between the junction’s critical currents to facilitate 
the writing, and chose r > 1 without loss of generality. This gives 


sin 9x 


tan g2 = —————— 
P2 r +cos px 


and we get for the overall current-phase relation 


I(gy) = Ic» sin (arctan nee) à (6.2) 


The integral of this equation can be written in a closed form and we find 


r? + 2rcosgy +1 
J (oxides > ee : = = -La yr? + 2rcos Qx +1. 
This results in a total system Hamiltonian of 
H = 4EcN? — Eja "E +2rcos gy +1, (6.3) 


where we have again used the charging energy Ec = e?/2C and the Josephson energy 
of the qubit junction Ej1 = I,,,;o/27t, and number and phase operators as in Sec.[.5] 


In analogy to Eq. (2.27) we perform a Taylor expansion up to fourth order in gy and find 


r 2 r(r-r+H1) 4 
7422 Zar ri) Pr er 


with 1 JA pt Í 


H = 4Ec Ñ? + Ej ( 


With the same approximations as for the simple transmon case, this gives 


T v^ r? —r +1 EC 3 
H=,/8EcE bb bti (bti +1). : 


In the limit of r = 1;5/1;1 — oco, where JJ; is dominating (i.e., limiting) the qubit’s 
behavior, this formula recovers the simple approximated transmon Hamiltonian from 
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Figure 6.3: Two-junction transmon. The analytic transition frequencies for a transmon qubit with two JJs 
in series are shown on top of the measured data. Here, JJ; gives the overall envelope and JJ? creates the 
periodic minima. For the measured data, the same normalization was applied as in Fig. [6.1] The analytic 
formula matches the measured data very well for wo, and w12, which appears in the measurement due to 
thermal population of the |1) state. Adapted from [7o]. 


Eq. (2.29). We also see from this formula that the anharmonicity changes with - sit 
which means that the effective anharmonicity of a transmon with two junctions in series 
is always lower than if only one junction was present. For the limit r = 1, where the two 


junctions in series have the same critical current, the anharmonicity is reduced by 75 %. 


When we assume an almost rectangular cross-section for the two junctions, each 

sinc I ^|. Besides the periodicity field Bo, 
we also introduced an offset field Ba here, which accounts for stray background 
magnetic fields. Fitting Eq. to our data, we find a good agreement for wo 
and w12, shown by the solid and dashed red lines in Fig. The extracted values 
for each junction are given in Table 


critical current should follow I. = I? 


From the periodicity of the pattern, we can calculate back to the junction's size with 
Po = Ba, A. The effective cross section of the junction is given by A = (d+2Ar)! with 
l being the junction length. Since the magnetic field reaches into the superconductor 


Table 6.1: Parameters of the two Josephson junctions. B4 is the offset field, Bo, the periodicity, and ! the 
length of the junction. 


Ej/h Ba Bo 1 


0 


Jı 1615GHz 18mT 300mT 209nm 
Jj  300GHz -0.2mT 25.5mT 246gum 
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on a length scale of AL above and below the junction, this has to be added to the 
thickness d, which is the oxide barrier. For our junctions, we take d = Inm and 
AL = 16nm for aluminum. With this we get the junctions’ lengths to be l4 = 209 nm 
and l) = 2.46 um, which agrees well with the design parameters. 


The increased sensitivity against magnetic fields and the abrupt drops in the qubit 
transition spectrum can be disadvantageous for many applications. To avoid this 
spurious influence, the parasitic junctions JJ? and JJ5 can be shunted with a bandage 
in an additional fabrication step [171]. Alternatively, other fabrication methods apart 
from shadow angle evaporation can be used, e.g., overlap junctions [172]. In this 
process, no additional junctions are created. 


Due to the fact that we cannot see the first minimum of the qubit junction in our 
measured data, not to mention further maxima, the extracted size of the qubit junction 
has to be taken carefully. What has been neglected in the treatment up to now is 
that the critical current of a Josephson junction depends also on the superconducting 
gap of the leads, which in turn depends on the magnetic field. This is modeled 


by the Ambegaokar-Baratoff relation [173, 174]: 


A(B) 


7T 
IcRn = —A(B) tanh SkaT 


(6.6) 


For our values of A and T, the tanh term is negligible and we see that I. « A. Together 
with A(B) = Ag4/1 — (B/ B.)?, where B, is the critical field, we get a relation for 
Icı(B) and can calculate the qubit transitions. Assuming a critical field of Be = 168 mT 
which gives the best fit to our data, we see no difference to the previously presented 
model within the range of our acquired data. This can be seen by the good match 
of the colored dashed lines in Fig. In this limit, we assumed that there is no 
Fraunhofer-type reduction of I, ;, i.e., the qubit junction is assumed to be point-like. 
In reality, both effects coexist at the same time. 


6.3 Resonator behavior in magnetic fields 


Before we come to time resolved data on the qubit coherence, we shortly take a look 
at a similar measure: the quality factor of the readout resonator. Since the transmon 
can be explained as a modified resonator, it seems likely that the transmon is affected 
by the same spurious mechanisms as the resonator and the resonator's quality factor 
Q = w/ ^o evolves similar to the qubit decoherence time with T = 2/A«q, where 
Aw, is the qubit linewidth. At the same time, the resonator acts as readout device 
and eventually limits the signal-to-noise ratio of our measurement within reasonable 
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Figure 6.4: Comparison between different I. (B) relations for JJı. The critical current of JJ; creates the 
envelope for the transition frequencies. It can be either suppressed by a Fraunhofer-type interference 
pattern as explained in Sec. [2.4.2] (red dashed line) or by a reduction of the superconducting gap energy A 
(blue dashed line). The narrow-dashed lines show the wo transitions, and the wide-dashed lines show 
w12. For the measured data in the background, the same normalizations apply as before. Within the range 
of acquired data, no difference between the to methods can be seen. Increasing the measurement range is 
also not possible with the used setup due to the strong reduction of the readout resonator's quality factor 
and therefore decreased signal-to-noise ratio. Adapted from [70]. 


acquisition time. This is already visible in Fig.(6.1| where the qubit signal fades out for 
high field values due to the decreasing quality factor of the resonator. 


Superconducting resonators in magnetic fields have already been studied in detail, 
with a focus on an alignment of the magnetic field perpendicular to the resonator 


[130 175], QP recombination after infrared pulses [116] or a comparison between 
zero-field-cooled and field-cooled data [115]. 


Our measured data are presented in Fig.[6.5|and match well to the previously published 
data. The data are extracted by means of a circle fit [136], where applicable. For very 
low quality factors, the circle fit did not converge due to a small phase contribution 
of cable resonances which could not be calibrated away. For these measurements, a 
Lorentzian fit to the amplitude was used to obtain the resonance frequency, which 
is shown as open symbols in Fig. 


The data correspond well to the in-plane data previously measured by other groups 
[131], where they compare it to a loss rate calculated from the classical Bean model 
133]. This simulation matches our data well for the case of a weakly inhomogeneous 
RF current distribution. Although Fig. |6.5|c) suggests that the resonator is completely 
interspersed with flux vortices at B ~ 100 mT, a closer look in the data does not 
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Figure 6.5: Resonator data in a magnetic field. Spectroscopic data of a TiN A/2 resonator in microstrip 
geometry. As errorbars, the uncertainty of the fit values is indicated. A decrease in resonance frequency a) 
can be seen as well as a strong decrease of its internal quality factor b). Both parameters are subject to a 
strong hysteresis. Red (blue) symbols are used for an upsweep (downsweep) of the magnetic field. Filled 
symbols represent data extracted from a circle fit. Where this fit did not converge, a Lorentzian fit was 
applied (open symbols) to extract the resonance frequency. The gray bars and symbols denote the field 
range of the corresponding sweep data. c) shows the inverse of the quality factor for a better comparison 
with the data presented in Ref. 1131]. Adapted from {70}. 


support this statement, since the data extracted by a Lorentzian still shows a difference 
between up and down sweep for the loaded quality factor (data not shown). 


For the coupling quality factor Q. we see no significant change for different values of 
B, which was expected since this quantity is defined purely by the geometric coupling 
of transmission line and resonator. We find an average value of Q. = (9.3 + 1.3) MHz. 


From the resonator measurements we already see an offset field of Bog, = 8.5 mT 
in both wr(B) and Q;(B), which we will find again for the later measurements of 
the qubit coherence parameters. We assume that this offset field appears due to the 
presence of stray fields from magnetized components around the qubit chip, which 
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are partially aligned perpendicular to the chip. If we assume small misalignment 
between the chip and the coil field, we find that an applied field of 8.5 mT would 
compensate a perpendicular magnetic field of B ~ 450 pT, which is well on the order 
of typical stray fields, for a misalignment of about « z 3°. 


6.4 Qubit coherence in magnetic fields 


6.4.1 Measurement scheme 


To acquire reliable data and account for the variety of changes which appear when 
changing the magnetic field, a measurement scheme was constructed which ramps the 
current of the solenoid to the specified field value and then recalibrates the necessary 
experimental parameters to start the measurement, as follows: 


1. The readout tone is scanned around the previously measured frequency of 


the resonator. The recorded amplitude signal is fitted to a Lorentzian and the 
resonance frequency wr taken as new readout frequency to get a good contrast 
in the dispersive measurements. 


. The manipulation tone is scanned across the previous qubit transition frequency 


to find the new qubit frequency woı by again fitting a Lorentzian to the phase 
data of the recorded dispersive shift. 


. Applying a Rabi sequence at the previously determined x; frequency and 


fitting the data to an exponentially decaying cosine yields the time t,, for a 7 
pulse. 


. A sequence of Tı, TÈ and/or TF is executed to get the desired coherence 


parameters. 


In the steps 1-3, selection mechanisms apply to sort out non-converging fits by 
detecting too high errors or fit results that leave a predefined range. The sequence is 
then repeated once for the same field value and skipped if it fails the second time. The 
number of averages as well as points in frequency (steps 1 and 2) or time (steps 3 and 
4) is reduced to be able to perform the whole sequence within 10 minutes while still 
maintaining a data quality good enough for the fitting, despite the decreased SNR. 
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Figure 6.6: Time-domain measurements of a transmon in a magnetic field. Red (blue) triangles show 
data acquired on an upsweep (downsweep) of the magnetic field. The colored lines are only guide to 
the eye to trace related data points. As errorbars, the standard deviation of the fits is indicated. a) 
The measured qubit transition frequency matches well the calculated transition frequency (gray dashed 
line) and shows no hysteresis. b) The resonator frequency shows in contrast a clear hysteresis and a 
small decrease with applied field. c) A hysteresis is also seen for the qubit decay time Tı, which shows 
a pronounced maximum at low fields. d) The resonator's quality factor extracted from time domain 
measurements can only be seen as an indicative value due to the non-linear amplitude detection scheme 
in the time domain setup. Properly measured quality factors are found in Fig. [6-5] Adapted from [70]. 


6.4.2 Measured data 


The acquired data for multiple measurement sequences in the range of B = +23.7mT 
is shown in Fig. 6.6] Again, all data points acquired on an upsweep of the magnetic 
field are marked in blue, where the data from a downsweep are shown in red. In 
Fig.[6.6]a) we see again the measured qubit transition frequency, which follows Eq. 
and shows no hysteresis. This means that the qubit transition frequency is dominated 
by the characteristic Fraunhofer pattern of the two serial junctions. 
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6.4.3 Separation of loss mechanisms 


Looking at the measured data for the relaxation time Tı, presented in Fig. 
we see a prononounced maximum at low fields and a distinct hysteresis between 
up and down sweep. We can now analyze this data phenomenologically by sep- 
arating the losses of the system as 

1 

Ti =i Thyst + Thon-hyst + Teonst- 
In this classification, Ihyst accounts for loss mechanisms showing a hysteretic field 
dependence, Tnon-hyst collects losses that depend directly on the magnetic field strength, 
and the losses associated with Icons: do not depend on the magnetic field. 


As already discussed in Sec. a variety of different loss mechanisms comes into 
play, especially when applying a magnetic field. For Ihyst we assume that the main 
contribution is coming from the creation, annihilation and movement of flux vortices. 
As described in Sec. 2.6.4] it is known both theoretically and experimentally that the 
vortex configuration in a superconductor exhibits a large hysteresis on its magnetic 
history, which can be explained by the Bean model [132/133]. Due to the large aspect 
ratio of the qubit island with a diameter of 554 um and a thickness of 40 nm, we assume 
that the small misalignment to the external field suffices in order to generate flux 
vortices perpendicular to the film. This also seems to be a realistic assumption when 
regarding other experiments, where much effort is required to get a field alignment 
without the creation of perpendicular vortices. Otherwise, these vortices dominate 
the loss and frequency properties of the resonator [176]. 


From a comparison to the resonator loss data we see that shapes and signs of 
Qı and T; are generally similar (compare Figs. [6.6|c) and d)), but the qubit loss 
data feature a much more pronounced overall envelope. We attribute this envelope 
- and therefore a contribution of Inon-hyst — to the presence of excitations in the 
superconductor, i.e., quasiparticles (OPs). From the discussion in Sec.D.6.3|we know 
that the losses due to OPs are proportional to the density of QPs and hence to the 
squared magnetic field, Top « xop « B?. The QP density is not reported to have a 
hysteretic dependence on the effective field and the relaxation to an equilibrium QP 
density is expected to happen within a few ps. The hysteretic vortex configuration 
however affects the effective field in the superconductor as well as the QP creation 
and recombination rate, and with that the OP density. 


Other losses like radiative decay and losses to two-level-systems or defects in the 
dielectrics, which are discussed in Sec. are assumed to be independent from 
the magnetic field and are therefore summed up in Iconst- 
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Figure 6.7: Purcell limitation of the qubit decay time. Calculated Purcell contribution to the T time of 
the qubit, following Eq. using the measured resonator data shown in Fig. [6.5]and the calculated 
qubit transition frequency for the corresponding field values. Regular spikes in the calculated T; time 
can be seen, coming from a strongly decreased qubit transition frequency due to the blockage of JJ and 
therefore an increased detuning between resonator and qubit. The resulting T; times are at no point 
limiting the qubit coherence and the Purcell losses can therefore be neglected or considered to be constant 
against the magnetic field for this sample. Note the larger range of the x axis compared to the T; data 


presented in Fig. 6.6] 


For the Purcell loss, as introduced in Sec. a variety of strongly changing 
parameters have to be taken into account. These are for instance qubit and resonator 
frequency and resonator linewidth. The Purcell loss is therefore also subject to changes 
in the magnetic field. Taking the spectroscopic data as presented in Fig. 16.5] and 
combining it with the analytic transmon transition formula from Eq. (6.5), we find the 
Purcell loss contribution depicted in Fig. 6.7] Since the Purcell contribution ranges 
from Ipurcei ~ 350 Hz to IPurce ~% 40 Hz for fields much stronger than the ones in the 
previously presented measurements, we conclude that the Purcell contribution does 
not change the qubit’s relaxation rate significantly for any applied magnetic field value. 


Summing up the argumentation in this section, we see that we can model the qubit 
relaxation rate by I = Iconst + C(B — Botts) (113/115)/116]. If we fit this envelope to our 
measured data, we find T;onst = 53.4 kHz, C = 0.785 kHz mT 2 and Borg = 2.25 mT. 
The fit is modified here such that only solutions which are below the measured data 
are taken into account. We assume the remaining hysteretic losses to come from the 
entering and movement of flux vortices. The different offset field Bos¢, compared to 
the previously stated offset is explained by the fact that in this evaluation, the data 
at high magnetic field and therefore low decay time are weighted more strongly 
compared to the previous case, where the maximum of T, was evaluated. Although 
the parabola shown in Fig. [6.8] is a proper envelope for the measured data, we do 
not want to make any claim that this is a proof for our chosen partitioning of the 
loss mechanisms. In fact, the different loss mechanisms are not distinguishable by 
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Figure 6.8: Modeling loss contributions. The decay rate I} = 1/Tı is presented, where red (blue) triangles 
again mark upsweeps (downsweeps). The black line is a lower-limit-fit including a constant relaxation rate 
and a parabolic contribution as described in the text. The colored lines connect sequentially taken data 
points and are a guide for the eye. Adapted from [70]. 


our measurement technique. Our partitioning only represents the most obvious loss 
channels and their dependence on magnetic fields. 


6.4.4 Increased magnetic field 


If we further increase the magnetic field, we observe quantum coherence of the qubit 
up to values of Bappı © 40 mT. This measurement is shown in Fig. Although 
Ti = 0.49 us is significantly reduced at this high field, we clearly observe Rabi 
oscillations and an exponential decay after a 7t pulse, as demonstrated in the inset 
of Fig. The quality factor Qj of the resonator is significantly reduced at this 
point, which explains the low SNR. For even higher values of B, we are not able 
to track the qubit transition due to this SNR and the decreased T1, T? times which 
result in a broadening of the qubit linewidth. 


Performing a down sweep from this field value (blue triangles) does not show a 
pronounced maximum as before but only a slight increase in Tj over a broader range. 
We see however a fine structure in the data, showing multiple drops in Tı which 
coincide with the onset of a deviation from the Fraunhofer pattern, followed by a 
jump in frequency. We attribute this effect to the presence of flux vortices in the qubit 
islands due to the previously applied high fields. Their local field influences the field 
seen by the junction and therefore qubit frequency and coherence. 
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Figure 6.9: Coherence properties up to higher field values. Increasing the magnetic field further, we can 
see qubit coherence up to B — 39.6 mT and a) still observe Rabi oscillations and b) measure a relaxation 
rate. On a downsweep of the field afterwards (blue triangles in c)), the Tı time does not recover as before 
but shows multiple abrupt changes. These jumps are attributed to the movement or annihilation of a 
flux vortex. The values shown for (cz) in a) and b) are nomalized to the accessible qubit values, i.e., an 
increased residual population is calibrated away. Adapted from [70]. 


6.5 Pure dephasing rates 


Up to now, we only considered the relaxation rate of the qubit, I1. To get information 
on the phase coherence of the qubit, which is typically the quantity of interest for 
quantum computation and many quantum sensors, we simultaneously measure 
Tı and If and relate them via 


1 — 
T3 


1 
m= a + Ip, 


where Ig is the pure dephasing rate. If we find this dephasing rate to be constant, this 
means that the qubit behavior is well described by the relaxation measurement alone. 


6.5.1 Measurement scheme 


In order to get physically connected I, and IX rates, we acquire the measurement data 
for both measurements in turn, such that temporal fluctuations of the qubit properties 
influence both measurements likewise. This technique was developed for the measure- 
ments presented in Ref. and is explained in more detail there. Without doing so, 
the two rates do not need to be result from the same perturbations and measurements 
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Figure 6.10: Pure dephasing rate. The pure dephasing rate IT, can be seen as the y-intercept in this 
visualization. For the areas with a higher slope of wo (B) (dark crosses) we could expect a higher pure 
dephasing rate due to a higher sensitivity to magnetic field noise and fluctuations of the solenoid current. 
A clear trend is however not visible and all data points match well to the linear regression. The standard 
errors from the fits for T; and TR are shown as error bars. Adapted from [7o]. 


with TR > 2T; can occur, which makes no sense from a physical point of view. Before 
each measurement, the re-calibration routine as presented in Sec. is executed. 


6.5.2 Measured data 


We plot TÈ against the simultaneously acquired I rate, as shown in Fig. and 
fit the data to a straight line. From the y-intercept we find a pure dephasing rate 
of I; = 93.9 KHz. Due to the steeper slope of «o(B) for the regions close to the 
blockage of JJ, we would expect an increased I, for data points with lower wo. 
For the measurements presented in Fig. 6-10] where wg is color-coded, we however 
see no clear correlation between I, and wo. 


The causality between noise in the solenoid current 5; and the resulting I5; 


is given by [177]: 
dwo1 V? 
D 01 
Baal 3l ) Silw « n) (6.7) 


where the relevant scale for w is the time between the Ramsey pulses, being on the 
order of w/27t ~ 100 KHz. The measured power spectral density of our current source 
is about Sz x 10715 A? Hz 1, with which we can calculate T (B) as shown in the upper 
part of Fig. In the lower part, T, »(B) is color-coded onto the qubit transition 
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Figure 6.11: Calculated dephasing rate. From the measured noise power spectral density S; of our current 
source we can calculate the current noise contribution to the pure dephasing rate and plot it against the 
applied magnetic field (black curve). The gray dashed line represents the extracted pure dephasing rate 
from Fig. [6.10] To indicate the areas, where the qubit dephasing rate is limited by the current noise, the 


calculated I is color-coded into the qubit transition frequency «gi (B). This shows that only in the regions 


of a steep slope of «o (B) the current noise exceeds the other pure dephasing channels. 


frequency. From this plot, we can see that I; 9 is mainly not limited by the current noise 
but by other factors which are constant against the magnetic field. Only in the areas of 
a steep slope of «91 (B), the dephasing rate is increased by the magnetic field noise. 


6.6 Conclusion 


In this chapter, we measured the behavior of a superconducting transmon qubit in 
a magnetic field, where we demonstrated quantum coherence up to a flux density 
of B — 40mT. We explained the observed decay rates by different mechanisms 
and saw a strong indication that the loss mechanisms can be split into constant, 
hysteretic and non-hysteretic losses. Studying the pure dephasing rate, we found 
it to be constant and not affected by the magnetic field, such that the dephasing 
is governed by the qubit's relaxation. 


From the observation of the qubit's transition frequency, we found the influence of 
a second, parasitic junction. This junction is typically present in circuits fabricated 
by shadow angle evaporation, but commonly neglected as it has a high critical 
current. Due to its large cross section, it has an increased sensitivity to in-plane 
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magnetic fields. To calculate the influence of this parasitic junction, we derived 
a theoretical description which resulted in an approximated Hamiltonian for a 
transmon qubit with two serial junctions. 


Analog to the qubit, we studied the behavior of a superconducting resonator 
in a magnetic field. With our results being comparable to previously published 
data by other groups, we find a strong increase of the losses in the resonator, 
expressed by an increased linewidth. This broad linewidth finally appears as bottle 
neck for the measurement, as it impedes a fast and efficient readout of the qubit 
at even higher field strengths. 


From our measurements, and the partitioning of the qubit losses, we can derive multiple 
handles for increasing the qubit’s stability in a magnetic field. Most importantly, traps 
for moving flux vortices and excess quasiparticle could be implemented in future 
designs to decrease the influence of a magnetic field on the qubit coherence. Shunting 
or avoiding the parasitic junction is also strongly recommended, if a monotonous 
decrease of the qubit in the magnetic field is requested. 


Despite these possibilities for optimization, we find the non-optimized version of a 
state-of-the-art transmon qubit to be well compatible with magnetic fields beyond 
the critical field of bulk aluminum. Even though the decay and coherence times are 
strongly suppressed compared to the unperturbed case, they are still on the order of 
many other superconducting qubit designs, which have been used for outstanding 
publications. With our findings, the usability range of superconducting qubits as 
versatile sensors is increased and a promising basis for future developments in 
superconductor-magnet-hybrid devices is created. 
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Combining superconducting qubits with magnetic oscillations promises versatile 
applications in quantum sensing and quantum computation, as will be presented 
in this chapter. So far, we have studied in this work how we can perform quantum 
sensing with a transmon qubit and how a superconducting qubit behaves in a 
magnetic field. Following this route, our natural next step is the realization of a 
coupling between qubits and a magnetic material. 


In this chapter, I will start with a description of oscillating magnetic systems and 
provide basic theoretical background to understand the origin and the mechanisms 
of these oscillations. We then go into detail on the coupling mechanism between 
the qubit and this system, and show an experimental approach to realize a direct 
coupling. Finally, we present first results on this approach, showing a rich set of 
interesting features which are still subject to further research. This can be seen as 
a starting point for further experiments with superconducting quantum bits in a 
magnetic field, coupled to magnetically tunable systems. 


7.1 Introduction 


After the revolutionary development in the field of electronics in the past 75 years and 
the emerging number of applications in spintronics, a very active field of fundamental 
research is constituted by the field of magnonics. Here, the information is not stored 
and transmitted by the charge or spin of single electrons. Instead, collective excitations 
of the electronic spin system are used, where the long-range order of a ferro- or 
ferrimagnetic material supersedes the need for a physical transportation of electrons. 
These collective excitations can be seen as synchronous precessional motion of the 
spins around an externally applied field. In a quasi-particle description, they are called 
magnons and are described in greater detail in the next section. 


In the field of magnonics, many basic computation elements and logic gates have 
been realized and demonstrated, for example a magnonic transistor [52], magnonic 
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AND and XOR gates as well as directional coupler and half-adder [178]. Due to 
the decreased propagation velocity of magnons compared to electromagnetic waves, 
magnonic resonators are much shorter in length and therefore have a smaller footprint. 
This is especially beneficial for applications where a dense packaging is required. 
Since the spin wave can propagate without a particle transport, the loss mechanisms 
known from electronics do not apply. 


Although these devices have been actively developed and engineered for several 
years, little is known about the quantum behavior of magnons. Many experiments 
up to now have been performed at room temperature and high powers to create 
a non-equilibrium spin distribution. Magnonic properties at low excitations and 
fundamental studies about the included damping mechanisms however have only 
recently become a part of fundamental research. While cavity magnonics already 
enables to study the magnonic system at low temperatures and excitation powers, 
single magnon operations can only be performed when including a single photon 
device, i.e., a qubit. This creates the basis for experiments in quantum magnonics. 


Apart from the fundamental interest in exploring the single-excitation behavior of 
magnons and the underlying effects, one more applied goal is driving the research of 
qubit magnon coupling: Due to the interaction of magnons with optical photons 
[179], a bidirectional transducer between superconducting qubits and optical photons 
may be realized. This would be one of the missing building blocks for a quantum 
computer in a quantum internet, since superconducting qubits are among the most 
promising candidates to build the computation cells, and optical photons are almost 
exclusively used to transfer quantum information between different places. 


7.2 Magnons and magnetic oscillations 


In ferromagnetic materials, the elementary magnetic moments preferably align along 
a common axis due to an exchange interaction, which forms a long-range order. Since 
this interaction is stronger than the repulsion of their dipolar magnetic fields, all 
spins within a domain spontaneously align in parallel below a critical temperature 
[61]. For the following discussion, we assume that an external field is applied 
to create a uniform bias field throughout the sample such that all spins align in 
the same direction. We also assume that the temperature is low compared to the 
magnetic field, such that kgT < ji. B, where kg is the Boltzmann constant and ü 
is the magnetic moment of the individual spins. 


One intuitive excitation of the ground state would be to flip exactly one spin in the 
system, as illustrated in Fig.[7.1]a). This is however energetically not favorable, since 
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Figure 7.1: Excitations in a ferromagnet. Under the influence of an external field B, all spins align along 
the z axis in the ground state. a) One form of excitation is to flip one spin. This is energetically not 
favorable since the exchange energy between the flipped spin and its neighbors is needed. Energetically 
more favorable is a small tilt of all spins in the system, which then start a precessional motion around the 
external field. This can be in-phase for all spins as shown in b), or wave-like as shown in c). 


this spin would stand contrary to all its neighbors. An energetically more favorable 
solution is a distributed excitation, where all spins are slightly tilted, which can be 
seen as a linear superposition of all possible single-spin flips, illustrated in Fig.|7.1|b). 


If we compare the spin in the magnetic field to a spinning gyroscope in gravity, 
we see that a tilt immediately results in a precessional motion around the direction 
of the external force to keep the angular momentum constant. This precession is 
illustrated by the small curved arrows in Fig.[7.1]b) and c). Depending on parameters 
such as shape, material and the uniformity of the excitation, the precession can be 
in-phase as in Fig. |7.1|b) or can have a wave vector K with |k| = 271/A as illustrated 
in Fig. 7.1|c) for A = 1, where A is the wavelength. 


7.2.1 Magnon frequency 
Instead of describing the precession of each spin individually, typically the Holstein- 
Primakoff transformation [180] is applied, which describes all excitations in the 


ferromagnet by a superposition of spin wave excitations, or magnons. The Hamiltonian 
for a ferromagnetic spin wave then reads: 


(7.1) 


where the sum goes over all possible wave vectors k. In the general case and 
especially for structured magnetic materials, the magnon dispersion w(k) does not 
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only depend on the absolute value |k|, but also on the direction of propagation, 
the sample dimensions, and the orientation of the external magnetic field. Analog 
to other bosonic systems, fi. and f are the magnon creation and annihilation 


operators for the corresponding K vector. 


If we assume that all spins are precessing in phase, k= 0, we can replace the individual 
spins by a large macrospin. The precession frequency of this macrospin is given by 


Wp = [v8 


where y is the gyromagnetic ratio. We know that the flux density B inside the magnetic 
ae, is given by the external field Hoa and the internal magnetization M, ie., 
— Ho (Flex + M). This gives for the resonance frequency 


wp = YHo | Bexi + M| . 


In this formula however, the shape of the sample was not included. If we assume 
a thin film in parallel to the external field, we obtain [181 


wp = YHo V Hea (Hex E M). (7.2) 


This equation for the spin-wave frequency wp only holds for an in-plane orien- 
tation, does not cover the influence of a magnetocrystalline anisotropy field and 
assumes K ~ 0, i.e., a homogeneous AC excitation. We therefore take Eq. as 
reference for an order-of-magnitude estimation of the precession frequency and 
keep in mind that other modes exist. 


7.2.2 Hysteresis 


When we ramp down the external field, the spins still stay aligned with their neighbors, 
but due to dipolar interactions, a zero net magnetization would be energetically 
favorable. For this reason, domains form in the ferromagnet, which have different 
magnetic orientations and therefore reduce the overall magnetization, while still 
preserving a parallel spin orientation within the domain. In reality, the domain wall 
motion and reordering is dominated by material defects, leading to a magnetization 
depending on the magnetic history of the sample and therefore a hysteresis [61 .A 
sketch of the magnetization with respect to the applied magnetic field is given in 
Fig.[7.2] Three quantities are here of particular interest: The saturation magnetization 
Ms which is observed when all domains are aligned parallel to the external field, the 
remanence magnetization Mg which remains when setting the external field to zero, 
and the coercive field Hc, which has to be applied to reach zero magnetization again. 
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Figure 7.2: Magnetic hysteresis. The demagnetized sample starts at the origin and the magnetization 
increases with increasing applied field, until the saturation magnetization Ms is reached. When the 
magnetic field is decreased, the magnetization shows a hysteresis and only slowly decreases again. For 
Hext = 0, a remanent magnetization of Mr remains, and the magnetization only reaches zero, when the 
coercive field Hex: = Hc is applied. The shape of this curve was calculated by the Jiles-Atherton formula 
182] with symbolic values and does not represent the measured samples. 


For the experiments presented in Sec. we keep in mind that the magnetization 
changes quickly close to the coercive field when the field direction is changed and 
therefore the oscillation frequency is very sensitive to the magnetic field. Due to the 
formation and switching of different domains in this area, we do not expect a well 
defined mode spectrum. It is therefore beneficial to change the magnetization as little 
as possible by magnetizing the sample with a strong magnetic field and performing 
only low-range sweeps afterwards without changing the sign of the external field. 


7.3 Hybrid systems 


If we want to couple a magnonic system with a cavity or a qubit, we need to align 
the fields such that the oscillating AC field of cavity or qubit is perpendicular to 
the static magnetic bias field. Then, the AC field can excite the precessional motion 
of the spins, leading to a coupling of the two systems. This can be realized in 3D 
cavities or planar geometries, which is detailed in the following 


7.3.1 Cavity magnonics 
Using a macroscopic microwave cavity for the coupling with magnetic oscillations is 


advantageous for the first experiments: The spatial extent of a 3D microwave cavity 
assures a homogeneous AC field distribution over the magnetic sample in terms 
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of field strength and field direction. This eases theoretical considerations as well as 
numerical calculations, and allows to give a clear number of participating spins. The 
magnetic sample is typically of spherical shape, as this reduces the demagnetization 
field and creates fewer spatial magnon modes than rectangular shapes. These effects 
are beneficial when exciting the k=0 mode, the so called Kittel mode, and therefore 
help to create a simple and predictable mode spectrum. Also, the modular system, 
where cavity and magnetic medium can be easily exchanged, enables to study the 
exact same samples under different configurations. 


Most research in the field of cavity magnonics has hence been performed with magnetic 
spheres placed in the anti-node of the cavity's AC magnetic field [|62}/64| |183| 184]. 
These experiments opened new insights, but also gave rise to new research fields 
and raised new questions in the field of cavity magnonics. 


7.3.2 Cavity quantum magnonics 


To study these systems in more detail at low excitation powers, a natural approach 
is to introduce a nonlinear element, which allows to resolve the number of quantum 
excitations. This is realized by a superconducting qubit, placed in an antinode of the 
cavity’s electric field. Via a common coupling to a cavity mode, a virtual qubit-magnon 
coupling is established [67], and single magnon resolution was demonstrated [68]. This 
magnon counting can be equally regarded as the dispersive shift of the qubit frequency 
by the population of the magnonic system or by an AC Stark shift ofthe qubit due to the 
AC signal of the magnonic excitation. These studies render the first quantum resolved 
measurements on a magnonic system. In follow-up experiments, also a time resolved 
control and readout of the magnon-qubit interaction was demonstrated [69}/140]. 


Especially with regard to the inclusion of a superconducting qubit, the 3D approach 
simplifies several experimental issues: First, the range in applicable magnetic fields is 
increased, since copper cavities can be used. Fabricated with a sufficiently high quality 
factor, these cavities do not alter in a magnetic field, in contrast to a superconducting 
resonator, as found in Sec. 6.3] The increased field range also expands the range of 
applicable magnetic materials, as discussed in|7.4.1| Second, due to the spatial distance 
between qubit and magnetic material, the bias magnetic field can be generated locally 
at the magnetic sample. The qubit can be additionally shielded from the magnetic bias 
field by permalloy and superconducting shields. And third, the modular approach 
again enables to characterize the qubit-cavity and magnon-cavity system independently 
before combining them. In addition, individual components can be exchanged, if their 
frequencies have to be adjusted, while the other components stay unmodified. 
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Although this 3D approach is useful for first experiments, it has several drawbacks: 
Especially with respect to scaling and future applications, cavities have a large spatial 
extent and cannot be easily integrated. In addition, due to the narrow AC field 
distribution of the qubit, only a virtual, cavity mediated coupling between qubit and 
magnon can be achieved and a direct coupling is not possible in a 3D geometry. 


7.4 Direct qubit magnon coupling 


These drawbacks can be solved by a 2D approach, where a planar qubit-resonator 
system is coupled to a magnetic system. In contrast to the 3D approach, multiple 
systems can be placed on one chip with a reduced cross-talk and the overall dimensions 
for one qubit-magnon-cell are drastically decreased. In planar geometry, it is also 
more straight-forward to create a tunable qubit, which can be tuned by fast flux 
pulses and thus allows for a dynamical coupling. With the possibility to apply 6z 
pulses to the qubit, direct entanglement between qubit and magnon states can be 
reached and the coupled system can be brought into a superposition state. The planar 
approach also opens the possibility to directly couple the qubit to the magnetic 
material and avoid the detour via coupling cavity modes. 


For a direct coupling, the magnetic material has to be placed in the maximum of the 
qubit’s AC magnetic field. This corresponds to the areas with the highest current 
density in the qubit and hence the leads to the junction. The alternating current across 
the JJ creates an oscillating magnetic field, which can in turn excite the precession 
of the spins in a magnetic system on top of the qubit. 


7.4.1 Choice of magnetic material 


One of the most frequently used magnetic materials for coupling experiments in 
the microwave regime is Yid!] which is a ferrimagnetic insulator. Its very low 
losses and the therefore small linewidth are beneficial for resonance experiments, 
making YIG a very promising material for coherent information processing. For 
2D experiments, it was recently investigated in pre-characterization experiments 


with superconducting resonators 185]. 


A disadvantage for YIG is the need of GGQ?lss a substrate for the crystal growth process. 
This substrate shows strongly increased losses at low temperatures [186], which 


1 YIG: Yttrium Iron Garnet 
? GGG: Gadolinium Gallium Garnet 
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counteracts the low-loss properties of YIG. To circumvent this problem, YIG spheres 
are used in 3D applications, which are fabricated from bulk YIG. Alternatively, the 
substrate can be polished away after fabrication, which is technologically challenging 
and influences the mechanical stability of the sample. 


One additional and much stronger constraint is given by the field-dependent resonance 
frequency from Eq. (7.2). If we take literature values of y = 28GHz/T and pugMs = 
175 mT [187,188], we see that we need external magnetic fields on the order of Bext > 
100 mT to create resonance frequencies at wp = 5 GHz. This is clearly incompatible 
with the maximum magnetic field for a superconducting qubit, studied in Ch. |6| 


To decrease the bias magnetic field while still reaching resonance frequencies in the 
gigahertz range, we can increase the saturation magnetization according to Eq. 
by choosing a different magnetic material. The increased saturation magnetization 
is typically accompanied with an increased spin density ps, which increases in turn 
the coupling strength after Eq. (7.5). Permalloy, an iron-nickel alloy, is reported to 
have a high saturation magnetization and was recently investigated in combination 
with superconducting resonators [189] [190]. From the Slater-Pauling curve, which 
plots the atomic saturation magnetization versus the number of electrons per 
atom [191], we find an even higher maximum between the ferromagnets iron and 
cobalt. More detailed studies of Fe;Co;. , alloys find a maximum of the saturation 
magnetization for x ~ 30% [182]. 


In pre-characterization measurements [188], we measured a sample fabricated by 
an external group, where unfortunately the stoichiometry and fabrication methods 
are widely unknown. For this sample, we found a saturation magnetization of 
woMs = 2.1T for FeCo. Together with a gyromagnetic ratio of y = 28.2 GHz/T we 
find that a magnetic field of Be; = 15 mT suffices to bring the resonance frequencies 
to wp = 5 GHz. The measured data are shown in Fig. 7.3|together with a fit to Eq. (7.2). 


From these data, we can also extract the linewidth of the magnetic resonance, which 
corresponds to the damping of the spin wave oscillation. In the relevant frequency 
range of 4-8 GHz, we find on average km/27t = (350 + 19) MHz at room temperature. 
In the literature however, Fe,Co;. , films with an ultra-low magnetic linewidth are 
reported for x ~ 25%, resulting in km/27 = 56 MHz at wm/27 = 10 GHz [192]. 


108 


7.4 Direct qubit magnon coupling 


FMR frequency wrmr (GHz) 


—60 —40 —20 0 20 40 60 
applied magnetic field Bext (MT) 


Figure 7.3: FMR measurements on an externally fabricated sample. Results from a pre-characterization 
study on a FeCo sample. The orange line shows the theoretical curve for y = 28.2 GHz/T and ju9Ms = 2.1T 
together with an offset field of Bots = 3mT. The sample used for this measurement is a thin FeCo 
film, fabricated in an external lab with unknown stoichiometry and fabrication procedure, and not 
lithographically structured. Measurement data and fit values are taken from Ref. i 


7.4.2 Coupling strength calculation 
To calculate the coupling of a spin ensemble to the qubit, we start with the individual 


coupling strength g; of each spin. By applying perturbation theory to the Tavis- 
Cummings Hamiltonian |193], we obtain for the total coupling of the system 


N 
28i- (7.3) 


For the case of equal coupling strength, g; = g, this results in the well-known relation 
gs = gVN. In our case however, we assume that all spins S are equal, but feel a 
different local magnetic field B,(7;) at their location 7;. The index 1 emphasizes that 
the field is created by one excitation in the qubit. Under these circumstances, it holds 
gi = YSBı(F;), where y is the gyromagnetic ratio. 


We can calculate the field distribution B,(7;) with numerical simulations from the 
known sample geometry. Due to the limited precision of numeric simulations, it can 
be disadvantageous to simulate B; (7;) at the qubit current Ig, which corresponds to a 
single excitation in the qubit and is on the order of several nanoampere. Instead, we 
assume B « I, since the AC field can be considered as small perturbation to the DC 
bias field. We therefore use a higher current Isim for the simulation, which also makes 
the numeric simulation independent from the qubit parameters. We therefore obtain: 


P I N 4 
yim) = N In Jimatan (7.4) 
sim 


109 


7 Towards coupling of qubits and magnetic materials 


where we replaced the sum by an integral over the volume of the magnetic sample, 
Vs. The current flowing in a transmon qubit corresponding to its first excited 


state is found to be 117]: 


Iq = woe y hag /2Ec, 


where Ec again is the charge energy of the transmon. Introducing the spin den- 
sity og = N/Vs, we finally find 


woe /ħw 
Be / ps | IP rof avs. (7.5) 


Isim 


When m Eq. (7.5) to analytic expressions for the coupling strength found in 
the literature [62 ei. we see that these formulas involve quantities such as the mode 
volume Va = dV/|Bmax|* or a spatial overlap y instead. These quantities are 
useful for the " in where the cavity volume is confined and the magnetic field 
can be regarded as being constant over the magnetic material. For our 2D case, the 
integration boundaries for f |B|?dV are in principal infinity (in contrast to the confined 
space of dVs), which makes a numeric integration slow, and the field distribution 
at the magnetic material would not be respected. 


7.4.3 Simulation of the coupling strength 


The magnetic field distribution is simulated by 2D finite element simulation software, 
where we used FEMM?|for our purposes. For the simulations, we varied the spacing 
between the two chips, given by the resist pad thickness, and the width and thickness 
of the magnetic film. The results are shown in Fig. a) and b). We find that 
increasing the width of the magnetic structure further than w = 15pm does not 
significantly increase the coupling strength, as long as the chip spacing is small. In 
general, we try to reduce the size of the magnetic sample, as a dissipative metal 
close to the qubit gives rise to microwave losses. For the experiment however, a 
structure with w = 40 um is used, which is easier to align onto the qubit junction 
area. Together with a chip spacing of d = 1 um and a film thickness of h = 400nm, 
we calculate a coupling strength of g = 171.2 MHz. 


3 Finite Element Method Magnetics, version 4.2,|femm.info 
8 
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Figure 7.4: Simulated coupling between qubit and magnetic material. The coupling strength is calculated 
a) for different chip spacings d and b) for different film thicknesses h. A schematics of the simulated 
setup is given in c). The magnetic field distribution of the qubit is calculated and used for calculating the 
coupling strength after Eq. (75). As expected, the coupling increases with increasing width and thickness 
and therefore increasing number of spins, and decreases with increasing distance between the chips. For 
w Z 15pm, the coupling only slightly increases for a narrow chip spacing. Note that the chip spacing is 
measured between the substrates, i.e., for d = h + hg, the films would touch, where hg = 100 nm is the 
thickness of the qubit metallization. 


7.4.4 Flip-chip approach 


As seen from the simulations, we want to reduce the spacing between qubit and 
magnetic material to increase the coupling strength. A fabrication of both structures 
directly on top of each other would in principal be possible, but is accompanied by 
several technical challenges, i.e., damage of the Josephson junction in the deposition 
process of the magnetic material or contamination of the superconductor fabrication 
facilities with magnetic residues. Therefore, we chose to realize the coupling of the 
two systems with a flip-chip approach. Here, the two structures can be fabricated 
independently on individual chips and then placed face to face in a sample box. This 
design was already used for other coupling experiments with transmon qubits [102]. 
The sample box used there includes two rectangular chips which are rotated by 90? 
against each other. The distance between the chips is fixed by posts inside the sample 
box, which have a pre-defined step height and result in a minimum chip spacing of 
d = 35 um, due to the limited resolution of the CNC mill in our workshop. As seen 
in Fig. this distance is still too big for our application. 


The second approach was therefore to reduce the size of the top chip and take 
microstructured pads of optical resist as spacers between the chips. This avoids 
electrical contact as well as the creation of scratches in the structures in the alignment 
process of the chips. Since transparent sapphire is used as a substrate for the 
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Figure 7.5: Flip-chip approach. a) Photograph of the sample box with the qubit chip (black) in the 
center and the top chip (transparent) with the magnetic material. The box features two microwave 
connections, used for microwave readout and flux tuning of the qubit. The two chips are held in place by 
copper-beryllium clamps, which stay flexible at low temperatures. Design b) and micrograph c) of the top 
chip. A zoom-in of the design view is given in Fig. [7.6] The magnetic material is not centered on the chip 
area but shifted to the left, such that the top chip covers all of the underlying readout resonator (not 
shown). Otherwise, the resonator’s impedance would abruptly change in its center and additional modes 
would occur. In the dicing process, the substrate breaks off at the edges, which can be seen in c), but does 
not influence the functionality of the chip. 


upper chip, the magnetic structure can be aligned under an optical microscope 
to be above the junction area. The resist pads are structured by an optical mask 
and the resist is baked on a hotplate after development, making it durable against 
organic solvents. From typical fabrication parameters, we get a resist film thickness 
of d = 1.2um. The pads are positioned to be as far away as possible from both 
qubit and resonator to reduce additional microwave losses caused by the resist. 
The top chip is therefore only supported at its edges (see Fig. 75), leading to a 
bending in its center. This is increased by the clamps, which hold the top chip in 
place when mounting the sample in the solenoid. 


A manual adjustment of the chip positions with an accuracy in the range of few 
micrometers with the help of tweezers and a microscope is possible, as shown in 
Fig.[7.6]a). Although sapphire is a hard material, the chips touch already at a bending 
of only 0.05 % for a chip length of 4.4 mm, which can lead to scratches on all structures 
and eventually a galvanic contact. Great diligence is therefore needed in the alignment 
process to avoid damage as depicted in Fig. [7.6]o), where clear scratches across the 
qubit islands and the junction leads are visible. 
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Figure 7.6: Alignment of the top chip. a) Micrograph of the chip stack. Since the top chip is transparent, 
the alignment of the magnetic material on the qubit leads can be adjusted under a microscope. b) Design 
schematics of the stacked system. The qubit islands are depicted in blue, the junction leads in pink, and 
the magnetic material in red. c) Great care has to be taken when aligning the chips, since the structures 
can be damaged if too much force is applied. A deep scratch in the center island and the lower junction 
lead is visible after an unsuccessful alignment approach. The scale bar is valid for all subfigures. 


7.5 Measurements 


We installed the flip-chip sample in the cryostat and measured its properties at 
millikelvin temperatures under the influence of a static magnetic field. In the 
following, the measured data on the resonator and qubit frequencies are presented, 
followed by a discussion about the origins of the observed spectra. 


When we apply a magnetic field to tune the resonance frequency of the magnetic 
material, we have to keep in mind that FeCo shows a hysteretic behavior. We therefore 
expect more reproducible and reliable results when we do not sweep through the 
whole hysteresis loop, since domain wall shifting and domain switching could still 
occur in the frequency range of interest. Instead, we do not change the sign of 
the magnetic field for the following measurements and try to keep the change in 
magnetization as small as possible. If not noted differently, for the measurements 
presented in the following, we ramp the magnetic field to a value of B — 81mT 
for a short time, ramp the field down to 0 and cycle the base temperature above 
T — 5K to release all trapped flux in the superconductor. 


7.5.1 Resonator measurements 
When ramping up the magnetic field, we see several avoided level crossings in 


the resonator frequency. They show different coupling strengths, different slopes 
of the magnetic system's frequency with respect to the magnetic field, and even 
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Figure 7.7: Avoided level crossings in the resonator spectrum. The plot shows color-coded the amplitude 
of the reflected signal to find the resonance frequency of the resonator for different applied magnetic 
fields. The data are taken on an upsweep of the magnetic field. The insets show zoom-ins on different 
features we observed in the resonance spectrum. Unfortunately, the measured signal is too weak to extract 
a coupling strength for the different features by means of a fit, but we can estimate the coupling strength 
to be on the order of a few megahertz. 


different signs of the slope. One measurement is shown in Fig. 7.7Jand additional 
data can be found in Appendix|C| The coupling strength is just large enough to see 
the gap in the spectrum and the onsets of the anticrossing, but not large enough 
to track the resonance frequencies of the hybridized system. Therefore, and due to 
sudden jumps in the resonance spectrum, we are not able to fit the anticrossing and 
extract a coupling strength. Tracing the lines by eye, we can estimate the coupling 
strengths to be in the range of a few megahertz. 


Although the resonator is not directly coupled to the magnetic material, we would 
theoretically expect a weak coupling between both systems. This is on the one 
hand due to the large field distribution of microstrip resonators and on the other 
hand through the qubit-mediated coupling. The latter can be calculated by gpm ete = 
(gorgom)/ (wr — wo) [195], where gor (Som) is the coupling strength between qubit 
and resonator (magnetic material), and is expected to be in the range of a few megahertz. 


These data show a variety of magnetically tunable resonances, with coupling strengths 
in the same order of magnitude as expected from a spurious coupling between magnon 
and resonator. The variety of different modes would be explainable by taking into 
account the mode spectrum of the magnetic sample, as discussed in Sec. [7.5.3] 
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Figure 7.8: Four qubit sweeps in the magnetic field. Color-coded is the normalized dispersive readout 
signal in arbitrary units. White points represent masked data, where the resonator signal was not reliable. 
The order of the subplots corresponds to the order of measurement. It can be clearly seen that several 
features re-appear at comparable fields (2.3 — 2.4mT), but are different on each sweep. Other features, like 
the anticrossing at 2.58 mT, appear only in sweeps c) and d). 


7.5.2 Qubit measurements 


Scanning the qubit frequency, we find several avoided level crossings in its spectrum 
already at very low magnetic fields. Again, they differ strongly in coupling strength, 
slope and the sign of the slope. In total, they do not correspond to the magnetic 
resonance we expected. Due to the slow drifts of the individual resonances as well 
as sudden jumps, it is not possible to predict the exact position of an anticrossing 
for a new measurement. Together with a long data acquisition time, we are hence 
not able to resolve the anticrossings with a higher precision. 


The data of several qubit scans are shown in Fig. The described sequence of 
ramping the solenoid field to Be = 81 mT and afterwards cycling the fridge was 
executed for subfigures c) and d). 


A fit of the two unperturbed anticrossings in Fig. |7.8|c) and d) is given in Appendix|C] 
This is however subject to large errorbars, as we cannot see both branches of the 
avoided level crossing simultaneously due to the weak coupling strength. From the 
visible data we can make a rough estimate for the prominent avoided level crossing 
at Bext = 2.59mT and find a lower bound of g/27z 2 4MHz. 


This order of magnitude estimation can be enhanced in future measurement, which 
include a drive tone close to the qubit frequency. Now that the position of the avoided 
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level crossings is approximately known, the drive tone frequency can be scanned close 
to the qubit frequency. As described in Ch. 5|we then would expect a stronger AC 
Stark shift of the qubit levels when the drive matches the resonance frequency of the 
magnetic system, due to the resonant amplitude enhancement of the microwave signal. 
This would enable to trace the resonance lines for larger detunings and subsequently 
narrow down the estimate of the coupling strength. 


7.5.3 Discussion 


The acquired data on resonator and qubit frequency show a clear interaction between 
resonator, qubit, and multiple modes of a separate system, whose resonance frequencies 
can be tuned by a magnetic field. 


From theoretical and measured properties of the FeCo film, we would expect the 
magnetic resonance to reach the qubit frequency only at stronger magnetic fields 
on the order of Be = 20 mT. Possible explanations for this discrepancy are the 
existence of additional shape and magnetocrystalline anisotropy fields, which we 
ignored in the theoretical considerations in Sec. In addition, the hysteresis 
curve and especially the remanent magnetization are not known for the magnetic 
films used for the experiment. 


Depending on the orientation of the external field to the film, the shape anisotropy 
can also counteract the applied magnetic field [181], which can result in a de- 
creasing resonance frequency with increasing field. This behavior can be used 
to explain the anticrossings for which the second system's frequency goes down 
with increasing magnetic field. 


In our theoretical considerations, we only took into account the uniform mode 
with k = 0. Due to the inhomogeneous field created by the qubit leads, it would 
also be possible to excite magnon modes with k Æ 0, which in general show 
resonance frequencies depending on their wave vector, wm(k). This may be an 
explanation for the great variety of different modes we observe in our spectrum, 
each with their own coupling strength. 


Efforts to further verify these explanations, such as simple FMR measurements and 
a SQUID magnetometer did not give new insights, since the sample size and the 
number of spins was too small for the sensitivity of the measurement devices. These 
measurements could therefore neither confirm nor refute the theory that the avoided 
level crossings are hybrid resonances with the magnetic system. 


From Sec.B.6]we know that an avoided level crossing is only visible when all included 
linewidths are similar to or smaller than the coupling strength. In the measured 


116 


7.5 Measurements 


data however we see the onsets of avoided level crossings, but cannot trace the 
individual branches far enough to see both branches simultaneously. We therefore 
can follow that the largest linewidth is on the same order of magnitude as the 
coupling strength. With a lower boundary for g/2rt = 4MHz, as estimated from 
the qubit spectrum, together with the FMR linedwith previously measured on a 
different sample, km/27 = (350 + 19) MHz [188], we clearly would have to follow that 
the measured anticrossings cannot result from magnon oscillations. If we however 
take the value of 2/27 ~ 110 MHz which we get from the avoided crossing fit in 
Appendix |C| and assume a high-quality FeCo film with tm/27 < 56 MHz [192], 
the anticrossing would be well resolvable. 


Since the data basis is rather weak and the observed crossings do not match our intuitive 
expectation of an anticrossing with one single straight line, as expected from the FMR 
signal shown in Fig. we want to emphasize that this chapter makes no claim about 
the origin of the observed anticrossings. Instead, we can only state that we see a variety 
of resonances coupling to our qubit, which can be tuned by a magnetic field. It is well 
possible that these resonances do not have their origin in the magnetic material we 
placed on top of the qubit, but result from other effects, like a changing flux vortex 
configuration or, e.g., magnetically tunable defects. We want to note the following 
observations from our experiments and leave the conclusion up to the interested reader: 


e We could not observe anticrossings when the magnetic field was swept in a 
bipolar sweep (i.e., from positive to negative values or vice versa). The crossings 
only reappeared after applying a strong magnetic field in one direction. This 
supports the idea of having a hysteretic system. 


e This hysteresis "survived" a thermal cycle of the sample above the supercon- 
ducting transition temperature, and we found anticrossings at comparable field 
strengths repeatedly. This negates the idea that the resonances come from an 
interaction of the superconductor itself with the magnetic field. 


e In the pre-characterization of the qubit chip without a magnetic material on top, 
we did not see any anticrossing. There however, we did not perform the sequence 
of applying a strong magnetic field, cycling the base temperature and scanning 
the magnetic field with a fine resolution. The absence of avoided level crossings 
without the magnetic chip on top support the theory that the resonances arise 
from magnon oscillations. 
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7.6 Conclusion and outlook 


In this section, we presented the basic mechanism to achieve the coupling between a 
superconducting qubit and a magnetic material. We derived an analytical formula to 
calculate the coupling strength for a given geometry from the results of finite elements 
simulations. In the last part, we presented measured data on the coupled system, which 
show a rich set of features with different magnetically tunable resonances. The origin of 
these resonances could not be traced back unambiguously to the magnetic alloy on top 
of the qubit, but opens a new door to further analyze these resonantly coupled features. 


Based on the knowledge acquired from the measurements in this chapter, several simple 
experiments can be performed that could help to find the origin of the magnetically 
tunable resonances. These experiments include the measurement of a bare qubit 
in a magnetic field following the described measurement scheme, additional FMR 
measurements of the magnetic sample at low temperatures and the fabrication of new 
magnetic samples with a larger size, which could be characterized by other methods. 


If these additional measurements confirm the appearance of multiple magnon modes 
at the qubit frequency in an external field of a few millitesla, and also demonstrate 
the low linewidth of the selected FeCo alloy at low temperatures, the measurements 
presented in this chapter would render one of the first experimental demonstrations of 
a direct qubit-magnon coupling. On this system, more measurements like coherence 
analysis and noise spectroscopy could be performed to find out more about the magnon 
system. Replacing the qubit chip by a tunable version, together with a stabilized 
magnetic bias field, would enable to study the magnon spectrum on a larger frequency 
scale and find more information about the frequency-dependent coherence of the 
magnonic system. Eventually, a fast tunable qubit would enable to create entanglement 
between qubit and magnon states and bring the system in a superposition. 


For the case that the additional measurements demonstrate that the magnetically 
tunable resonances are independent of the magnetic sample on the qubit, the presented 
data would open the door to a new field of research. One likely option in this case 
would be that defects on the qubit, which have been studied in great detail for 
electric fields and strain, are also sensitive to magnetic fields. This effect can then 
be used to study these defects in more detail and learn more about one of the most 
prominent loss channels in superconducting qubits. 
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The goal of this thesis was to utilize a superconducting quantum bit for quantum 
sensing experiments and pave the way for the investigation of magnetic excitations in 
the quantum regime. By developing a new sensing scheme for oscillating microwave 
signals, studying a superconducting qubit in a magnetic field and interfacing the qubit 
with a ferromagnetic material, we created a solid foundation for future experiments, 
which allow to study magnons with a quantum sensor. 


For our implementation of a quantum sensor, I started this thesis with the development 
of a new sensing scheme, based on the AC Stark effect of an anharmonic multilevel 
quantum system [71]. By comparing the shift of higher excited qudit states with a 
model system, this scheme allows to characterize a microwave signal simultaneously 
in amplitude and frequency. Utilizing a fixed-frequency superconducting qubit, we 
demonstrated the sensitivity in a remarkable frequency range of more than one 
gigahertz and for a great range in signal strength. Outstanding in this sensing scheme 
is the locality of the measurement, which allows to measure signals arriving at the 
qudit position in situ. This is of unique advantage for the development, characterization 
and optimization of quantum circuits, as it enables to measure crosstalk and spurious 
couplings between different circuit elements. As demonstrated, our scheme allows 
to characterize the frequency-dependent microwave transmission function from the 
room temperature measurement equipment to the qubit at low temperatures. This 
calibration can be executed under the exact same ambient conditions as the actual 
experiments, giving valid calibration data. This includes temperature-dependent 
cable resonances as well as the frequency-dependent coupling between qubit and 
microwave transmission line, which is otherwise inaccessible. The presented sensor 
is also helpful for the development of hybrid quantum systems, where it allows for 
the characterization of excited resonance modes coupled to the qudit. 


The measurements with this sensor have been carried on in our group with a focus 
on increasing the sensor precision [72]. Utilizing time-resolved data acquisition, 
the change in qudit frequency was determined with higher accuracy, leading to an 
improved resolution of the sensor in both amplitude and frequency. 


Continuing the way towards enabling quantum sensing of magnetic excitations, we 
fathomed the sensor's limits in compatibility with a magnetic field. To this end, we 
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studied both frequency and coherence properties of a superconducting transmon 
qubit in an in-plane magnetic field. Although superconductivity suffers under this 
influence, we were able to demonstrate a quantum coherent behavior up to remarkably 
high magnetic field values of B ~ 40 mT, which was beyond the expected regime 
and is four times larger than the critical field of bulk aluminum [70]. These findings 
build the foundation for quantum sensing experiments with superconducting qubits 
in magnetic fields and shine a promising light on future applications. 


Analyzing our measured data, we found a second, spurious Josephson junction in 
the qubit to be relevant. This stray junction exists due to the fabrication process 
and is commonly neglected. Having a large cross section in the magnetic field, 
this additional junction gives rise to magnetically induced fluctuations of the qubit 
transition frequency. Although our calculations and measurements show that this 
spurious junction does not limit the dephasing time of our sample in the case of small 
applied fields, the situation changes for samples with larger junction cross sections and 
larger coherence times. Meanwhile, these spurious or parasitic junctions have gained 
additional research interest in independent studies, where they were found to host 
atomic defects, which increases decoherence independent of the magnetic field [37]. 


In the last experimental chapter, we combined the previous findings to build a 
hybrid chip, consisting of a superconducting quantum bit and a ferromagnetic 
structure. By analytical considerations and numeric simulations of the geometry, we 
estimated the coupling strength and performed pre-characterization measurements 
of similar magnetic materials. In the final experiment, where we directly interfaced 
qubit and magnetic material at millikelvin temperatures, we observed a rich set of 
interesting features when tuning the magnetic field. Several onsets of avoided level 
crossings between qubit and magnetically tunable resonances were visible, suggesting 
that only minor adjustments are needed to reach the strong coupling regime. The 
origin of this great variety in different magnetically tunable modes is still subject to 
investigation, but the presence of different magnetostatic modes is seen as a very 
plausible explanation. These promising data inspire the further investigation of 
qubit-magnon hybrid quantum systems and demonstrate the feasibility of a direct 
coupling approach in planar geometry. 


Outlook 


Quantum sensing with superconducting qubits has many near- and mid-term 
applications, especially with respect to improving quantum computation devices 
and developing new hybrid quantum systems. Already now, the sensing scheme 
presented in this thesis can be used to quantify microwave crosstalk in quantum chips 
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with multiple qubits, resonators and microwave feed lines. Based on the AC Stark 
shift, this sensing scheme can be applied without the need for pulsed microwave 
measurements or tunable qubit samples. Due to this astonishing simplicity, the 
transmon sensor is an ideal device to enable, e.g., quality assurance in the rapidly 
growing business field around cryogenic quantum technologies. Compared to a 
standard cryogenic characterization of microwave components using a vector network 
analyzer, only a transmon and two low-quality microwave sources are needed, 
without any additional cryogenic wiring. This would enable the growing number of 
suppliers for high-quality cryogenic microwave components, i.e., microwave cables, 
attenuators, circulators, and amplifiers, to provide a full cryogenic calibration at 
acceptable cost, which renders a unique selling point. 


Based on the broad field of current applications for SQUID sensors, which are 
used already now for the sensing of magnetic fields created by the human body 
in biomedical applications or the detection of cosmic radiation by magnetothermal 
effects, many sensing experiments are conceivable. The scope of these applications can 
be enhanced by replacing the SQUID with a superconducting qubit, and eventually 
introduce the sensitivity of quantum coherence and entanglement. Using methods 
like dynamical decoupling and spatial entanglement, the sensitivity to temporal 
and spatial fluctuations can be dynamically adjusted, opening new prospects on 
measurement accuracy and flexibility. 


The promising results we found in the experiments towards a coherent coupling 
of superconducting qubits and magnons allow for a positive outlook on future 
hybrid quantum systems. Enabling a coherent energy transfer between qubit and 
magnon system, the quantum mechanical dissipation mechanisms in magnons could 
be studied in-depth by relaxometry studies and phase coherence measurements, to 
find reasonable countermeasures. The entanglement of qubit and magnon states, 
and qubit-assisted relaxation time measurements for single magnons would pave the 
way for further investigations of magnonic systems with respect to applications in 
quantum storage devices or the realization of a quantum bus. In combination with 
the large variety of classical magnonic computation building blocks and systems 
like magnonic metamaterials, a versatile toolbox can be created to advance the 
implementation of quantum effects in computation. 


The great scientific interest in enabling a quantum-resolved study of magnon 
excitations can be seen by the number of experiments published during the course 
of this thesis. Studying magnons in a 3D cavity at low temperatures, we and 
other groups investigated the temperature and power dependence of internal 
losses in YIG, a commonly used ferrimagnet. A superconducting qubit was soon 
added to this system [66], but only very recently, a coherent control of the magnon 
states was achieved in this cavity-mediated 3D approach [69]. Following a planar 
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realization, combinations of planar superconducting resonators with ferromagnets 
have been studied at the same time [185} 194], and recently also other magnetic 
materials, which are more promising for a planar combination with a superconducting 
qubit, reached the focus of several research groups [189] [190]. 


Simultaneously, the magnetic field compatibility of superconducting qubits has gained 
more research attraction. Here, new materials are investigated, such as semiconducting 
nanowires [167] and disordered superconductors for the junction [196] as well ab fot 
the whole qubit [197]. Both approaches promise to increase the range in magnetic 
fields, in which superconducting qubits can be used and therefore lower the constraints 
for a direct coupling between qubit and different magnetic materials. 


Based on these developments and our promising findings, we are convinced that a 
coherent coupling between superconducting qubits and ferromagnetic excitations 
is possible in a planar geometry and will be realized in the near future. Finally, 
the results of these quantum sensing experiments contribute to the development of 
quantum computers and enable new communication interfaces. 
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Appendix 


A Perturbation theory 


The full expansion of Eq. up to fourth order in perturbation theory is given 
by the following formulas: 
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Appendix 


B Fabrication parameters of the magnetic sample 


Table 1: Magnetic material 


adhesion agent 
hotplate 110?C 50s 
HDMS in desiccator 24-5 min 


resist application 
resist AZ5214E 


ramp 500 rpm 5s 
spin 6000 rpm 60s 
hotplate 110°C 50s 
exposure 
dose 13 mW cm” 2s 
hotplate 110°C 50s 
flood exposure 13 mW cm? 30s 
developing 
AZ developer, H20, 1:1 33s 
sputtering 37 min 
Argon pressure 8.5 x 107? mbar 
Fe gun power 100W 
Co gun power 37W 
stripping 
NEP in ultrasonic bath level 1 5 min 
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B Fabrication parameters of the magnetic sample 


Table 2: Resist pads 


adhesion agent 


hotplate 110°C 50s 
HDMS in desiccator 2+5min 
resist application 
resist AZ5214E 
ramp 500 rpm 5s 
spin 6000rpm 60s 
hotplate 110°C 50s 
exposure 
dose 13mW cm~? 2s 
hotplate 110°C 50s 
flood exposure 13mW cm? 30s 
developing 
AZ developer, H20, 1:1 33s 
glassing 
hotplate 230 °C 6min 
protective coating 
resist AZ5214E 
ramp 500 rpm 5s 
spin 1000 rpm 60s 
hotplate 110°C 50s 
dicing 
stripping 
IPA ~20s 
NEP 1s 
IPA 
H20 
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C Additional measurement data on the coupled system 


Additional measurement data on a resonator scan are shown in Fig. 1| For the 
measurements presented in Fig. [7.8] c) and d) in the main part, we performed an 
anticrossing fit as included in our software package [141], shown in Fig. and Fig. B} 
The results are however only reasonable for the smaller splitting. 
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Figure 1: Avoided level crossings in the resonator spectrum. The amplitude of the reflected signal is 
shown as color against the VNA probe frequency and the applied magnetic field. Data was taken on 
an upsweep of the magnetic field. The insets show zoom-ins on different features we observed in the 
resonance spectrum. The experiment was performed similar to Fig. [7.7Jand shows a slightly different 
configuration of resonances. 
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C Additional measurement data on the coupled system 
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Figure 2: Avoided level crossing in the qubit spectrum. The measurement shown in Fig.[7.8]o) is fitted 
with an anticrossing fit. Although the fit results seem to accord with the data, the errors on the extracted 
values are larger than values themselves since the splitting of the curves is large enough that we cannot 
see both peaks at the same time. 
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Figure 3: Avoided level crossing in the qubit spectrum. The measurement shown in Fig. d) is 
fitted with an anticrossing fit. From the fit, we find values of g/27t = (3.30.3) MHz and a slope of 
(3.1 + 0.8) GHz/mT. 
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